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Abstract 

It is well known that free (unconstrained) initial data for the gravitational 
field in general relativity can be identified on an initial hypersurface consisting 
of two intersecting null hyper surf aces. Here the phase space of vacuum general 
relativity associated with such an initial data hypersurface is defined; a Poisson 
bracket is defined, via Peierls' prescription, on sufficiently regular functions on 
this phase space, called "observables" ; and a bracket on initial data is defined 
so that it reproduces the Peierls bracket between observables when these are 
expressed in terms of the initial data. The brackets between all elements of 
a free initial data set are calculated explicitly. The bracket on initial data 
presented here has all the characteristics of a Poisson bracket except that it 
does not satisfy the Jacobi relations (even though the brackets between the 
observables do). 

The initial data set used is closely related to that of Sachs |Sac62j . However, 
one significant difference is that it includes a "new" pair of degrees of freedom 
on the intersection of the two null hypersurfaces which are present but quite 
hidden in Sachs' formalism. As a step in the calculation an explicit expression 
for the symplectic 2-form in terms of these free initial data is obtained. 



1 Introduction 

Canonical formulations of general relativity are usually couched in terms of initial 
data on spacelike hypersurfaces [ADM62j |Sen82] [Ash86] . However it has been widely 
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known since the early 1960's that the initial value problem on piecewise null hyper- 
surfaces (hypersurfaces with null normal) is technically simpler than the spacelike 
initial value problem. The constraints reduce to ordinary differential equations, and 
free initial data, that parametrize the solutions to the constraints, can be identified 
explicitly [Sac62] [Pin63] |BBM62] [Dau63] . 

A canonical formalism in terms of the free null initial data would be of interest for 
many reasons. It might offer insight into the properties of classical solutions of general 
relativity, it should provide a powerful framework for the numerical computation of 
these solutions and finally, in the context of quantum gravity, it promises to shed light 
on black hole entropy and the entropy bounds proposed by Beckenstein and Bousso 
|Bec73j[5us95]|Bou99j . 

The main missing ingredient of such a formalism is the Poisson bracket on the 
space of free initial data. In 1978 Gambini and Restuccia |GR78j . using a judiciously 
chosen coordinate system, were able to obtain a perturbation series in Newton's 
gravitational constant for the bracket between the main free data. The main data is 
a specification of the conformal equivalence class of the (degenerate) 3-metric on the 
initial data hypersurface. In Gambini and Restuccia's work, and in the present work, 
this is determined by the "conformal 2-metric" , a field given on all of the initial data 
hypersurface, in contrast to the remaining data which lives on the intersection of the 
two null hypersurfaces only. 

Here I will present closed form expressions for the brackets between all elements of 
a complete free initial data set for vacuum general relativity. Although the approach 
taken is entirely different from that of Gambini and Restuccia, the bracket between the 
conformal 2-metric at distinct points obtained is equal to that obtained by summing 
their series. 

The first step in my approach is to define a Poisson bracket on a sufficiently 
rich family of sufficiently nice diffeomorphism invariant functions of the spacetime 
geometry of the Cauchy developments of the the initial data surface, which will be 
called "observables" . These Poisson brackets are defined by Peierls' prescription, 
which, unlike other definitions of the Poisson bracket, involves no arbitrary elements 
of the classical description of the system (such as the choice of time, or the choice of 
fields), and which is directly related to the quantum commutator. The next step is 
to define a bracket on a suitable free initial data set so that it reproduces the Poisson 
brackets between the observables. This requirement is fulfilled if the bracket is a 
certain generalized inverse to the symplectic 2-form. It turns out to be possible to 
calculate such a generalized inverse explicitly, yielding explicit closed form expressions 
for the brackets between all data. 

Unfortunately (at least from the point of view of a quantization project) the 
bracket obtained does not satisfy the Jacobi relations. The requirement that the 
bracket reproduce the Peierls brackets between observables seems absolutely neces- 
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sary. Nevertheless there still seems to be some freedom in the choice of bracket. 
Though the condition used here to determine the brackets on initial data is clearly 
sufficient to guarantee that the consequent brackets of observables reproduces the 
Peierls bracket, it is not clear that it is necessary. Some relaxation of this condition 
seems possible. Moreover, the brackets compatible with this condition are almost, 
but not quite, unique. The author is currently investigating whether a true Poisson 
bracket on initial data which reproduces the Peierls bracket on observables can be 
found. 

Null hypersurfaces are always swept out by congruences of null geodesies, and 
these generically form caustics. It seems preferable, at least at this stage, to work 
with smooth and structurally simple null hypersurfaces, so the hypersurfaces have 
been truncated such that they contain no caustics. This requires us to work with 
initial data on hypersurfaces with boundary, which brings with it new problems. 
On the other hand, it seems very natural and appropriate to general relativity to 
describe it via an "atlas of phase spaces"0 rather than to attempt to force it into a 
single phase space - thereby imposing global causality which, though comfortingly 
familiar, is entirely alien to the local nature of the theory. 

Be that as it may, the presence of boundaries has frustrated all my attempts to 
obtain the brackets by more conventional means, and is the reason that the bracket 
is a generalized inverse of the symplectic 2- form instead of a true inverse. Here a 
generalized inverse of the symplectic 2-form inverts the action of this 2-form on phase 
space vectors corresponding to perturbations of the gravitational field that vanish in 
a neighbourhood of the boundary of the initial data hypersurface. 

A variation of a solution can be represented as a variation of the initial data. 
However, for some initial data variables a variation of the solution metric that van- 
ishes near the boundary of the initial data hypersurface may produce variations of 
the initial data that do not vanish near the boundary. This occurs for instance for 
Sachs' initial data. It depends on how the initial data at point is related to the met- 
ric field. However, to calculate the generalized inverse of the symplectic 2-form it 
is convenient to use free initial data fields that are invariant near the boundary un- 
der variations that leave the solution metric undisturbed in a neighbourhood of the 
boundary, or indeed under variations that simply leave the 4-geometry undisturbed in 
such a neighbourhood. This requirement has led me to define a new free initial data 
set, differing somewhat from that of Sachs. It contains a new field on the intersection 
of the two null hypersurfaces. The information in this field is contained in Sachs's 
data but rather implicitly. The inclusion of this field is quite important in the phase 
space formalism since it is essentially the canonical conjugate to another datum on 
the intersection 2-surface (which is included explicitly in the Sachs data). 

1 I thank J-A Zapata for this evocative phrase. 
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In the following section I define the initial data hypersurface and the associated 
phase space, observables and the Peierls bracket between them, and the auxiliary 
pre-Poisson bracket on initial data. The relation of the Peierls bracket and the sym- 
plectic 2-form is found, and also the sufficient condition ensuring that the bracket 
on initial data reproduces the Peierls bracket on observables. The symplectic 2-form 
corresponding to the Einstein-Hilbert action is calculated explicitly in terms of the 
spacetime metric. Finally, the meaning of gauge invariance for initial data on a 
hypersurface with boundary is considered. In section [3] some convenient spacetime 
coordinate systems adapted to the geometry are presented and an initial data set 
suitable for our purposes is proposed and demonstrated to be free and complete. 
This section also contains a treatment of the relation between the variations of field 
components referred to different coordinate systems that are adapted to the fields. 
Then, in section H] the symplectic 2-form is expressed in terms of the initial data 
chosen. Finally, in section [5] the brackets between the free initial data are calculated 
in closed form. 

I close the paper with some concluding remarks about open issues and potential 
directions of research in section EJ 



2 Covariant canonical theory 
2.1 The problem 

General relativity is a local theory in the sense that initial data in a bounded region 
of space - a bounded achronal hypersurface in spacetime - suffices to determine com- 
pletely the gravitational field in a contiguous spacetime region, the causal domain of 
dependences of the hypersurface. It does not matter what is happening in the rest of 
the universe. We may thus focus attention on this domain of dependence and regard 
the space of values of the initial data on the hypersurface as a space of states, or 
phase space, of the gravitational field in the domain of dependence. 

Now in General Relativity diffeomorphic solutions to Einstein's field equations are 
equivalent, i.e. they represent the same predictions concerning measurements. Thus 
the phase space, which we shall denote $, actually consist of the equivalence classes 
of initial data sets corresponding to diffeomorphic solutions. Equivalently, it may be 
identified with the space of diffeomorphism equivalence classes of maximal Cauchy 
developments^] of data on the given hypersurfaceQ 

2 Achronal sets and domains of dependence are denned in appendix^ see also |HE73j or [Wald84j . 

3 The, slightly generalized, definition of Cauchy developments, used here as well as that of 
maximal Cauchy developments, is given in appendix [XJ For the sake of simplicity we shall only 
consider C°° developments. 

4 $ is often called the "physical" or "reduced" phase space to distinguish it from the larger 
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The aim of the present work is to find the Poisson bracket on the phase space 
$ corresponding to an initial data hypersurface consisting of two intersecting null 
hypersurfaces. 

What will actually be obtained is the Poisson bracket on a particularly nice subset 
of functions on this phase space, which in turn suffices to determine, almost uniquely, 
the brackets between initial data. The brackets between the initial data obtained 
define a pre-Poisson bracket on the algebra of phase space functions, that is, a bracket 
having all the properties of a Poisson bracket save that it does not necessarily satisfy 
the Jacobi relation. The lack of a Jacobi relation is of course an important issue when 
considering the quantization of the system, since commutators satisfy this relation 
identically. 

Let us define precisely the initial data hypersurface to be considered. It is a 



compact manifold (with boundary and edges) p A/", consisting of two null hypersurface 
branches Ml and Mr joined on a smooth spacelike 2-surface Sq diffeomorphic to a 
disk (See Fig. [[]). Ml and Mr are swept out by the two congruences of future directed 
null geodesies (called generators) emerging normally from Sq, and are truncated at 
spacelike disks - Sl and Sr respectively - before these generators form caustics. 

From the point of view of the initial value problem M is simply a compact 3- 
manifold with boundary, consisting of two 3-cylinders (Ml and Mr) joined end to 
end on the 2-surface Sq. The initial data consists of smooth 3-metrics specified on 
Ml and Mr, with a continuity condition at Sq, and some further quantities specified 
on So only (see subsection 13.21) . In order that the branches Ml and Mr be null in 
the solution matching the data, the 3-metric on each branch must be degenerate. 
That is, there must exist at each point a degeneracy vector, a vector tangent to 
the branch that is orthogonal to all tangent vectors, which corresponds to the null 
direction tangent to the embedding of the branch in spacetime. Moreover, the space 
of such degeneracy vectors must be one dimensional at any point at which the branch 
is embedded smoothly (which requires that the tangent vectors to the branch are 
mapped to a three dimensional subspace of the spacetime tangent space). In fact 
the metric on a cross section transverse to the degeneracy vectors must be positive 
definite. Requiring these conditions to hold at all points of Ml and Mr is equivalent 
to requiring that the generators form no caustics on M. See appendix [B] 

Now let us consider the maximal Cauchy developments of data on M. The rep- 
resentation of a initial data set, or phase point, by its Cauchy development will play 
a central role in the definition of Poisson brackets for the following reason: Recall 

"kinematic" phase space K, constructed without use of the field equations, (see for example jHT92j . 
A geometrical construction of K for asymptotically flat and for compact boundaryless spacetimes 
is given in [LW 90J). The Poisson bracket we are looking for - the Poisson bracket on the physical 
phase space <£>, expressed directly in terms coordinates on $ - is often called the Dirac bracket. 
5 Manifolds with boundary and edges are defined in appendix El 
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Figure 1: The upper diagram shows an initial data surface like M embedded in a 
2+1 dimensional spacetime. In 3+1 dimensional spacetimes each of the components 
■A/l, A/}?, and 5o, of course have one dimension more than shown in this diagram. 
This is shown in the lower diagram, which depicts a three dimensional hypersurface 
N corresponding to a 3 + 1 spacetime. Unlike the upper diagram, the lower diagram 
represents only the intrinsic differential topology of A/", and does not indicate it's 
embedding in spacetime. 
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that the bracket of two phase space functions is again a function on phase space, and 
therefore takes as its argument a phase point. It turns out that the value of the Pois- 
son bracket at a phase point can be very naturally defined in terms of the spacetime 
geometry of the maximal Cauchy development corresponding to that phase point. 
The fundamental definition involves only retarded and advanced Green's functions 
for the linearized Einstein equation on this Cauchy development. 

But how do we know that Cauchy developments of a given set of data on M 
actually exist? Rendall |Ren90] has proved for a particular set of initial data variables, 
equivalent to that we shall use, that if these data are smooth on TV they determine 
a C°° Cauchy development of at least a neighbourhood M C M of So in N '. In 
fact Rendall's proof, which reduces the initial value problem on Af to an equivalent 
initial value problem for data on a spacelike hypersurface, combined with standard 
results on the spacelike initial value problem ( |Wald84] Theorem 10.2.2), implies that 
the data on M determines a smooth maximal Cauchy development uniquely up to 
diffeomorphisms. 

I shall assume that Rendall's results in fact apply to all of Af, i.e. Af can be 
taken to be all of Af, assumption which I consider very likely to be correct. If it 
is not then the results are still valid if Af is replaced by a subset intj\f - which is 
entirely equivalent, by diffeomorphism invariance, to placing additional conditions on 
the initial data on Af. 

This assumption implies that there exists a smooth boundaryless, globally hyper- 
bolic0 metric manifold M containing Af and its Cauchy development isometrically in 
its interior0 

Indeed, in the present work a Cauchy development shall be defined as the domain 
of dependence of the initial data hypersurface embedded achronally in a manifold with 
a metric that satisfies Einstein's equation, at least within the domain of dependence 

6 Global hyperbolicity is defined in appendix O 

7 Rendall embeds Af in IR and then shows that there exists a smooth Lorentzian metric g' on 
IR 4 such that g' to the past of Af smoothly matches any Cauchy development (on IR 4 ) of the initial 
data to the future of Af. g' and a Cauchy development spliced together in this way defines a smooth 
metric on an open connected subset U of IR 4 excluding dAf but including intAf such that U — Af 
is not connected. Now the initial data used by Rendall is free. The only condition it must satisfy 
is smoothness. Thus it may be extended to a larger hypersurface Af + , of the same form as Af, 
but containing Af in its interior. There thus exists a boundaryless metric manifold U + containing 
all of Af which extends the Cauchy development of Af. The results of appendix [Bj which take as a 
starting point that Af is embedded achronally in a boundaryless Lorentzian metric manifold therefore 
apply. In that appendix it is shown that Af, and its maximal Cauchy development, are contained 
isometrically in a smooth, boundaryless metric manifold M which is globally hyperbolic (sec prop. 
lB~20l) . 

8 dX will always denote the boundary component of X viewed as a manifold with boundary. See 
appendix [XJ This boundary will in general differ from the topological boundary of the embedding 
of X in another manifold. 
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Figure 2: The domain of dependence of a pair of plane null rectangles in 2 + 1 
dimensional Minkowski space is shown. The straight lines ruling the rectangles are 
the generators. The domain of dependence is bounded to the future by the light cones 
of the two ends of So and the null surfaces orthogonal to Sr and Sl- 



(see def. lA.lip . Note that the domain of dependence D[M] of M is not in general 
a manifold, not even a manifold with boundary. On the other hand T> = intD[Af), 
the interior of D[Af], and intM together do form a manifold with boundary, 2). 'D 
with its metric is essentially the whole Cauchy development. The underline in V is 
meant to represent the fact that V consists of T> and its past boundary. (Prop. IB. Ill 
indicates that D[N] consists entirely of the future domain of dependence of M.) 'D 
may also be thought of as D[J\f] minus its future Cauchy horizon (prop. IB.22[) J^l 

By prop. IB. 231 the subsets V of all Cauchy developments, of distinct data, are 
diffeomorphic as manifolds (though of course not in general isometric). They are 
thus all diffeomorphic to a particular such manifold, T>q, which we may take to be 
the union of the interior of the hypersurface Ao — {t — \x\, \x\ <l,y 2 + z 2 < 1} and 
the interior of its domain of dependence in IR 4 equipped with the Minkowski metric 
ds 2 = -dt 2 + dx 2 + dy 2 + dz 2 - see Fig. 

[How can one be certain that A/", equipped with given initial data is indeed 
achronally embedded in the solution space. This question also arises in the spacelike 
Cauchy problem, but it becomes especially urgent when the initial data hypersur- 
face is piecewise null, like A/", because the generators will generally cross eventually, 
spoiling the achronality of the hypersurface they sweep out (by the corollary to the- 
orem 8.1.2. of |Wald84j ). Truncating the generators before a caustic is formed does 
not exclude the possibility that non-neighbouring generators intersect. (See Fig. [3]). 
Thus it seems one would have to truncate the generators prior to any crossing, and 

9 Cauchy horizons are defined in appendix [XI 
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Figure 3: An example of caustics and intersections of generators. Sq is a spacelike 
curve in 2 + 1 dimensional Minkowski spacetime having the shape of a half racetrack) 
- a semicircle extended at each end by a tangent straight line. The congruence of 
null geodesies normal to So and directed to the future and inward - the generators 
shown in the diagram - sweep out a null surface having the form of a ridge roof, 
terminated by a (half) cone over the semicircle. The generators from the semicircle 
form a caustic at the vertex of the cone. There neighbouring generators intersect. On 
the other hand generators from the two straight segments of So cross on a line (the 
ridge of the roof) starting at the caustic, but the generators that cross there are not 
neighbours at Sq. 
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the question arises how this requirement can be expressed as a restriction on the 
admissible initial data. 

It turns out that that is not necessary. Suppose M is swept out by the two future 
directed null geodesic congruences normal to a spacelike disk in a spacetime X, as per 
its definition, and that it contains no caustics of the generators. (Suitable restrictions 
on the initial data do assure that - see prop. IB. 41 ) Then if the generators of M 
cross, there always exists another spacetime, forming a locally isometric covering of a 
neighbourhood of M in X, in which the generators do not cross and the M is achronal. 
See appendix [Bj and in particular prop. IB. 61 Thus no restriction on the initial data 
is necessary to ensure the achronality of M. Moreover, the solutions in which M has 
self intersections and is not achronal are represented in the phase space of Cauchy 
developments by covering manifolds. 

A subtlety arises in connection with the Poisson brackets, which are defined in 
terms of small perturbations of the spacetime geometry. The small perturbations 
of the geometry of a non simply connected spacetime X are distinct from the small 
perturbations of its covering manifold X. Nevertheless it seems that if one is careful 
to identify functions of the geometry of X with functions of the the geometry of X 
that are invariant under the action of the fundamental group of X on X then the 
brackets between such functions calculated in X are equal to those calculated between 
the corresponding quantities in X.] 

All of the Cauchy developments can therefore be represented up to diffeomor- 
phisms by metrics on the single, fixed manifold V . Mapping the metric on V in each 
Cauchy development to T?o y i a the diffeomorphism of prop. IB. 231 one obtains a space 
So of solutions to Einsteins equations representing the Cauchy developments. The 
diffeomorphism equivalence classes of solutions in So form the phase space $. 

The fact that each metric in So is (diffeomorphic to) a maximal Cauchy develop- 
ment of null initial data on M is equivalent to the requirement that it be a solution 
such that intJ\fo C X>o is null in the solution metric and the interior of P is the in- 
terior of the domain of dependence of intMo in V , and also in any smooth extension 
of T>.o and its metric that satisfies the field equation. (The final condition ensures 
maximality. ) 

There is a further condition the metrics in So must satisfy. As we have seen, 
our assumptions regarding the initial value problem on M imply that the metric of 
the Cauchy development can be smoothly extended to a boundaryless manifold M 
containing A/", and in particular dj\f. This gives rise to a restriction on the metrics 
in So. If V U dJ\f in M is diffeomorphic to V Q U <9A/o in IR 4 then the restriction is 
simply that the solutions in <S must be smoothly extendible to dAfo in IR 4 . However 
we have not proved that these sets are diffeomorphic, only that V and V are, so the 
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form this restriction takes in So is not known with certainty at present E 

The requirement of smoothness at ON is relevant to the analysis of gauge trans- 
formations in subsection 12.31 In the definition of Poisson brackets, which is our main 
focus of interest, the issue can be sidestepped. The value of a Poisson bracket {A, B} 
at a given solution, or phase point, g involves only the variations of the quantities 
A and B under perturbations about the given fiducial solution. It is thus natural to 
work with metric fields living on the spacetime manifold M g of this fiducial solution, 
or on the domain of dependence of M in M g , rather than on Minkowski space or on 
the domain of dependence of A/"o in Minkowski space. S shall be the space of maximal 
Cauchy developments (or "solutions" for short) mapped diffeomorphically to V g , in 
complete analogy with S . T shall be the vector space of all smooth metric fields on 
M g . 

Of course, basing the description on 22 g and M g , rather than on V and Minkowski 
space, does not by itself resolve the issue of boundary conditions at ON for metrics 
other than the fiducial solution metric. One is just replacing one preferred solution, 
Minkowski space, by another, g. However, it does allow the issue to be resolved for 
the small perturbations relevant to the definition of the bracket. We shall see that the 
Poisson bracket {A, B} of real valued functions A and B on T involves the variations 
of A and B under perturbations of the metric in the space, L g , of smooth solutions 
on M g to the field equation linearized about gEE (see © or (HUD). 

In the sequel the only M and V considered will be those of the fiducial solution 
(whichever solution this is), so the subscript g will be dropped. 

It should be emphasized that while diffeomorphic metric fields may be considered 
physically equivalent they are not the same in our mathematical model. A manifold 
X consists of distinct a priori identifiable points (see the definition of manifolds in, 

10 What conceivably might happen, though I believe it cannot, is that on some generator from 
8Sq the future boundary of V in M is not only tangential to Af on the generator, as it is in flat 
spacetime, but it meets Af so "softly" that also the second derivatives (of spacetime coordinates) 
are equal on the two hypcrsurfaces. 

11 The linearized field equation has well defined solutions on all of M g even though the metric g 
need be a solution to the (vacuum) Einstein equation only on D[Af\. 

12 Note that if g is linearization stable then the space, L g \D_, of solutions to the linearized field 
equations in L g , restricted to V_, is identical with the space T g of tangents to the space of solutions 
to the full field equations on V_. At a solution that is not linearization stable L g \D_ is larger than T g . 

T g is essentially the tangent space to S at g £ S. The solutions in S are restricted by the 
requirement that the metric on 2? be a maximal Cauchy development of null data on Af, but in fact 
any variation of the solution about g £ S (i.e. any variation in T g ) can be made to respect these 
restrictions by adding suitable diffeomorphism generators, so the tangent space to S is T g with a 
certain "diffeomorphism gauge fixing" . 

It seems likely that most Cauchy developments of data on Af are linearization stable. Mon- 
crief |Mon75j has shown that Cauchy developments of compact spacelike Cauchy surfaces without 
boundaries are linearization stable iff they have no Killing vectors. 
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for example [Wald84] ) and it makes sense mathematically to speak of the value of 
the metric at the same point for different metric fields. This means that one may 
also distinguish adapted coordinate systems, such as Riemann normal coordinates, 
which depend on the metric placed on X, from fixed coordinates, which do not and 
are so to speak "painted on the points of X" . It also means that the variation Sg^lp) 
of the metric at a point p is meaningful, a fact which we shall use implicitly when 
formulating general relativity in terms of a variational principle and in the associated 
definitions of the symplectic structure and the Poisson bracket of the theory. 

Free initial data on M provide a convenient coordinate system on the phase space. 
By definition each valuation of such data corresponds to a solution - that is, the 
data are not subject to constraints - and this solution is uniquely determined, up to 
diffeomorphisms, on D[M). The algebra of functions on the phase space $ is thus 
represented by the algebra of functions of the initial data@ and the (pre-)Poisson 
bracket can be specified explicitly by giving the brackets between the initial data. 

Free initial data on piecewise null hypersurfaces have been found by several authors 
|Sac6 2j |Pen63) [BBM62] [Dau63] . Sachs [Sac62j in particular has proposed a set of free 
initial data suitable for M , and has argued (somewhat heuristically) that the values 
of these data determine the solution in a neighbourhood of M in P uniquely up to 
diffeomorphisms. Rendall [Ren90j has proved that these data, if C°° on each of the 
null branches Ml and Mr of M and continuous in a suitable sense at So, indeed 
determine a C°° solution uniquely up to diffeomorphisms in a neighbourhood of So in 
'D. This is somewhat weaker than what is needed, but it does support Sachs' original 
claim. Here, as already stated earlier in other words, it will be assumed that Sachs' 
data determine the solution, modulo diffeomorphisms, on all of 22- The brackets will 
be calculated between the elements of an initial data set, defined in subsection 13.21 
that is equivalent to, but slightly different from Sachs' data. 

To be exact, in the present work the pre- Poisson bracket on initial data will be 

13 Some distinct data sets define diffcomorphic solutions, so the initial data will be multivalued as 
functions on <f>. A function of initial data that represents a function on $ must of course have the 
same value on all data sets that correspond to a given diffcomorphism equivalence class of solutions, 
that is, on all equivalent initial data sets. When using the ADM spacelike initial data |ADM62j 
or the Sen-Ashtekar variables |Sen82| |Ash86j to represent functions on the physical phase space 
this requirement - the requirement that the functions be gauge invariant - is highly non-trivial, 
since finding the class of initial data sets equivalent to a given set requires solving Einstein's field 
equations. Here the freedom in the data corresponding to a given diffcomorphism equivalence class 
of solutions amounts only to the freedom to change coordinates on the two spacclike disks Sl and 
Sr, and in this way relabel the generators of Ml and Mr. The action of these transformations on 
the initial data is simple and explicitly known - it does not require the solution of the field equations 
- so the multivaluedncss of our initial data should not represent a problem. It is an interesting and 
open question whether this freedom in our data should be understood as a gauge freedom (generated 
by some sort of constraint?) or rather as a global symmetry. 
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evaluated at C°° solutions without Killing vectors such that the generators of a given 
branch, Ml or Mr, of M are either everywhere converging or everywhere diverging. 
The last condition ensures that the cross sectional area of a bundle of neighbour- 
ing generators is either monotonically increasing or monotonically decreasing, which 
makes it possible to use the so called area distance, a quantity proportional to the 
square root of this cross sectional area, to parametrize the generators. 

The corresponding limitation on the applicability of our formulae for the bracket 
on initial data is far weaker than one might think. Causality requires that the Poisson 
bracket between two observables is non-zero only if these depend on fields at causally 
related points (points connected by a causal curve), and it is natural (as will be 
discussed further on) to require the same for the brackets between initial data. On 
M causally related points must lie on the same generator (by prop. IB. Tj) . so the only 
non-trivial brackets between data on M are those between data on the same generator. 
Now suppose that the monotonicity conditions on the cross sectional area of bundles 
of generators does not hold for M, and suppose p G intSo is the origin of a generator 
on which we wish to evaluate a bracket. If the generators of both Ml and Mr have 
non-zero expansion rates at p, which generically is the case, then we can find a disk 
S' C Sq containing p in its interior on which these expansion rates are everywhere 
non-zero and of uniform sign. Our results then apply to a hypersurface M 1 swept out 
by the generators emerging from S' , for if the expansion rcitc £it p is negative then the 
Raychaudhuri equation guarantees that it remains negative along the generators until 
a caustic, or conjugate point, is reached, while if the expansion rate is positive the 
generators may be cut off before they begin to reconverge. Thus our results give the 
bracket between the data except when the expansion rate of the generators of Ml or 
Mr vanish at p, or when one of the initial data in the bracket pertains to a point which 
lies beyond the cutoff of M' because the neighbouring generators, initially diverging 
at Sq, have begun to reconverge. In fact, in the case of re converging generators the 
formulae we will find for the brackets of initial data are still meaningful and can be 
used to assign values to these. If the expansion rate of the generators vanishes at p the 
expressions obtained for the brackets of some of the data are still meaningfully but 
other brackets are singular. Whether any of these extrapolations are correct remains 
to be determined since these cases are not treated in the present derivation of the 
bracket. 

As the fundamental definition of the physical Poisson bracket an adaptation to 
general relativity of the Peierls bracket [Pei52j will be adopted. Let g be a solution 
to the field equations and let A and B be diffeomorphism invariant functionals of the 
metric, and suppose that the functional gradient 5A/5g^ v of A has compact support 
in spacetime. [The support of 8A/5g^ u will be called the domain of sensitivity of A, 
and will be denoted by sa] The value of the Peierls bracket {A, B} at g is determined 
by the first order perturbation of the solution, and thus of the value of B, caused by 



13 



the addition of a term proportional to A to the action, /. Of course the perturbation 
of the metric is not determined uniquely by the change in the action. Aside from 
the freedom to apply diffeomorphisms, one may choose a retarded perturbation, an 
advanced perturbation, or a combination of these. Let g + be a solution to the field 
equations corresponding to the modified action I\ — I + XA (with A a real constant) 
which equals g outside the causal future (domain of influence) J + [s A of sa and 
g~ a solution that agrees with g outside the causal past J~[sa] of saO (That is, 
g + is the solution with a retarded perturbation, while g~ incorporates an advanced 
perturbation.) Then the Peierls bracket is 

{A, B} = A A B = -^B(gl)\ x=0 - J^(<7a)U=o. (1) 
As the action for gravity we will use the Einstein-Hilbert action^ 

where Q is the spacetime domain of integration of the action, R the scalar curvature 
defined by the spacetime metric g^, g the determinant of this metric, and d A x the 
coordinate volume form in a coordinate system x fixed in the manifold Q. 

No boundary terms need be added to the action because boundary terms do not 
affect the Poisson bracket implicit in the action. (See [LW90] for a discussion of the 
effect of boundary terms on the presymplectic 2-form, and sub-subsection 12.2.11 for 
the relation between the presymplectic 2-form and the Peierls bracket.) 

This definition of the Poisson bracket has the virtue of being manifestly covariant 
and, moreover, is very directly related to the commutators of quantum theory. Peierls 
showed that for field theories in Minkowski space his bracket is equal to that computed 
from the conventional Poisson bracket on initial data. Further on I will show that this 
bracket has all the properties of a Poisson bracket, including antisymmetry and the 
Jacobi relation, which are not obvious from the expression ([T]). In fact it will be argued 
that the Peierls bracket between diffeomorphism invariant (and sufficiently smooth) 
functions of the spacetime geometry is equal to the conventional Poisson bracket 
between these functions implied by the Poisson brackets between the canonical ADM 
variables on a spacelike Cauchy surface |ADM62j . 



14 The causal past and future are denned in appendix 1X1 

15 The normalization of the action is chosen so that when a particle of mass m is coupled to 
the gravitational field with it's action normalized so that the canonical momentum of the particle 
is mdx^/dr, with r the proper time of the particle, then Newton's law of gravitation (with G = 
Newton's constant) is recovered in the appropriate limit. The sign conventions for the curvature 
tensor and scalar arc those of Wald84 , that is, R = R^ 1 ' with [V /I5 V,,]^ = R^ va p uj p for any 
1-form lo. 
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Unfortunately the initial data on M used in the present paper are too singular as 
functions of the spacetime geometry to allow an unambiguous, direct computation of 
the Peierls bracket between them. My approach is therefore to use Peierls' definition 
to define the bracket on a set of sufficiently regular diffeomorphism invariant functions 
of the metric on P, which will be called "observables" , and then to look for a bracket 
on the initial data which reproduces the Peierls bracket on the observables when these 
observables are expressed as functions of the initial data. 

I shall call a functional A of the metric on V an observable if it is diffeomorphism 
invariant and its functional gradient a^ u = "7=^^^ is C°° and of compact support 
contained in the interior of V. The support of a is of course the domain of sensitivity 
of A. 

If all measurable functions of the metric are observables in this sense then our 
approach is certainly justified because the only predictive content of a field theory 
is what it has to say about measurable quantities. The predictive content of a clas- 
sical field theory consists of the relations between measurable functions of the fields 
implied by the field equations. In addition if the field theory is derived from a local 
action, like general relativity is, then the action defines a Poisson bracket between the 
measurable functions of the classical fields which reflects the commutation relations of 
the corresponding quantum observablesl^l A classical field theory contains no further 
information about the classical world it models or the quantum world which might 
underlie this classical world. In particular the initial data and their brackets are no 
more than a means to represent the measurable quantities and their brackets. 

Thus, in a mathematical formulation of a field theory in terms of observables these 
should be nice enough, as functions of the field, that the Poisson bracket between them 
can be defined, and yet rich enough so that all measurable quantities are represented 
by observables, or at least the limits of sequences of observables. Furthermore it is at 
least convenient and probably necessary that the Poisson brackets of observables are 
themselves observables, so that the observables form a closed Poisson algebra. 

It is not clear that the class of observables we have defined above meet these 
criteria^ but it is reasonable to expect that our observables form at least a subset of 
an algebra of observables of general relativity that does, and, as we shall see, there 

16 The same classical theory can sometimes be derived from distinct actions that define distinct 
Poisson brackets, for instance in the case of the harmonic oscillator. This presumably means that 
there are also distinct quantum systems which behave according to this classical theory in the 
classical limit. The Poisson structure is merely a convenience at the classical level which allows 
one to put theories into a universal (canonical) form. It assumes physical significance only at the 
quantum level. 

17 It is not hard to see that the domain of sensitivity of the Peierls bracket of two of our observables 
is compact and contained in the interior of V_. However it is not obvious that the functional gradient 
of the bracket must be C°° in spacetime. 
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are enough of them so that their Peierls brackets, together with a reasonable causality 
condition, determine almost unique brackets between our null initial data. 

Do observables as defined above in fact exist? An example shows that they do: 
For a generic C°° metric g on T), and a point p in the interior of P, C°° coordinates 
consisting of scalars formed from the curvature can be found which are non-singular 
on a neighbourhood of p. Choose a compact domain cq C V, with non-empty interior, 
in this neighbourhood and let [c] cE be the coordinate range corresponding to cq. 
Then for all metrics in a suitable neighbourhood U C T of g, in which all components 
of the metric with respect to the curvature scalar coordinates are nearly equal to those 
of g throughout the coordinate range [c], the same curvature scalars continue to be 
good spacetime coordinates on the pre-image c of [c] and this pre-image will be a 
compact subset of V diffeomorphic to cq. 

Now note that the metric components taken in these coordinates are diffeomor- 
phism invariants. Suitable functionals of these metric components on [c] define ob- 
servables. The domain of sensitivity of such a functional is necessarily a closed subset 
of c, and thus a compact subset of the interior of V. [The functional gradient of such 
a functional, A, vanishes outside c because a smooth infinitesimal variation of the 
metric which vanishes on c leaves the curvature scalars invariant on c and changes 
them only smoothly and infinitesimally outside c. This implies that A is unchanged 
because the preimage of [c] remains unchanged (i.e. equal to c) and the metric com- 
ponents with respect to the curvature scalar coordinates are unchanged in [c].} It is 
not difficult to construct explicitly a functional A of this form so that its functional 
gradient exists, for metrics in U, and is C°° in spacetime. For instance, if x^ are 
manifold fixed coordinates, z a are curvature scalar coordinates, and / is a smooth 
function from H 4 to IR that vanishes outside [c] then 

A = f fXz)^-det[g af} \(z)d A z= f f(z(x))yj ' -det[ 9lxv \{x)d A x (3) 

J[c] Jc 

satisfies all requirements. (It is less obvious how to construct observables that are well 
defined on all of JF, or S, and still satisfy our requirements, but that is not necessary 
here.) 

For non-generic metrics g which have Killing vectors the situation seems to be 
different. The construction of functionally differentiable observables which distinguish 
between metrics with isometries (Killing vectors) seems to be problematic. Perhaps 
this should not surprise us since the quotient manifold of the space of metric fields T 
by the group of diffeomorphisms of V is singular at metrics with isometries. Notice 
also that observables that are actually measured are always defined using some sort of 
inhomogeneity or asymmetry, if not in the metric then in the matter fields, to single 
out points or regions of spacetime. Recall that metrics with isometries are explicitly 
excluded from the scope of the present work. 
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The definition (CD) of the Peierls bracket has one serious problem: It is in general 
not obvious that the exact stationary points g + and g~ of the modified action I + \A 
actually exist. Fortunately this problem is easily overcome. The exact stationary 
points g + and g~ are not really necessary in order to define the Peierls bracket. The 
first order (in A) approximations to these are sufficient. The Peierls bracket is a 
creature of first order perturbation theory and is well defined whenever this theory 
is: Exact stationarity of the modified action requires 

+ 16ttGAcv = (4) 

where G^ = — l/2g )lu R is the Einstein tensor. This is just Einstein's field 
equation with stress energy tensor — 2Aa MJ/ . [The diffeomorphism invariance of A 
implies that is divergenceless, as the stress energy tensor in Einstein's equation 
must be: For any C°° vector field \ 

= £ x A = -2^ cvVW^ d 4 x = 2 J^W^a^V^Q d A x (5) 

and thus Vtv = oQ 

On the other hand, if g is a solution to Einstein's vacuum equation G^ = 0, then 
the modified action is stationary to first order in A at g + A5^g provided 

SaG^ + 16-irGa^ = 0. (6) 

Here 5a acting on a functional of the metric is the linear order variation of the 
functional due to the variation 5a9 of the metric. 

The Peierls bracket of two observables A and B may be defined by 

{A,B} = A A B (7) 

with 

*A = 5t- 5a (8) 

where 5\g and S^g are solutions to (jSJ) that vanish outside J + [s^] and respec- 
tively. 

When the exact solutions g + and g~ exist then dg ± /d\ satisfies the conditions 
defining 5^g, so the definition (pQ) of {A, B} is consistent with (|7j) and (JSJ) - The 
advantage of the definition via (ED) and (jHJ) is that one can show that 5\g and 5^g 
always exist, provided A is an observable (a is C°°, divergenceless, and of compact 
support contained in the interior of 22), and is unique up to diffeomorphisms, assuring 

18 Conversely, if A is defined on a diffeomorphism invariant domain of metric fields in T and a is 
divergenceless on all of this domain. A is diffeomorphism invariant. 
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the existence and uniqueness of {A, B} for any observable B. All that is needed 
are Green's functions for linearized general relativity, so that the advanced (— ) and 
retarded (+) perturbations of the metric S^g can be expressed aS 

Sig ar (y) = -16nG J a^(x)G% aT (x,y)V^9 d 4 x, (9) 

where G~ and G + are respectively the advanced and the retarded Green's functions 
for the linearized Einstein equation in the gauge used. In terms of these Green's 
functions the Peierls bracket is 

{AB} = 16vrG J a^(x)A^ ar (x,y)f3° T (y)V=gd 4 x y^d*y, (10) 

with A = G + — G~ and (3^ = ^7=^7 the functional gradient of B. Note that (fTUj) 
not only provides a quite explicit definition of the Peierls bracket, but also allows the 
bracket to be extended to a much wider class of functionals than our observables. 

To define the desired Green's function it is necessary to fix the freedom to change 
solutions to dHJ) by linearized diffeomorphisms, ^aQ^v ^aQ^v + 2V( At £ I/ ). We shall 
adopt the "transverse gauge" in which it is required that the trace reverse, 7^ = 
Ifiv — \9fj,u(g crp 1a P ), of the perturbation 7^ = 5g^ of the metric satisfies the condition 

X u = vv = o- (11) 

The linearized Einstein tensor may be written as 

SG^ = S7 (M Xu) - l/2g^ aX a ~ 1/2(V CT V CT 7^ " 2R% U T % T ), (12) 

so if ( TlTT) holds for &Ag then ([6]) is equivalent to the linear, diagonal, second order 
hyperbolic system 

V CT V CT V - 2i?%/7<x P = 32ttGV (13) 

(see |Wald84] section 7.5.) This equation defines unique advanced and retarded 
Green's functions |GDD82j . with which one may calculate the unique advanced and 
retarded solutions, 7 ± ; to ( fTBl) . Moreover, these Green's functions are regular in the 
sense that the solutions obtained by integrating them against a C°° source distribution 
a are also C°°. 

To show that the 7 correspond to solutions of the linearized Einstein equations 
([H]) it is sufficient to show that they satisfy the transverse gauge condition. Taking the 

19 As already noted — 2Xa acts as a stress energy tensor generating the perturbation of the metric. 
The first order (in A) expression given in © for the retarded perturbation should be compared 
with formulae for the gravitational radiation produced by a given stress energy, (see for example 
|Wald84] 1. 
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divergence of (fl3l) and using the Bianchi identity and the fact that the unperturbed 
metric g^ v satisfies Einstein's vacuum field equation one obtains 

V'V.VV = 327rGV"ev, (14) 

a linear, diagonal, second order hyperbolic system for V M 7 M ^, with the divergence of a 
as the source. Since this source term vanishes, because A is diffeomorphism invariant, 
and 7 =b vanishes either to the past or to the future of the support of a, V M 7^, = 
everywhere. The corresponding metric perturbations, S^g^ = yt v — l/2g fjLV g aT ^ T , 
which may be expressed in the form ([9]), thus are transverse gauge solutions to the 
linearized Einstein equations with source a, ([6]). They may therefore be used to 
evaluate the variations 5^B of B. 

Any variation of the metric may be transformed to one satisfying the transverse 
gauge condition by a suitable linearized diffeomorphism - see |Wald84j . The unicity 
of the advanced and retarded solutions to (1131) therefore implies that all retarded and 
advanced perturbations satisfying ([6]) are equal, modulo a linearized diffeomorphism, 
to S^g^u and S^g^u respectively. Since B is diffeomorphism invariant, they all define 
the same variations of B, 5^B, as these. {A, B} is thus well defined and unique. 

The difficulty in applying Peierls' definition of the bracket directly to null initial 
data can now be understood. The basic problem is that the Green's functions G (p, q) 
and A(p, q) are discontinuous in q on the light cone of p (the boundary of the causal 
domain of influence J + [p] U J~[p] of p in M). On a flat spacetime G + and G~ 
are distributions supported entirely on the light cone itself. On curved spacetimes 
additional contributions should appear that are smooth functions on the interior of 
J + [p] but discontinuous on the boundary@ 

This means that the Peierls bracket ( fTUl) between null initial data is ambiguous: As 
will be evident from the definitions of our initial data, the domain of sensitivity of an 
initial datum associated with a point p £ J\f consists of (a portion of) the generator, 
or generators, through p. The Peierls brackets between initial data at causally related 
points of A/", which necessarily lie on the same generator (by prop. IB. 71 in appendix 
IB1 . is thus ambiguous because the integral ( ITUl) is ambiguous in this case. It would 
be interesting to find out whether the bracket on null initial data found in the present 
work can be obtained from a simple disambiguation of (FlUj) . 



2 The retarded Green's function in curved spacetime for the masslcss Klein-Gordon equation, 
which is very similar to the transverse gauge linearized Einstein equation, can be estimated in the 
limit in which the two points, p and q, are close using an adiabatic expansion (sec Birrel and Davies 
|BD82j p. 74). In normal coordinates about p the zero order term is the Minkowski spacetime 
Green's function, and the lowest order correction (in a vacuum Einstein spacetime) is 

j^R xp \R^(q - pT(q - p) v o(p, <?), (is) 

where 6(j>, q) is +1 when q is in the causal future of p, and zero otherwise. 
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2.2 The Peierls bracket, the presymplectic 2-form, and the 
bracket on initial data 



We have defined the Peierls bracket between observables, our aim is now to find a 
bracket {•,•}* on the null initial data which reproduces the Peierls bracket on observ- 
ables when these are expressed in terms of this data. In the following I will show 
that the definition (0 [H]) of the Peierls bracket in terms of first order perturbations 
governed by equation © implies that it is, in a certain sense, inverse to the presym- 
plectic 2-form on the space L g of solutions to the linearized field equations, an object 
which can be computed explicitly in terms of initial data from the action. We shall 
see, moreover, that to ensure that the bracket on the initial data on M reproduces 
the Peierls bracket on observables it is sufficient that this bracket is also inverse (in 
a slightly different sense) to the same presymplectic 2-form. This condition will be 
used in section |5] to calculate the bracket on initial data. 

2.2.1 The Peierls bracket and the presymplectic 2-form 

The relation between the Peierls bracket and the presymplectic 2-form is valid for 
any local, Lagrangean field theory - that is, for any theory derived from an action 
which is the spacetime integral of a Lagrangean density which is a function, at each 
point, only of the fields and their derivatives up to a finite order (see |LW90j for a 
precise definition of this class of theories). The specific features of general relativity 
are therefore not needed in the demonstration of this relation, and a notation which 
pushes these features into the background will be adopted. The key result, equation 
flHHj) . may also obtained less abstractly working with the explicit form of the linearized 
Einstein equations, as is shown in appendix where (IHU]) is obtained from a Kirchhoff 
type formula for the solution to the linearized field equations in terms of initial data. 

Our first step will be to demonstrate a relation between variations of an observable 
A and the perturbation A a produced by that observable. Specifically it will be shown 
that the contraction of A J 21 l with the presymplectic 2-form is equal to the gradient 
of A when both are pulled back to the space L g of perturbations about a solution g 
which satisfy the linearized field equations!^! 

We shall work with the action on a domain Q bounded to the past by M and to the 
future by a smooth closed Cauchy surface E J^j of T>, the interior of which is contained 

21 Like any smooth variation of the metric field, can be regarded as a tangent vector to T . 

22 What this says is essentially that is (minus) the Hamiltonian vector field, va generated by 
the Hamiltonian A and the presymplectic 2-form Q. In general va is defined by the condition that 
5A = —Cl[vA,S] V<5 tangent to the phase space. 

23 Cauchy surfaces are defined in appendixEl The existence of Cauchy surfaces of 2? is guaranteed 
by the global hypcrbolicity of T> (prop. IA.4I or |HE73j prop. 6.6.3.) and Geroch's theorem, prop. 
6.6.8. of [HE 73] . That these may be taken to be smooth follows from the result of |BS03j . 
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in V. Specifically Q is the closure in M of the open set I + [Af] PI 7 _ [£ + ] C S + is 
chosen to lie to the future of the domain of sensitivity, sa, of A so that sa is contained 
in the interior of Q. (See Fig. 12.2.11 ) 

By prop. IB. 151 S + is compact, so prop. 6.6.1 of |H E73j and its corollary imply 
that, in the globally hyperbolic spacetime M, J + [Af] H J~[E+] = /+ [Af] H I [£ + ] D Q 
is compact. Since Q is closed it is also compact. Integrals over Q are therefore defined 
and the Stokes theorem applies to Q and dQ. 

Af and £ + form the entire boundary of Q in M: Suppose p G Q, then, since 
Q is the closure of a subset of V, p e P. Since V = intD[Af] = m£Z?[E + ] |HE73j 
prop. 6.5.1 requires that every inextendible timelike curve 7 through p crosses both 
Af and £ + . 7 must cross TV at p or at a point r to the past of p, because Q C / + [A/] 
which is disjoint from I~ [Af] since is achronal, and similarly 7 must cross S + at 
p or at a point s to the future of p. Thus if p does not lie in either Af or S + then 
j> G J + (r) n/~(s), and therefore p lies in the interior of Q. (This is of course consistent 
with prop. IB. 151 which shows that <9£ + = dAf.) 

For later convenience S + and Af will both be oriented with the positive side to the 
future@ Thus, as oriented manifolds, dQ = S + U —Af and <9£ + = dAf. 

At an exact stationary point of the modified action I + XA the variation of this 
action is a boundary term] 26 ! 

SI + \6A = (j)[5}, (16) 

where is the integral of a local functional of the fields and their variations over the 
boundary of Q @ 

24 The chronological past, i" - [S], and future, I + [S], of a subset S of spacetime are defined in 
appendix El 

25 The space of n-forms at a point of an n dimensional manifold is one dimensional. Thus there 
are two equivalence classes of such n-forms under multiplication by positive real numbers. The 
orientation of a manifold is a continuous choice of one of these equivalence classes as the "positive" 
one. A vector transverse to an n — 1 dimensional submanifold emerges on the "positive side" of the 
submanifold if the contraction of the vector with a positive n-form is a positive (n — l)-form of the 
submanifold. 

26 By a stationary point of an action we shall always mean a field configuration at which the action 
is stationary with respect to all variations that vanish in a neighbourhood of the boundary. That is, 
we mean a solution to the field equations corresponding to the action. 

27 We will always take the view that 8 represents the derivative along an unspecified one parameter 
family of field configurations, and is thus a vector in the tangent space to the space of fields being 
considered. Another point of view, adopted in [CW87] . is that S is the exterior derivative operator 
on the space of fields, and it's contraction with a vector is the derivative along the vector. Though 
this second point of view is not adopted in the discussion (and we denote the exterior derivative 
operator on the space of fields by dl) our formulae, with slight modifications, may be interpreted 
in this way, if one takes S without subscripts (and dl) to be the exterior derivative, and 6s with 
subscripts to be directional derivatives along vector fields. 
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Figure 4: The diagram shows schematically (or literally in 1 + 1 dimensional space- 
time) the initial data hypersurface J\f, its maximal Cauchy development (minus future 
Cauchy horizon) P, a boundaryless manifold M containing V and A/", a Cauchy sur- 
face S + of the interior of V, the domain Q between M and S + , and the domain of 
sensitivity sa of an observable A. The support in V of the retarded and advanced 
perturbations of the metric generated by A are also indicated - as white regions where 
the vertical or diagonal hatching is absent. 
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As we have seen already, whether or not exact stationary points of / + XA exist, 
approximate stationary points where the action / + XA is stationary up to first order 
in A can be constructed by adding to any exact stationary metric g of the unperturbed 
action / a suitable correction Ao-^g linear in A: At the perturbed metric g + \5a9, 

5I + X5A = <j)[5]+0{X 1+ ), (17) 

where 0(A 1+ ) vanishes faster than A as A — ► 0. 

Note that A does not contribute to 0, because its domain of sensitivity sa does 
not intersect the boundary of Q, and consequently is independent of A. When I is 
the Einstein-Hilbert action (j2J) on Q: 

= -idL^"^^ (19) 

where V is the metric compatible covariant derivative, V is the Christoffel symbol, 
and gE^ =1/3! e^ pT dx v A dx p A dx T is the coordinate volume 3- form on <9Qo 

The linearized field equation that governs &Ag is obtained by differentiating ( flTl) 
with respect to A and setting A = 0. Note that the derivative receives contributions 
both from the explicit A dependence in (TTTI) and from the A dependence of the metric 
g + Ao-^g. The result is 

5 A 5I + 5A = 5 A (f>[5], (22) 

where, as in ([6]), 5a acting on a functional of the metric is the linear variation of the 
functional due to the variation 5a9 in the metrical The contribution from the term 
0(A 1+ ) vanishes because dO(X 1+ )/dX\ x =o = lim A ^ 0(A 1+ )/A = 0. 



28 The exterior, or wedge, product of a sequence of forms, of valence n\, n,2, .... adding up to a 
total valence of n, is defined to be n\/(n\\ri2\...) times the antisymmetric component of their tensor 
product. For example if a is a 1-form and b a 2-form, then 

[aAb] lu , <7 =3a,[ IJ b 1 , a .]. (20) 

The integral of an n-form a over an rt-manifold X is defined by 

1 

ix n - Jx 



a= n\ j^' ^^^ -^^x, (21) 



with e antisymmetric and g 12 - " = l. 

29 Note that 8 and 8a in (|2"2"|) may be viewed as vector fields on the space Tq of field configurations 
on Q, but they may equally well be taken to be vector fields on the space T of field configurations 
on all of M since no term in (|2"2")l depends on their action on the fields outside Q. We shall take the 
latter point of view. 
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In general relativity the validity of ( l22l) for all 5 with support contained in the 
interior of Q is equivalent to the validity of (J6j) there. Of course (1221) also has a 
boundary term, probed by 5s supported on the boundary. This term provides ad- 
ditional information, not contained in which is the focus of our present interest 
because, as we shall see, it determines the symplectic structure and Peierls bracket 
of the theory. 

Perturbations g + r5og of g which are first order (in r) stationary points of the 
unperturbed action / obey the same equation, but without the source term SA: 

S SI = 5 Q <f>[5}- (23) 

In general relativity the bulk term of this equation corresponds to the linearized 
vacuum Einstein equation on the interior of Q. 

Setting 5 = 5$ in (1221) and 5 = 5a in ( 1231) one obtains 



5 A = 6 A <f>[6 ]-6 A 6 I (24) 
= 6 A <i>[6 ]-6o<j>[6 A ]-[6 A ,6 ]I (25) 

The commutator [5 A , 5o] of the two vector fields 5 A and 5o on T is again a vector field 
(the Lie bracket of the two), so, since g is an exact stationary point of /, [5 A , 5q]I = 
(f)[[5 A ,5 ]\, and 

5 A = 5 a <I>[5o]-5o<I>[5a]-<I>[[5a,5 q ]\ (26) 
= &A<f>[5 A ,5 ]. (27) 

Here dl is the exterior derivative operator on JF, so the right hand side of f[2"Tj) is the 
curl of evaluated on the pair of vector fields S A and *H Its value at S depends 
only on the values of these vector fields at g and not on how they vary away from 
g. This is because the commutator term in (I2T|) cancels the contributions to the first 
two terms coming from the action of 5 A on 5q and vice versaWl 

Now consider the case in which 5a is the retarded perturbation 5\, which vanishes 
outside J + [sa] in the globally hyperbolic spacetime M. The achronality of Af requires 



30 In general, if (5 is a 1-form and v i and V2 are vector fields then 

dAP(v 1 ,v 2 )=d Vl 0(v2)-d V2 /3(v 1 )-P([v 1 ,v 2 ]) (28) 

where V2] is the Lie bracket of the vectors defined by [^1,^2] = d Vl V2 — d V2 vi or equivalently by 
the requirement that it's action (on a scalar) is equal to the commutator of the actions of V\ and v 2 '- 

31 If the reader prefers, the action of 6a on 6q and the action of <5o on 6 a, and thus also the 
commutator [6 a, <5o] , can be set to zero by taking the variations 8a9p.u and Sog^ to be independent 
of the metric on which they are evaluated, that is, by taking Sa and 5q to be constant vector fields 
on the vector space T . 
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J + [sa] to be disjoint from M (otherwise a future directed timelike curve from M to a 
point in sa and a future directed causal curve from that point to M could be deformed 
to a timelike curve from M to M). Furthermore, the facts that sa is compact and M 
globally hyperbolic imply that J + [syi] is closed ( |HE73j prop. 6.6.1.) so it is disjoint 
from an open neighbourhood of M . The boundary integral (p hi f[2"T|) may therefore 
be replaced by its restriction to £+, denoted Os + . Then 

M = <ffAe E+ [a+ £ ]. (29) 

(Since is an integral of a local functional of the fields - see (fl9l) . its functional 
derivatives are local functionals of the fields, so dl A 0[<5j[, 5q] receives no contribution 
from the parts of dQ outside the (closed) support of <Jj[.) In principle the restriction of 
the boundary integral to just the portion S + of the boundary is not unambiguous. 
The integrand of is defined only up to total derivatives, that integrate to zero 
because ddQ = 0, but which can contribute boundary terms to 6e + (see |LW90] ) . 
However, these terms do not affect (]29j) because 5\ vanishes in a neighbourhood of 
<9£+ = eW C Af. 

Now, we also have - for the same reasons - that the advanced perturbation 5^ 
vanishes on an open neighbourhood of the future boundary E + . Thus, on E + , the 
variation 5\ may be replaced by 

Aa = 5\ - 5~ A . (30) 

Equation fl2Ul) then becomes 

M = dlAe E+ [A jl ,5o] = nx + [A A ,6 ], (31) 

where is the presymplectic 2-form on JF associated with the hypersurface E + . 
For any pair of variations, 8\ and 82 , 

^[SiM = dl A9 s+ [5i,5 2 ] = 5i6 s+ [(5 2 ] -8 2 Qx + [5 1 } - e s+ [[5 h 5 2 ]\. (32) 

(pTL]) holds for any 5 in the space of solutions to the unperturbed linearized field 
equation on Q, ( l23i) . All the more so it holds for a variation 5q G L 5 which satisfies 
the linearized field equation on all of M. 

Notice that the variation A A g also satisfies the linearized field equation (j23|) : 
Subtraction of equation (I2"2l) with 5a = 8^ from the same equation with 5a = 5\ 
yields AaSI = Aa4>[8] for all 5, which is just (1231) . Indeed, since this is true also 
when Q is replaced by any compact subset of M with smooth boundary, G L g p 

32 If g is linearization stable this means that , with a suitable "diffcomorphism gauge fixing" , is 
tangent to S. Inj^Uand the following developments S A can be any retarded/advanced perturbation 
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(This last result also follows directly from (E]).) Indeed, since A A B = {A, B} for any 
observable B, Aa must be the Hamiltonian vector field of —A. 

Thus, if Q is defined to be the restriction of the presymplectic 2- form to L g , 
then dSTj) implies tha1@ 

6A = n[A A ,6) V8eL g . (33) 

f)33p is the equation we were looking for. It is the key equation relating the Peierls 
bracket (defined by A^) to the presymplectic 2-form. 

f determines A^ up to degeneracy vectors of Cl, that is, up to vectors 7 such 

O[ 7 ,5] = V6eL g . (34) 

This is sufficient to determine uniquely the Peierls bracket, {B, A} = AbA, of A 
with any observable B with domain of sensitivity contained in Q: (I33p holds for any 
5 £ L g . It therefore holds in particular when 5 = As- Thus 

{B,A} = A B A = n[A Al A B }. (35) 

Clearly the right side of this equation is unaffected by the addition of degeneracy 
vectors of Cl to A^ and A B . 

fl need not be defined using S + . One obtains the same 2-form on L g if one defines 
it using Af in place of £+. 

Proof, (ft = @s+ — @Af- Thus, if 61 and S 2 belong to L g , then@ 

nv + [5i,5 2 ] -OvPi.fc] = 5 1( j ) [5 2 ]-5 2( f>[5 1 ]- ( l>[[5 1 ,5 2 ]] (36) 

= 5 X 5 2 I - S 2 SJ - [6 U 6 2 }I = 1. (37) 

Actually, since 0£ + [5] does not depend on the shape of the remainder, dQ — E + , of the 
boundary of Q, any piecewise differentiable hypersurface E which together with E + 
encloses a compact spacetime region on which Einstein's field equations hold can be 
used instead of N or E + to evaluate VL[5\, 5 2 ], and gives the same result. In particular 

to the metric generated by A (i.e. by the source IQnrGa in the linearized field equation). It need 
not satisfy any further restriction, such as the transverse gauge condition, and so neither does A^. 
Diffeomorphism generators with support contained in the causal domain of influence of sa may be 
freely added to A^- If this freedom is exploited to ensure that A^, viewed as a perturbation of the 
metric, does not move the boundaries of the domain of dependence of J\f in M then it is tangent to 
S provided g is linearization stable. 

33 Restricting So in pip to lie in L g implies no loss of generality since any solution of the linearized 
field equations on Q can be extended to a solution on all of M. 

34 In fact the ambiguity in A^ is further restricted by the definitions of 8\ and , which require 
that A a vanishes outside of the causal domains of influence J + [s y i] and J^s^] of sa.. 

35 Recall that the commutator of the actions, as derivative operators, of two vector fields is equal 
to the action of their Lie bracket, which is itself a vector field. 
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any smooth Cauchy surface of T> will do. [This does not exhaust the possibilities. 
In principle the hypersurface need not even be achronal, or contained in T>. On the 
other hand our proof only guarantees the validity of (1331) if sa C T>.] 
From (fT9"j) one immediately obtains Q for the Einstein-Hilbert action: 

n^h] = 5 1 Qj,[5 2 ]-5 2 Qj:[5 1 ]-Qj:[[5 1 ,5 2 ]] (38) 

/ S^Mg^V^dZ, - (1 <- 2) (39) 

07TG JS 

This agrees (up to an over all factor) with the expression proposed by Crnkovic and 
Witten |CW87j . Moreover, when S is spacelike (]3"9"j) reduces to the standard expres- 
sion in terms of ADM variables |ADM 62J [LW90] : In Gaussian normal coordinates 
r a based on X r° = on X and d/dr° is the unit future directed tangent to the 
congruence of geodesies normal to X. Consequently g 00 = goo = — 1, g° l = go% = 0, 
and gij = hij where indices k, ... range over {1, 2, 3} and h is the induced 3-metric 
on S. Moreover r° = P 00 = and T°. = K ij} Y) Q = ti k K kj where K {j = l/2d h i:j is 
the extrinsic curvature of S. The presymplectic 2-form is thus 

nv[6 u 6 2 ] = j^Jj 2 K5 l Vh + 5 2 K t Mh l3 ^h}d 3 r - (1 <-> 2), (40) 

where K = Kijh 1 ^ and h = det[hij]. In terms of and the ADM momentum 
variabld^l 

n ij = \fh[K ij - Kh ij ]/WirG (41) 
the presymplectic 2-form becomes 

tov[5i t S 2 ] = J Stf'Szhijtfr -(1^2). (42) 
2.2.2 The Peierls bracket is a Poisson bracket 

(135p allows us to verify that the Peierls bracket {•, •} really is a Poisson bracket on 
the observables. In order that {•, •} be a Poisson bracket it must |Mar92j 

i be skew symmetric, 

ii be linear in each argument (under addition and multiplication by a constant), 

iii satisfy the Leibniz rule {A, BC} = {A, B}C + B{A, C}, 
and 

36 The normalization of ir corresponds to the normalization of the action © used here. 
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iv satisfy the Jacobi relation, 

{A, {B, C}} + {B, {C, A}} + {C, {A, B}} = (43) 



immediately shows that i, ii and iii are satisfied. The last requirement, the 
Jacobi relation, follows from the fact that Q, being the curl of 0, is itself a closed 
2-form. For any 2-form a and triplet of vector fields u,v,w 

Cyc d u a(v, w) = d A a(u, v, w) + Cy^c a([u, v],w) (44) 

(where Cyc denotes the sum over cyclic permutations of the sequence u, v, w). From 
(1551). (gtj)' and the fact that cl A ft = it follows thaS 



X = Cjc {A, {B, C}} = -Cyc A A n[A B , A c ] = -Cyc ti{{A A , A B ], A c ] (47) 
and therefore 

X = Cyc [A A ,A B ]C = Cyc {A, {B,C}}-{B, {A, C}} (48) 

c (X -{C,{A,B}}) = 3X -X = 2X (49) 



A,B,C 

Thus X = 0, which is precisely the Jacobi relation 
2.2.3 Definition of the bracket on initial data 

We turn now to the problem of defining a bracket on the initial data on A/", to be called 
the auxiliary bracket and denoted by {•,•}*, which reproduces the Peierls bracket on 
observables when these are expressed in terms of the initial data. 

A skew symmetric bracket, {<fi,<fj}*, defined on all pairs pi, ipj of initial data 
determine a bracket 

{^GK = £-— {v^,}*-— . (50) 



37 Actually it is not necessary to invoke a differentiable structure on T to obtain the Jacobi 
relation. The crucial relation 

C^c 6 3 ] = Cyc n[[S 1 ,S 2 ),6 s ], (45) 

that we have obtained from (T4"4"|) and cK A = 0, follows algebraically from (f32l) - 

n[^,^] = 5ie E+ [5 2 ]-<5 2 e s+ [<5 1 ]-0 s+ [[<5 1 ,5 2 ]] V6 u 6 2 eL g . (46) 
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on all pairs F, G of differentiable functions of the initial data. Here the indices /, J 
comprehend all continuous and discrete labels of the initial data components (includ- 
ing position in A/"), YjI,j represents the integral and/or sum over the set of such labels 
and 3/ dip i represents ordinary partial differentiation or functional differentiation as 
appropriate. {•,•}* as defined by f[5*Ul) is skew symmetric, linear in its arguments, and 
satisfies the Leibniz rule. This makes it a pre-Poisson bracket - a Poisson bracket 
save that it may not satisfy the Jacobi relation. 

This bracket is defined on observables only if they are differentiable functionals of 
the initial data. That they are can be demonstrated if it is assumed, as we have done 
elsewhere, that Rendall's [Ren90j results on the Cauchy development of null initial 
data apply to all of Af, and not just to a neighbourhood of So in Af. Rendall shows 
that in the Cauchy developments corresponding to a differentiable one parameter 
family of Sachs initial data the metric at a given harmonic coordinate point is also 
differentiable in the family parameter. The same is then easily seen to be the case 
for the closely related initial data we shall use (defined in subsection 13. 3j) . 

The derivative, 5g^, of the metric field in the family parameter is a solution to the 
linearized field equation, so the corresponding derivative of an observable A satisfies 
5 A = Q_^[Aa, S}. The linear operator Q_^[Aa, •] is a distribution on the derivative 5ipi 
of the initial data. Indeed it can be seen at once from the explicit form (11951) of Qjy 
in terms of our initial data, and the smoothness of A^, that Qj^[Aa, ■] is a regular 
distribution, equivalent to integration of the variations of the initial data against 
smooth functions. Since Q/j-[Aa, •] depends only on the definition of A and on the 
unperturbed metric, and not on which one parameter family of initial data is being 
considered, it is in fact the functional derivative of A by the initial data. dA/difj 
thus exists and is a smooth function of the the continuous position label in the index 
I. 

Note that A^ vanishes in a spacetime neighbourhood of dAf (because both 5^ and 
5\ do so). We will take care in the definition of our initial data to make sure that 
variations of the metric that vanish near dAf also leave the initial data undisturbed 
near dAf, so Aj^pj = for data in some neighbourhood of dAf. Since D,j^[Aa,S] 
is an integral over the initial data and their variations under A^ and 5 on Af (see 
( 11951) ) the variations Sipj corresponding to points sufficiently close to dAf will not 
contribute. That is, the derivative dA/dtpi vanishes for / corresponding to points in 
a neighbourhood of dAf. 

A sufficient condition ensuring that {•,•}* reproduces the Peierls bracket on ob- 
servables is that the motion {A, generated by any observable A via the bracket 
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satisfie@ 

5A = ti[{A, ■}*,§] (51) 

where L° g is the space of variations satisfying the linearized field equations and the 
requirement that dg^ vanishes in some neighbourhood of dAf in M. Then 

{A, B} = A A B = Q[{B, ■}*, A A ] = -Q[A A , {B, •}*] = -{B, A}, = {A, B},. (52) 

for all observables A and B. The first equality is just the definition of the Peierls 
bracket; the second follows from (1511) . using the fact that A^ vanishes in a neighbour- 
hood of dj\f; the third from the antisymmetry of 0; the fourth from (13"3l) ; and the 
last from the assumed antisymmetry of {-,-}*■ 
Condition (15TT) is satisfied if 

6tp I = n Ar [{<p I ,-}„8\ (53) 

for all initial data tpi living on the interior of M . Indeed it is sufficient that ( 1531) holds 
when both sides are integrated against smooth functions that vanish in a neighbour- 
hood of dj\f, for (1511) is obtained from (1531) by contracting both sides of ( 1531) with 

dA/d<pj. 

(153"]) even implies that {A, -}* = A^ distributionally on initial data on the interior 
of N '. If 9 is initial data smeared with a smooth function supported in intM then 
{9, •}*, like {A, •}*, corresponds to an element of L g . One may therefore repeat the 
arguments of equation ( 1521) . from the second equality on, with 9 in place of B, and 
conclude that A A 9 = {A,^}*. (This result also follows from (I353p and (1531) .) 

Notice that both (15T|) and (153"]) . by ensuring that {•,•}* reproduces the Peierls 
bracket on observables, also ensure that the Jacobi relation holds for {■,■}* on ob- 
servables. Nevertheless neither ( 15T1) nor (1531) imply that the Jacobi relation holds for 
brackets between arbitrary functions of the datao Indeed the bracket solving (15T1) 
that we find does not satisfy the Jacobi relation. As long as one is only interested in 
reproducing the Peierls Poisson algebra of the observables the validity of the Jacobi 
relation for brackets of the initial data is evidently not necessary, but if one wishes to 
quantize the system by replacing these brackets by commutators this relation becomes 
indispensable, since it holds identically for commutators. 

38 {^4, •}* is a variation of the initial data. That it necessarily corresponds to a solution to the 
linearized Einstein equation of course follows from the extensions of Rcndall's results Ren90J that 
we have assumed, but it can also be seen directly from the explicit formula for the solution (|367|) 
obtained in appendix [Cl 

39 One might think that the Jacobi relation follows from (j55)) . which essentially requires that {■,■}* 
be inverse to Cl. As we saw in sub-subsection [?.2.2[ the analogous condition ([55)1 docs imply the 
Jacobi relation for the Peierls bracket between observables. However (|55)) applies only to variations 
d vanishing in some neighbourhood of dAf, so {■,■}* is only almost inverse to Q, which turns out not 
to be sufficient. 
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Our task will therefore be to find a skew symmetric bracket on initial data, 
{<fii,<fj}*, satisfying (1531) . Two further conditions will be imposed on this bracket. 
The first of these relates to constraints. We shall calculate the bracket on a set of 
initial data which are almost, but not entirely, free. They are subject to one con- 
straint at the intersection 2-surface Sq of the two null branches of A/". The brackets on 
these almost free data then allow us to evaluate the brackets of a second set of data, 
which are completely free. As is always the case for constrained data the bracket on 
the almost free data defines an equivalent bracket on the completely free data iff it 
respects the constraints. That is to say, we must require that the functions of the 
data that are constrained to vanish also have vanishing brackets with all data. 

The second, and last, condition is a requirement of causality. We require that the 
auxiliary bracket vanishes between data at points of M that are not connected by 
a causal curve. Since disjoint generators are not connected by causal curves (prop. 
IB. 71) it follows that the only non-zero brackets are those between data on the same 
generator, which facilitates considerably our calculations. 

Whether the causality requirement is independent of (|53|) or implicit in it is unclear 
at present. The causality requirement seems quite reasonable in view of the causality 
property of the Peierls bracket - namely that the Peierls bracket between observables 
with domains of sensitivity unconnected by any causal curve vanishes. As will be 
clear from the definitions of our data, the domain of sensitivity of a datum at a point 
p E A/", as a functional of the spacetime metric@ is contained in the generator through 
p, or the pair of generators through p in the case that p G Sq. Thus the Peierls bracket 
between data on distinct generators, if it can be defined, should vanish. Indeed this 
can be done. ( fTUl) allows us to define the Peierls bracket between initial data that 
do not live on the same generator, provided a suitable fixing of the diffeomorphism 
freedom in the data and the Green's functions is adopted, and this Peierls bracket 
must vanish because the advanced and retarded Green's functions that enter ( TlOl) 
vanish between points that cannot be connected by causal curves. That the causality 
requirement is at least compatible with the other requirements follows a posteriori 
from the fact that a bracket satisfying all the requirements is found. 



2.2.4 Degeneracy vectors of the presymplectic 2-form and the consistency 
of the definition of the auxiliary pre-Poisson bracket 

The consistency of (I53p will be fully established post hoc by the fact that a solution to 
this condition is found. Nevertheless it is interesting to attempt to verify consistency 
already at this stage, before entering into the details of a particular representation 
of the solutions (i.e. of the initial data). We shall succeed in doing so in the generic 

40 Quite generally the domain of sensitivity of a functional F of the metric will be defined to be 
the support in spacetime of its functional gradient SF/Sg^ . 
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case in which a neighbourhood of each point can be coordinatized by scalars defined 
by the geometry. 

The consistency of (153j) requires that all data tpj be invariant under the action of 
degeneracy vectors 7 of such that the corresponding variation of the metric, 7$, 
vanishes in a spacetime neighbourhood of dM. At least on generic solutions these 
turn out to be diffeomorphism generators, i.e 7 = £^ for some vector field £ on 
In fact the following stronger result holds 

Proposition 2.1 In the generic case, in which there exist at each point of V four 
smooth functions of the metric g and a finite number of its derivatives, having smooth, 
linearly independent spacetime gradients, any variation 7 G L g that satisfies the gen- 
eralized degeneracy condition 

fi( 7 , 5) = V5 G L° g (54) 
is a diffeomorphism generator on V. 

Proof. Any divergenceless, C°°, symmetric, bivalent tensor field a of compact support 
contained in T> defines, via ([H]) and ([9]) a vector G L g which, by ( |33|) . must satisfy 

J bg^aTyf^gtfx = Cl[A A ,S\ V5 G L g . (55) 

If S is a generalized degeneracy vector 7, satisfying f[5"4"|) then the right side vanishes 
and we have 

/ 19^, a^V^g d 4 x = (56) 

for all a satisfying our conditions. 

Let p be a point of V, and z a , a = 0, 3 four smooth functionals of the metric 
with linearly independent gradients dz a . By the inverse function theorem these form 
a smooth chart on a neighbourhood U of p in M D V. The components of the 
metric in this chart are diffeomorphism invariant, so differentiating any functional of 
these components by the components of the metric in a manifold fixed chart yields 
a symmetric, divergenceless tensor. Generalizing this observation we shall consider 
linear functionals of the form 

F[5g] = J ^r P ^g~ z d A z (57) 

where g a p are the components of the metric in the z chart and f a/3 is an arbitrary 
symmetric matrix of functions with support contained in UT\T>. The variation of g a p 

41 The tangent to a one parameter family of diffeomorphisms is the Lie derivative £^ along the 
vector field £ of tangents to the orbits of the manifold points. 
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is related to the variations of the components, g^, of the metric in a manifold fixed 
chart x by 

5g aP = dza [Sg^ - £ m g^} (58) 



dx^ . . f 5z a (q) 



where 

= ^(9) y ^W^Jy^W (59) 

at each q E U. (q is fixed in the spacetime manifold, ^"(g) thus varies when the 
metric varies.) Thus 

F[Sg] = J[Sg,v-2V^[5U}rV^d A x (60) 

= J{5g,,r + 2m p 9 P ^jnv^d A x (ei) 

= /WV^^ (62) 



with 



and 



1" = (63) 



/""CO = /^(r) + -= / -^—^-g^f-r^d'x (64) 
y^gxJ 5g^{r)dz a 

This tensor is smooth because the functional derivative is a sum with smooth coeffi- 
cients of derivatives, up to finite order, of delta distributions, and this distribution is 
convolved with a smooth function. Clearly the support of / is contained in that of /, 
which in turn is contained in U D T>. Finally, / is divergenceless. This can be verified 
by letting 5 = £( with £ smooth and of compact support contained in U: 

F[£(] = - J C, u V,r V=g~d 4 x (65) 

vanishes because £[£C] = C> since z a are functionals of the metric only, and thus the 
variation of the z components of the metric vanishes. ( 0651) can also be verified by 
direct calculation.) (1651) for all smooth, compactly supported £ implies V^/^ = 0. 

Thus F[y] = for all f a/3 , (which are smooth and have compact support contained 
in U n T>) which implies that 

o = jgap = ig^v + £tYi\g»v (66) 

That is 7 = a diffeomorphism generator, on U D T>. Clearly this implies that 7 is 
a diffeomorphism generator throughout D. 
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In fact 7 is a diffeomorphism generator also on intM. By the Whitney extension 
theorem (see |AR67j ) 7*7 can be extended smoothly to a neighbourhood of intM in 
M. £[7] is defined on the intersection of U with this neighbourhood. Therefore, since 
both 7 and define smooth variations of the metric, jg = £fy\g on U H intM 
which is a subset of the boundary of D in the topology of U. I 

The genericity assumption does not actually seem to be necessary. In particular, 
the result holds in a flat spacetime, which does not satisfy this assumption. 

In any case, we may conclude that in the generic case that a generalized degeneracy 
vector is a diffeomorphism generator on V. Conversely, a direct calculation (which 
will be useful for another purpose further on)@ shows that 

Proposition 2.2 Any diffeomorphism generator is a generalized degeneracy vector 

Proof. Suppose £ is a C°° vector field on P, 5 is a variation satisfying the linearized 
vacuum field equation, and E is a smooth Cauchy surface of T>, contained in T>. Note 
that a smooth Cauchy surface exists by [BS03] . and that by prop. IB. 151 the closure 
E of E in M is an embedded compact 3-manifold with boundary dM. For the sake 
of simplicity (it makes no difference to the result) suppose £ is independent of the 
metric, and thus 5£ = 0. Then 

Q[£^,5] = £^[5}-5Q s [£^]. (67) 

Now 

©sW = [ srtg^V^gcfc, (68) 

a (69) 



where a is the 3-form -1/(3! 8nG) ST^g^ 'e^y/^g dx x A dx p A dx T . Thus 

£^[5}= f £^a= f £j[dAa]+ f dA[£j«]. (70) 

But d A a = 1 / {l&nG)8R^/—g dx° A dx l A dx 2 A dx 3 , the divergence term in the 
variation of the Einstein-Hilbert Lagrangean density, which vanishes because g G S 
and 5 £ L g and thus R = 5R = 0. Therefore 

£&z[8\ = (71) 

eST^g^V^gdE^ (72) 

/_ + Jv CT ^^v^) (73) 



/as 
1 



12 



8nG JdT. 

A more general argument is given in subsection 
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Here dY*^ v = 1/2 e M ^ pr dx p A dx T is the coordinate area 2-form, and in the last line the 
identity 

w% = ~<T{v^ + v^r, - v r ^}. (74) 

has been used. 

The second term in ( 16T1) is the 5 variation of 

But (1T1"]) and Einstein's field equation, which is satisfied by g, imply that 

£Ag^V=g = (v*v ff r - iv,v*e - ~v,v*o^ (76) 

= V CT V^r ] v^ (77) 
= d a [VH a] ^- (78) 

Taking the integral one obtains 

Summing ( 1731) and minus the S variation of ( 1791) one obtains 

to[£t,S\ = [_e{W»5^g- + VJg™V^} + 5{V^^}dZ^ (80) 

1d7tG Jas 

= TTrV / - ^T^ 1 V=0 + 5[y/=g^] V^dE^, (81) 

lD7rG JdT, 

Since this is a boundary integral it certainly vanishes if S vanishes in a neighbourhood 
of <9£ = dj\f, showing that any diffeomorphism generator is a generalized degeneracy 
vector of Cl. I 

In order that £^ be a true degeneracy vector the boundary integral (]8ip must 
vanish for all 5 G L 5 . A sufficient condition is that £ and it's gradient vanish on <9£. 
It is not clear to the author whether this is a necessary condition. 

According to prop. 12.11 when the metric g on V is generic a degeneracy vector of 
Cl which preserves the metric in a neighbourhood of ON is a diffeomorphism gener- 
ator which vanishes in this neighbourhood. Our initial data will be invariant under 
all diffeomorphisms which preserve cW, so they are certainly invariant under such 
degeneracy vectors, implying that (15*3"]) is consistent in this case. 

On the other hand on any solution, even one with a non-generic metric, (J53l) only 
defines {<pi, •}* up to a generalized degeneracy vector. Thus, according to prop. 12.21 
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one is free to add a diffeomorphism generator. This ambiguity is to be expected if 
{■,■}* is determined only by the requirement that it reproduces the Peierls bracket 
on observables, because observables are diffeomorphism invariant. 

The consequent ambiguity of the brackets between our initial data is not very great 
however, because the data are invariant under a large subgroup of diffeomorphisms. 
Nevertheless, the freedom to add diffeomorphism generators to the solutions of ( |53l) 
will be used to simplify calculations further on. 

Finally, the fact that any diffeomorphism generator that vanishes in a neighbour- 
hood of dM is a true degeneracy vector of Q implies that we may, without changing 
the consequences of ([SSI), add a suitable diffeomorphism generator to each variation 
in L° g such that the total variation preserves the null character of the branches of M. 
This partial "gauge fixing" of L° g , which will be adopted from here on, simplifies our 
subsequent considerations, principally because the fixed hypersurface M then carries 
null initial data for the varied solutions and the variation of this data is simply a 
function of the variation of the metric (and of its derivatives) on M. 

It is tempting to try to go further and restrict the variations so that they leave 
the generators of M fixed (i.e. the null generators for a varied solution are the same 
as in the fiducial solution). This can be done by adding suitable diffeomorphism 
generators £^ to the variations in L g , but in general £ cannot be chosen to vanish in a 
neighbourhood of dM, so the resulting variation will not lie in L° g : Imagine following 
the generators back in time from the boundary dM to So, that is, from the disks Sl 
and Sr where the generators of Ml and Mr, respectively, cross the boundary. (See Fig. 
CD). A perturbation which vanishes in a neighbourhood of dM will leave the generators 
fixed in that neighbourhood. However, it will in general affect which generator from 
Sl meets a given generator from Sr at Sq. So no diffeomorphism acting only in 
the interior of M can restore the generators to their unperturbed positions. Aside 
from the exceptional cases in which the the geometry near dM admits non-trivial 
isometries, these are the only diffeomorphisms that leave the metric invariant near 
dM. Thus, while a perturbation of the metric that preserves the geometry near dM 
can always be reduced by the addition of a suitable diffeomorphism to a perturbation 
that vanishes (i.e. preserves the metric tensor) in a neighbourhood of dM, it is not in 
general possible to add a diffeomorphism so that the total perturbation leaves both 
the metric invariant near dM and the generators invariant on all of M. 

We shall continue to use variations 5 which vanish near dM in (153jl . and the 
generators will not in general be fixed under under these variations, even though the 
hypersurface M they collectively sweep out is. 
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2.2.5 The auxiliary bracket for spacelike initial data 

Our conditions also serve to define a bracket on spacelike initial data, which turns 
out to be the same as the Poisson bracket used in the ADM formulation of canonical 
GR |ADM62j . This implies that the Peierls bracket between a pair of observables 
of the domain of dependence of a Cauchy surface E is equal to the bracket between 
these in the ADM theory. Our theory is thus equivalent to the ADM theory. 

The auxiliary bracket is a solution to ( 1531) . but it is convenient to first solve the 
simpler, but analogous condition 

5(pi = £lx[{<pi, -}a> $] V<5 vanishing near <9E. (82) 

This condition is stronger than (1531) because it must hold for all 5 vanishing in a 
neighbourhood of the boundary, not just those satisfying the linearized field equations. 
On the other hand {pi, -}a will not be required to satisfy the linearized field equations 
either (for this reason f2 must be used instead of 0)@ When E is a spacelike Cauchy 
surface of V one may use the ADM variables, consisting of the induced 3-metric h 
and the momentum variable n defined in ( 14TT) . as the initial data. In terms of these 
variables the presymplectic 2- form is given by the simple expression ( 1421) . It follows 
immediately that the unique solution to (1H2"]) is the standard ADM kinematic Poisson 
bracket: 

{h ll (x),7r kl (y)} A = 5f l 6 l j) 6 3 (x,y) (83) 

{M*)> h k i(y)U = 0r <J », *t*'(2/)}a = (84) 

for all x and y in the interior of E. 

Now let's return to the auxiliary bracket {•, •}*, which satisfies (|5"3"|) . (1531 requires 
that 5ipj = fin [{</?/, •}*, 5} holds for all 5 that vanish near <9E and satisfy the linearized 
field equations (that is V<5 G L°). Moreover {<fi,-}* is also required to satisfy the 
linearized field equations. These conditions are equivalent to requiring that the actions 
of 5 and {ifi, ■}* annihilate the constraints on the data on E (for if a perturbation of 
the initial data satisfies the linearized constraints then there always exists a solution 
to the linearized field equations matching the perturbation) @ If 6 annihilates the 

43 '}a is (a generalization of) the kinematic Poisson bracket defined in |LW90j . 

44 Proof: The sourceless transverse gauge linearized field equation 

V ff V CT7Al „ - 2i?V T 7 CTT = (85) 

for the trace reversed metric perturbation = bg^ v — \g^u{g' Tf '5g CTp ) determines a unique solution 
on the entire domain of dependence of £ given any smooth symmetric tensor fields 7^ and V n7At „ 
on E, where n is the unit normal to S. (see (fTH|) .) If these data satisfy the transverse gauge condition 

V"V ee X u = (86) 



37 



constraints then it follows from (1821) that for any constraint C 

n E [{C,-} A ,5] = 0. (88) 

Thus, if we could add some combination of actions of constraints {Ci, -} A to {<pi, -}a 
so that the sum annihilates the constraints then we would have a solution {ipj, ■}* to 
(loll]) . This can be done, at least formally: A gauge fixing is chosen that completely 
fixes the gauge transformations generated by the first class constraints, and 5 is 
required to respect this gauge - that is, to annihilate the gauge fixing constraints. 
Then the Dirac bracket, 

{-,-}* = {-,-}a - £{-, aUX^iCj, -} A , (89) 

fulfills our requirements. Here the set {Ci} includes all constraints, also the gauge 
fixing constraints, and X is the inverse of the matrix {Ci, Cj}&. The Dirac bracket 
is clearly a pre-Poisson bracket. In fact it also satisfies the Jacobi relations^ Note 
however that several issues arise when this construction is attempted for general 
relativity, which posses an infinite number of constraints. (A finite set of constraints 
at each point of £). An inverse X to {Ci, C.,}a may not exist even if the latter matrix 
is non-degenerate - i.e. has no degeneracy vectors; the sum in (1891) may not converge; 
the sums in the proof of the Jacobi relations may not converge. 

All the constraints on the ADM variables h and it are first class and generate, via 
the kinematic bracket {•, -}a, diffeomorphisms of the Cauchy surface £ in the solution 
spacetimc [ADM62J. The effect of such a diffeomorphism of £ on the the data on 
£ is of course equivalent to that of leaving £ fixed and applying a suitable diffeo- 
morphisms to the spacetime metric field. Since the observables are invariant under 
such diffeomorphisms the kinematic bracket between an observable and a constraint 
must vanish. Since X has non-zero matrix elements only between the original first 
class constraints and the corresponding gauge fixing constraints it follows that the 
kinematic bracket between two observables, is equal to their Dirac bracket, which in 

on £ and the linearized constraints 

= n^SG^ = n^V^Xu) ~ l/2.g M „V ffX CT - 1/2(V ct V ct7 ^ - 2i?V T 7ffT )] (87) 

then in the corresponding solutions to (|85[) also the normal derivative of the transverse gauge con- 
dition holds: V„Xi/ = 0. But taking the divergence of (|55|) we find that this equation implies 
V^Vo-Xj, = 0. The fact that \ an d its normal derivative vanish on £ thus implies that it vanishes 
on all of the domain of dependence. But then (|85p is equivalent to the linearized Einstein vacuum 
field equation. 5g^ v = 7 M „ = 7 M „ — l/2g^g' JT ^/ crT thus solves this equation. 

45 To verify the Jacobi relations for the Dirac bracket note that the kinematic bracket {•, - }a docs 
satisfy the Jacobi relations, and that {A, = {^4*, £?*}a where A* = A — ^^{A, Ci] /\X l ' j 'Cj and 
B* is defined similarly. 
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turn, by (152j) . is equal to their Peierls bracket. In other words the Peierls bracket of 
the observables is equal to the standard ADM bracket between them. 

2.3 Gauge transformations 

Degeneracy vectors of the presymplectic 2-form are normally associated with gauge 
transformations. Though it is not strictly necessary for our purpose of obtaining 
the auxiliary pre-Poisson bracket between initial data it is interesting, and ultimately 
important for the theory, to see how this works out in the situation we are considering, 
in which the initial data is specified on a hypersurface with a finite boundary. Here 
a definition of gauge transformations is proposed and some of its consequences are 
worked out. 

Gauge transformations are transformations of the solutions that affect only those 
degrees of freedom of the field which are in principle unpredictable. A more precise 
definition (inspired by that of [LW90] ) would be the following: 

A map 9 from the space of solutions to itself is a gauge transformation iff for any 
pair Ci, C*2 of disjoint closed sets in spacetime there exists another map 9' from the 
space of solutions to itself that reduces to 9 on the field on C\ and to the identity on 
the field on C 2 . 

Put another way, the gauge transformed solution on any closed portion C\ of 
spacetime can always be joined to the untransformed solution on any disjoint closed 
region C2, forming another solution. Thus knowledge of the field on all of C2 does 
not suffice to determine whether or not the field is transformed by 9 on Ci. 

Here a slightly different definition will be proposed. As discussed earlier, the 
developments of smooth initial data are expected to be smooth at ON . However 
disjoint subsets of V which are closed in V may have closures in M which are not 
disjoint, containing common points on dM which is outside V. The requirement that 
gauge related solutions can be spliced on all pairs of disjoint sets that are closed 
in 22 combined with the requirement that the resulting spliced solutions be smooth 
at dM puts restrictions on the gauge transformations near dN '. Whether gauge 
transformations should be restricted in this way or not depends on the usefulness 
of the notion of gauge obtained. Here I opt to remove these restrictions on the 
transformations by placing further restrictions on the pairs of subsets C\,C2- The 
resulting set of gauge transformations is certainly simpler as a consequence. 

Definition 2.1 A map 9 from the space of solutions S to itself is a gauge transfor- 
mation iff for any pair C\, C2 of subsets of V having disjoint closures in M there 
exists another map 9' : S — > S that reduces to 9 on the field on C\ and to the identity 
on the field on C2. 
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Note that the manifold V here is the union of intM and the interior of its domain 
of dependence in some, arbitrarily chosen solution which we take as fiducial. (For 
instance V formed from Ao and the interior of its domain of dependence in Minkowski 
space could be used - see subsection 12.11 ) As in subsection 12.11 ^ consists of the metric 
fields of the maximal Cauchy developments, mapped via the diffeomorphism of prop. 
[R23l to V. 

According to the definition def. 12.11 (and also the previous one) the gauge trans- 
formations of general relativity are diffeomorphisms, as one would expect: 

Proposition 2.3 A gauge transformations of general relativity acts on any given 
solution as a diffeomorphism of V to V @ 

Proof: Consider a particular solution g to Einsteins field equations on V and a 
gauge transformation 9. Let C\ be a compact set in the interior T> of V and C2 the 
causal past in D of a spacelike Cauchy surface £2 of T> to the past of C\. C2 is closed 
in T> since it is the complement of the open set I + [T,2, T)\ in T>. Its closure in M thus 
intersects D only on C2 itself, and thus does not intersect C*i, which, being compact, 
is already closed in M. 

The well known uniqueness theorems for the Cauchy problem of general relativity 
(see |HE73j ) then show that the solution on C2 determines that on G\ up to a diffeo- 
morphism, which necessarily maps V — C2 to itself (since S contains only solutions 
in which V — C2 is the spacetime of the maximal future Cauchy development of E 2 ). 
It follows that on C\, on the given solution, 9 is equivalent to a diffeomorphism that 
maps C\ to a compact subset of V. 

In fact 9 is equivalent on all of T> to a diffeomorphism from T> to itself: If C[ is a 
compact subset of T> which contains C\ then the diffeomorphism on C\ corresponding 
to 9 may be chosen to be the restriction to C\ of the diffeomorphism corresponding 
to 9 on C[. If the metric on C\ has no isometries this follows from the fact that 
the mapping 9 of the metric determines a unique diffeomorphism. When there are 
isometries we are still free to impose this condition. Now, because T> is paracompact it 
may be covered by a sequence of such nested compact sets. If there are no isometries 
of all of T> then some compact set in the sequence, and all following it, will also be free 
of isometries so 9 defines a unique diffeomorphism on each set in the sequence, such 
that the diffeomorphisms on the larger sets are extensions of the diffeomorphisms on 
the smaller sets. 9 is therefore equivalent to a diffeomorphism : T> — > T>. If T> 
does have isometries a compact set in the sequence may be chosen that has only the 
isometries of D as a whole. The diffeomorphisms equivalent to 9 on the smaller sets 

46 A gauge transformation may act as one diffeomorphism on one solution and a different dif- 
feomorphism on another. Thus as a map from S to S it is not necessarily equivalent to any one 
diffeomorphism. 
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in the sequence as well as on the larger sets may then be chosen to be consistent 
with the diffeomorphism on this set, and again 9 is equivalent to a diffeomorphism 
(f) : T> — > V, although in this case it is not unique. 

This result extends to the bounding hypersurface intM. That is, on all of V, 
including intM, 9 is equivalent to a diffeomorphism from V to V. To show this we 
shall demonstrate, using the fact that both g and 9g are smooth on V, and thus can 
be smoothly extended to a boundaryless manifold Z G M containing V D intM, that 
a suitable extension of the diffeomorphism <p to intM can be constructed. 

Let p be a point of intM and let U be a convex normal neighbourhood of p in Z 
according to a smooth extension g of the metric g to Z. Let q be a point of U fl T> 
from which p can be reached via a past directed timelike geodesic £, and suppose that 
a Riemann ^-normal coordinate chart, y, covering U and based at q is set up. The 
map (f) takes q to q' G T> and the restriction of y to U fl T> to a new chart y' with 
origin at g' which is normal with respect to the metric 9g. An extension of the chart 
y 1 to a slightly larger domain defines, together with the chart y, the extension of <ft 
that is needed. 

First we need to understand a little more about the image under the diffeomor- 
phism of the domain U fl T>. Let £° be the segment of the geodesic I which lies in 
T>, and let 1°' = 4>{i°). £°' is timelike according to 9g and past directed, if the time 
orientation on V is always chosen so that intM lies to the past of T>, because with 
this convention preserves the time orientation of timelike curves. Since q' li es i n T^i 
and thus in the future domain of dependence, according to 9g, of intM in V, but i°' 
does not cross intM £° must be an extendible curve - it must have a past endpoint, 
p 1 , in V. p' cannot lie in D since 1° would then have a past endpoint in T>. p' thus 
lies in / M. (p' will turn out to be the image of p under the extension of <fi.) 

The same argument applies when p is replaced by any other point p in the open 
subset of intM fl U that lies in the chronological past I~ [q] of q. That is, the image 
under of the geodesic segment from qtop (not including p itself) has a past endpoint 
on intM. 

Let g' be a smooth extension of 9g from V to Z, then the past directed, time- 
like geodesies from q that define the normal coordinates y' near p' can be extended 
smoothly through their past end points on intM C Z. These extended geodesies can 
then be used to extend the domain of the normal coordinates y'. The geodesies do 
not form any caustics on intM near p': In terms of the metric components 

d_ d_ 

dy^'dy 7 ^ 



9U = 9'(^,^jp) (90) 



in normal coordinates a caustic occurs at a point if and only if det[g' a/3 \ approaches 
zero as the point is approached. But the y' components of the metric 9g at a point 
r G T> are the same as the y components of the metric g at _1 (r) G T>. There can 
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thus be no caustic at p' since the y chart does not break down at p, being good on all 
of U. The inverse function theorem then implies that the y' coordinates constructed 
from the extended timelike geodesies form a good chart of a neighbourhood of p' in 
Z. Since the diffeomorphism y'^ 1 o y coincides with (f) on the portion of its domain V 
lying in T>, and furthermore it maps no point outside D into T> and y'(p') = y(p), the 
definition = y' -1 o y on V extends as a diffeomorphism to a domain including a 
neighbourhood of p in Z. 

Of course <p(p) = p'. Moreover, the limiting values of the y components of g as p 
is approached from T> are the same as the limiting values of the y' components of g' 
as p' is approached from T>, which implies that the actual values of the y components 
of g at p and the y' components of g' at p' are equal, because the smoothness of 
the metrics g and g' imply that the charts y and y ; , and the corresponding metric 
components, are smooth. Thus 9g = <p*g in a neighbourhood of p in intN ' . 

<fi can be extended in this way to all of intN. If 4>i an d 4>2 are two extensions of <ft 
such that a point t G intN lies in the intersection of their domains of definition then, 
because 0i = 2 on V and both are continuous 4>\(t) = 4> 2 (t). I 

Now consider linearized gauge transformations. The tangent 7 G L g at a solution 
g to a one parameter family of gauge transforms of g clearly satisfies the following 
condition: For all pairs C\, C 2 CP having disjoint closures in M there exists 7' G L g , 
such that 7' = 7 on C\ and 7' = on C%- I shall call any 7 G L g that satisfies this 
condition a linearized gauge transformation. These can be characterized as follows: 

Proposition 2.4 Linearized gauge transformations are generalized degeneracy vec- 
tors of the presymplectic 2-form. 

Proof: The hypersurface independence of Q implies that any linearized gauge trans- 
formation 70 which vanishes in a neighbourhood of ON is a degeneracy vector of Cl: 
Let Ei, E 2 C V be smooth Cauchy surfaces of V, such thal0 

E 2 C/ + [E i; P], (91) 

and let C[ = J~[T,i,V] and C 2 = J + [E 2 ;£>]. These are disjoint closed subsets of 
T>. However their closures, C 1 , C 2 in M intersect on ON. A point p G C l D C 2 
lies in D[N; M] = D[E i; M] = D[E 2 ;M], and thus^by prop. 8.3.2. of |Wald84j any 
inextendible timelike curve through p must cross Ex, at p or later, and E 2 , at p or 
earlier. But E 2 C / + [S!; M] so the crossing of E 2 must be after or simultaneous with 
the crossing of Ei. The only consistent possibility is that p G Ei D E 2 , which by (|9T|) 
and prop. IB. 151 implies that p G dN. 

47 Time ordered sequences of Cauchy surfaces exist by Geroch's theorem, prop. 6.6.8. of [HE73] . 
and may be taken to be smooth by the result of BS03] . 
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Let C\ be the support of 70 in the causal past of Si. The closure C\ of this set in 
M is a subset of C x which excludes a neighbourhood of dj\f. C\ therefore does not 
intersect Ci. By the definition of linearized gauge transformations there must thus 
exist 7q G L g such that 7q = 7 on C\ and 7o = on C 2 . Therefore 

%o,-] = ^[7U = £M7o,-] = 0. (92) 

Now consider an arbitrary linearized gauge transformation 7 and a solution 5 
to the linearized field equations vanishing in a neighbourhood U C M of <9A/". By 
definition there exists another linearized gauge transformation 7' which is equal to 7 
onP-(/ (which contains the support of 5) and which vanishes on a closed subset of 
U which contains ON in its interior. Thus 

%,£]=%', ,5] (93) 

and, since 7' vanishes in a neighbourhood of cW, this vanishes by ( 1921) . This is valid 
for any 5 G L° so 7 is a generalized degeneracy vector of Cl. I 

This last result and prop. I2.1l imply that, at least when the metric on V is generic 
(in the sense of prop. 12.11) . a linearized gauge transformation must be a smooth 
diffeomorphism generator - that is, a Lie derivative £^ along a smooth vector field £. 
This can in fact be demonstrated directly, without the genericity assumption, by an 
argument analogous to that used above to characterize the gauge transformations of 
full general relativity, from the fact that data on a Cauchy surface of T> determines 
the solution to the linearized field equations up to the addition of terms of the form 
£i9iiv = 2V(^^) (see Appendix[C]). Conversely, it is easy to see that any Lie derivative 

along a smooth vector field £ is a linearized gauge transformation according to 
our definition. 

The linearized gauge transformations are thus precisely the diffeomorphism gener- 
ators, they are always generalized degeneracy vectors of Q, and, at least for a generic 
metric on V all degeneracy vectors are diffeomorphism generators and thus linearized 
gauge transformations. 

It would perhaps be natural to impose a further condition on the linearized gauge 
transformations that all tangents to full gauge transformations satisfy, namely that 
the integral curves of £ do not leave V. This would restrict the linearized gauge 
transformations to only a subset of the generalized degeneracy vectors of f2, although 
it does not require them to be true degeneracy vectors. 

The best - most natural and most useful - choice of definition of gauge transforma- 
tions in the context of a spacetime with boundaries is not clear to the author. In this 
section some possible definitions have been presented. Note however that nothing in 
the remainder of the paper depends on the choice of these definitions. Indeed the 
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general notion of gauge transformations is not used in the paper, precisely because 
no clearly preferred definition was found. Instead the classes of transformations con- 
sidered are identified explicitly in each instance (e.g. diffeomorphisms that preserve 
a neighbourhood of dM). 



3 Free initial data and adapted coordinates 
3.1 The area parameter 

The tangents tia to the generators that sweep out the branch A/"J3 of M are also 
normal to Ma (prop. IB.2[) . As a consequence Ma is a null hypersurface (its normals are 



null vectors) and the 3-metric induced on Ma by the spacetime metric is degenerate, 
as there is a vector in the tangent space - namely ua - which is orthogonal to all 
tangent vectors. Put another way, if one uses a tangent to the generators as the first 
vector in a basis of the tangent space of Ma then the 3-metric takes the form 




, 

Q n ij 



(94) 



where is a symmetric 2x2 matrix. 

The degeneracy of the 3-metric implies that the distance between neighbouring 
points p and q on Ma does not depend on the component of the (infinitesimal) dis- 
placement vector p~q parallel to the generators. Similarly the area, 



A=y/t 1 >t 1 t 2 >t 2 -(t 1 -t 2 ) 2 , (95) 

of an infinitesimal parallelogram in Ma spanned by vectors t\ and t 2 based at p G Ma 
does not depend on the components of t 1 and t 2 along the generators. It depends 
on p and on the set of generators cut by the parallelogram, but not on the "angle" 
at which it cuts these, that is, on the tangent plane of the parallelogram. This fact 
makes possible a very useful (non-affine) parametrization of the generators. It implies 
that the cross sectional area of an infinitesimal bundle of generators neighbouring a 
given generator 7, and terminating on a surface of area A in Sa, depends only on the 
point p G 7 at which the cross section is taken, and (linearly) on A. In other words, 
the factor A(p)/A by which the area expands following the generators from Sa to 
p depends only on p. The expansion factor can therefore be used to parametrize 7, 
provided it increases or decreases monotonically along 7. 



IS 



Here and from here on subscripts A, B, ... represent indices ranging over L and R. 
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We shall use the square root of the expansion factor, 



r = y/A/A, (96) 

to parametrize the generators. The cross sectional area of a bundle of neighbouring 
generators is then proportional to r 2 , which means that r is an area parameter. That 
is, r(p) is proportional to the distance from p to the caustic point, where neighbouring 
generators focus and A = 0, inferred from the area subtended at the caustic by a 
standard disk at p by assuming that this area falls as 1/ distance 2 as in flat spacetime. 

H 

As mentioned in subsection 12.11 we shall consider only initial data such that the 
generators of a given branch, Ml or Mr, of M are either everywhere diverging (A 
increasing) or everywhere converging toward the future (A decreasing). This ensures 
that r is a good parameter on each generator of the branch. 

Parametrizing the generators, and the initial data on them, by r makes it easy 
to incorporate the condition that Ml and Mr are truncated before their generators 
reach a caustic. One simply requires that each branch corresponds to an interval 
of positive, non-zero, values of r. @ The use of the area parameter also simplifies 
the mathematics of the null initial value problem. This was noted and exploited by 
Gambini and Restuccia [GR78] to obtain an expansion in powers of Newton's constant 
for the brackets between the main null initial data. 

49 In vacuum general relativity a bundle of neighbouring null geodesies is either incomplete or has 
at least one caustic, except in the special case that its cross sectional area is strictly constant along 
the bundle. (This follows from the Raychaudhuri equation discussed in the next paragraph.) Such 
constant cross section bundles do not exist in generic vacuum solutions (see prop. 4.4.5. of |HE73j 
and the discussion which follows it). If a bundle of null geodesies has several caustics, it is most 
natural to r to be a measure of distance to the nearest caustic in the direction along the bundle in 
which r is decreasing. Of course this is just an issue of interpretation. (|96|) defines r regardless of 
the interpretation adopted. 

50 This does not in itself guarantee that the generators do not cross in spacetime, just that 
neighbouring generators do not (more precisely, that the separation of neighbouring generators does 
not become zero to first order in their separation on So). But recall that we are considering Cauchy 
developments, and in a Cauchy development of initial data on AT, the generators of N never intersect 
(except of course at their starting points on So). [If N intersects itself in a solution to Einstein's 
field equations then (a portion of) the domain of dependence of Af in this solution is represented in 
the catalogue of Cauchy developments of N by a covering manifold (see appendix IBI prop. IB.6|) .] 

If the reader prefers not to represent spacetime by a covering manifold an alternative argument 
is that in order to calculate the brackets of the data along a given generator 7 one may use instead 
of Af an initial data surface Af' swept out by the generators emerging from a sufficiently small disk 
.Sq C Sq containing the base point of 7 in So such that these generators do not cross. 
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3.2 A geometrical description of the initial data 



We now turn to the initial data, beginning with a coordinate free description. All the 
data have an interpretation in terms of the geometry of the corresponding solution, 
but first let us consider M and the initial data on their own, prior to embedding in 
the solution spacetime. 

The set of initial data that we will consider consists of a 3-metric on M and further 
data specified on So only. The 3-metric is specified by two degenerate, positive semi- 
definite 3-metrics gi and qr defined on Ml and Mr respectively, gi and gn are 
required to be smooth and to match at Sq in the sense that they induce the same 
2-metric on S . The degeneracy subspace N£ of qa at every point p e Ma is one 
dimensional, and the integral curves of these subspaces, which will be called the 
"generators" of Ma, form a congruence that threads the cylinder Ma along its axis. 
That is, each integral curve originates on So and terminates on Sa- The congruence 
as a whole defines a (trivial) fibration of Ma with the generators being the fibres and 
S the base space. 

The data on So are most easily described in terms of two fields, x an d oj, that 
depend on a choice of parametrization of the generators. We shall therefore suppose 
that a good parametrization has been (arbitrarily) chosen on the generators. The 
parametrization defines fields of tangent vectors on Ml and Mr which will be denoted, 
like the tangent vectors to the parametrized generators in the spacetime context, by 
til and ur respectively. 

X is a smooth scalar field on Sq. uj is a smooth 1-form field on So, co-tangent to 
So- Under a change of parametrization such that ul and ur are rescaled by factors 
oll and (Xr 



where a -<= beneath a form indicates that the pullback of the form to Sq is to be 
taken. 

In the corresponding solution to Einstein's vacuum equations the 3-metric on M 
matches the 3-metric induced on the embedding of M by the spacetime geometry (this 
means of course that the "generators" defined as the integral curves of the degeneracy 
subspaces of the 3-metric coincide with the generators of the embedding of M) and 
X and uj are related to the spacetime metric and connection via 



X 



ollOlrX 

UJ + d log«L — d log Or, 



(97) 
(98) 



X = n L -n R 



(99) 



and 



uj = 



(100) 
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X thus measures the component of the spacetime metric not captured by the 3- 
metrics gi and gR of M at So, while ui is somewhat akin to an extrinsic curvature, 
since it measures an aspect of the derivatives along So of its two null normals, Ul and 
n R . 

Note that the sign of x is determined once the sense of and of ur are chosen. 
If both point away from So, or both toward So, then x must be negative, for then til 
and ur are both null future directed or both null past directed when j\f is embedded 
in a solution spacetime, so their inner product must be negative. 

The initial data x an d ui have been postulated to transform under a change of 
parametrization of the generators in the same way as the solution spacetime quantities 
that they correspond to, so the correspondence does not depend on the parametriza- 
tion chosen. 

An interesting alternative expression for uj in terms of the spacetime geometry can 
be obtained if one extends both m and ur to all of Af in such a way that both are 
orthogonal to the cross sections of M cutting the generators at a constant parameter 
value. Then 

uj = ' , (101) 

n L ■ n R 

where the Lie bracket, a vector, is mapped to a 1-form using the spacetime metric. 
(The right side is called the normalized twist |Epp95| but is not to be confused with 
the "twist" of a single congruence [Wald84] which always vanishes for the generators 
of a null hypersurface.) 

To see this suppose a tangent vector s of S is Lie dragged along ha to all of 
the branch Ha- This defines a vector field s orthogonal to both til and tir, since s 
is everywhere tangent to the isoparametric cross sections. The Lie derivative of the 
orthogonality relation s^ur^ = together with the fact that s satisfies £ nA s = 
implies 

= s»£ nA n Bli = 8"V nA n Bli + n Bfl V s n A - (102) 

Thus for all s tangent to So, uj ■ s takes the same value whether computed using 
the expression ( 11001) or the projection of the normalized twist f llOll) . proving their 
equivalence. 

From (11 01 1) it is clear that uj measures the non-integrability at So of the normal 
planes to the isoparametric sections of M. Integrability cannot in general be achieved 
by reparametrizing the generators because the transformation law (I98j) implies that 
the curl of uj on So is independent of the parametrization. 

Sachs [Sac62] and Rendall [Ren90j have given arguments that imply that any 
choice of these data (satisfying certain inequalities)!^] on j\f define a solution to Ein- 

51 The vacuum Einstein equations imply, via the Raychaudhuri focusing equation ( Wald84 eq. 
9.2.32), that if the cross sectional area of an infinitesimal bundle of neighbouring generators is 
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stein's field equations that is unique up to diffeomorphisms. More precisely, Sachs 
considered these data expressed in a special coordinate system adapted to the data - 
in which the generators lie along coordinate axes and each generator is parametrized 
by a coordinate which becomes an affine parameter in the corresponding solution 
spacetime - and argued that the solution metric exists and is unique in a particular 
coordinate system adapted to the solution which extends his coordinates on M to 
spacetime. Rendall then proved that Sachs' claim is true, at least in a neighbourhood 
of So in the causal domain of dependence D [Sf] of N '. 

The chief advantage of using coordinates adapted to the data, such as Sachs' coor- 
dinates, is that it makes it easier to identify freely specifiable data. The data g^, g^, 
X and uj are not themselves freely specifiable since gi and g^ are constrained to be 
degenerate, with restrictions on the degeneracy subspaces. Taking the generators as 
coordinate axes reduces this constraint to the requirement that certain metric compo- 
nents vanish (see (19~4")1 ). The use of a coordinate adapted to the data to parametrize 
the generators further reduces the freely specifiable data. Sachs normalizes the pa- 
rameters on the generators so that x = — 1 always. Moreover, the requirement that 
the parameters become affine parameters in the matching solution implies a constraint 
on the parameter dependence of the 3-metric on M . The cross sectional area of an 
infinitesimal bundle of neighbouring generators must satisfy the Raychaudhuri focus- 
ing equation corresponding to an affine parametrization ( [Wald84] eq. 9.2.32) with 
the Ricci curvature term set to zero, which substantially reduces the independently 
specifiable metric data. 

Another way to view the constraint on the metric components just mentioned 
is the following: Given any 3-metric on M satisfying the conditions of smoothness 
and degeneracy we have required of the initial data, and satisfying the inequalities 
explained in footnote EH there is a parametrization of each generator, unique up 
to an affine transformation, such that the metric satisfies the constraint. If this 
parameter, which I will term a "proto affine parameter" is adopted the metric must 
satisfy the constraint, but only because of the choice of parameter, and not because 
of Einstein's field equations. The field equations imply that in the matching solution 
the proto affine parameter becomes an affine parameter of the spacetime geometry. 
This argument will be developed in detail in subsection 13.51 

We will define both an almost freely specifiable data set and a completely free 
data set using charts similar to those of Sachs but with the proto affine parameter 
replaced by the area parameter r defined in (13.11) . The area parameter provides a 

decreasing in one direction along the bundle at one point then it must continue to decrease in the 
same direction until it reaches zero. That is to say,if the derivative of the area by any diffcrcntiablc 
parameter of the generators is negative at a point then it must remain negative until the area 
vanishes. This imposes inequalities on the 3-metric on J\f, although it does not imply any restriction 
on the variables Sachs used to describe the initial data. 
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good parametrization of the generators of a branch of M iff the cross sectional area 
of infinitesimal bundles of neighbouring generators in the branch is either everywhere 
increasing or everywhere decreasing toward the future. This is not the case in all 
solutions so our formalism cannot be applied to all initial data compatible with the 
field equations, but only to data satisfying this restriction on the 3-metric on 
However, as explained in subsection 12. 11 this leads to relatively little loss of generality 
in our results. The causality requirement implies that only data living on the same 
generator have non-zero brackets. In the calculation of the brackets on the generators 
through a point p G So we may therefore replace M by M' C M swept out by the gen- 
erators through a compact subset S' of Sq containing p in its interior. If the expansion 
rate of the congruence of generators of Ml and that of the generators of Mr are both 
non-zero at p (the generic case) then S' may be chosen so that these expansion rates 
are non-zero and of uniform sign on S' . If the expansion rate at S' of the generators 
of a branch is negative then by the field equations and the Raychaudhuri equation it 
will be negative on the whole branch, and r will decrease monotonically on the gen- 
erators. If the expansion rate is positive at S ; then r increases monotonically along 
each generator until it leaves the branch or until r reaches a maximum. In the latter 
case our results apply until the maximum is reached. Only if the expansion rate of 
the generators at p of one of the branches of TV is zero are our methods inapplicable. 

In the following subsection the coordinates that will be used to define our almost 
free and completely free initial data will be defined precisely. The relation between 
these data and Sachs' data is treated in detail in subsection 13.51 

3.3 Coordinate systems and free initial data 

Three special coordinate charts will be used to represent the different elements of the 
initial data. The clr chart, with coordinates a R = (ur^r^^i/^), covers Mr and a 
spacetime neighbourhood of intMR, is adapted to the generators of Mr, and is fixed 
in the manifold (i.e. solution independent) on Sr C dM. The ai chart, which will be 
defined in strict analogy with the aR chart, covers Ml and a spacetime neighbourhood 
oiintML, is adapted to the generators of Ml, and is fixed to the manifold on Sl C dM. 
Finally, the b chart, with coordinates fe M = (v + , v ~, 1 , 9 2 ), covers all of M and a 
spacetime neighbourhood of intM and is adapted to the generators of both branches 
of M, but is not a priori fixed anywhere. 

Greek lowercase indices a,(3, ... from the beginning of the alphabet will always refer 
to an a coordinate basis (of the branch Ma under consideration). Greek uppercase 
indices A, S... refer to the b coordinate basis. Lowercase Latin indices from 

52 This restriction is not a constraint because it implies only that the 3-metric satisfies inequalities, 
and not equations. 
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the latter part of the alphabet range over {1,2} and will be used to identify the 
two adapted coordinates y\ and y\ and their corresponding coordinate basis vectors, 
while lower case Latin indices a,b, ... from the beginning of the alphabet, which also 
range over {1,2}, identify the coordinates 8 1 and 9 2 and their basis vectors. 

The chart a A = (ua, ^a, y\, Ua) adapted to Ma is defined in three steps. First 
smooth, fixed (not solution dependent), coordinates y\ and y\ are chosen on Sa- 
Next this coordinate system is extended to a chart on all of Ma by convecting the 
y\ along the generators, i.e. setting them to be constant along the generators, and 
introducing the coordinate ta, taken to be the area parameter along the generators 
defined in subsection 13.11 Thus on Ma the two y coordinates identify the generators 
and ta identifies the points on each generator. Finally coordinates are defined on a 
spacetime domain containing intMA in its interior: Each constant ta section of M 
is a two dimensional spacelike disk and thus posses two null normals at each point, 
and two corresponding normal null geodesies, One is the generator of Ma through 
the point and the other is transverse to Ma- To extend the coordinate system to 
a neighbourhood of Ma the three coordinates y\, y A , and r A are convected along 
the transverse normal null geodesies and a fourth coordinate, ua, which parametrizes 
these geodesies is defined by the requirements Ua = on Ma and d UA ■ d rA = — 1. 

To lighten the notation the subscript ("A") identifying the branch Ma to which 
the coordinates are adapted will be dropped from here on when there is little risk 
of confusion. 

The components of the metric on a branch Ma take a very restricted form in the a 
coordinates adapted to that branch. The coordinate basis vectors d r and d u are null, 
and have inner product d r ■ d u — — 1. Moreover both are orthogonal to the constant 
(u, r) surfaces. Thus the line element takes the very simple form 



on Ma, where hij is a symmetric 2x2 matrix which must be positive definite in 
order that the spacetime metric be Lorentzian with signature — h + + . The induced 
3-metric on Ma is seen to be degenerate and takes precisely the form (194j) obtained 
in subsection 13.11 The line element on any 2-surface cutting a cross section of the 
congruence of generators is ds 2 = hijdy % dyK This means in particular that area 
density of the cross section is, 



the root of the determinant of the induced 2-metric. (As expected this depends only 
on position in Ma, and not on the tangent plane of the cross section.) The value of 
the coordinate r, the area parameter, at a point p e Ma is thus 



ds 2 = —2dudr + h i jdy t dy : > 



(103) 




(104) 




(105) 
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where p(j) is the value of p at the endpoint on Sa of the generator 7 through p. 
Equivalently 

p(r,y l ) = r 2 p(y i ). (106) 

The only freedom in p is thus its boundary value p on Sa- 

The remaining freedom in hij resides in the conformal^l 2-metric, a unimodular 
(determinant equals 1) 2 x 2 matrix defined by 

e i j = h ij /p. (107) 

Information about the metric geometry of Ma is also found in the range of values 
the a coordinates take on Ma- (y 1 ,!/ 2 ) ranges over a fixed (solution independent) 
domain diffeomorphic to a disk in IR 2 . r on the, other hand, ranges on each generator 
from 1 on Sa to a value r (y) on S (depending on the generator). The function 
r (y) depends on the solution metric, since a variation of the metric will in general 
change the rate at which the generators are focused, contracting or expanding cross 
sections of bundles of neighbouring generators, and thus affects tq. [ro may be either 
greater than 1, (generators converging toward the future) or less than 1 (generators 
diverging), and though r typically varies over S , it never passes through 1 because 
of the restrictions that have been imposed on the initial data to ensure that the area 
parameter is good on every generator.] 

Further information is contained in the map <j) from the coordinates yn to the co- 
ordinates ul on So- A perturbation of the solution metric, for instance a gravitational 
wave pulse passing through Mr but not Ml, will in general affect which generator 
from Sr intersects a given generator from Sl, and thus which values yR and of yi 
correspond to the same point on So- The map <fi is thus solution dependent. 

and the functions r (y), p(y), and eij(r,y) on each branch determine the 3- 
metric on M modulo the action of a group T>°^ of homeomorphisms M — > M which 
map each branch diffeomorphically to itself and reduce to the identity on Sr and 
Sr.ri The functions r M (|/i?), ^lo(i/l) an d 4> alone, without further specification of 
the 3-metric on M, determine the cir and coordinates on M up to the action of 
T>°j^. To see this consider the transformation of the a charts due to a change of the 
metric, within the class of metrics we allow on M, which preserves rLo(yii), tlo (vl) 
and <p. (Recall that once the y coordinates are fixed on Sl and Sr the a charts 
are completely determined by the metric, so this type of transformation is the only 
freedom that remains in these charts once t^o^r), r^o^i) and <f) are given.) The 

53 It is invariant under rescaling of the metric, hence the name. 

54 These homeomorphisms could perhaps be called diffcomorphisms since they respect all the 
natural differcntiable structure of Af. Furthermore, by the Whitney extension theorem the homeo- 
morphisms in Dj[f can always be extended to diffcomorphisms in a spacetime neighbourhood of Af 
taking each branch of Af to itself and leaving Sl and Sr invariant. 
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displacement of the a coordinate points induced on each branch by this change of 
metric leave Sl and S R invariant, since the coordinates are fixed there. Furthermore, 
the invariance of ro(y) on each branch (and the smoothness of the metric) implies 
that these displacements define a diffeomorphism of the branch to itself, with So 
mapped to So- Finally, the invariance of <ft implies that the diffeomorphisms on the 
two branches induce the same diffeomorphism on So, assuring that the transformation 
of the cil and a R charts is an element of T>°^-. 

If in addition to r^iva), fLoiui) and the non-zero a components of the 3-metric 
(i.e. hij = r 2 p(y)eij(r,y)) are given on each branch, then the 3-metric is determined 
up to homeomorphisms in T>°^. This is the most one could expect, since all these 
quantities are invariant under the action of T>°^ on the 3-metric. 

In the calculation of the brackets we shall use essentially these functions to rep- 
resent the 3-metric part of the initial data. In particular the conformal 2-metric on 
M, specified by the functions e Ri j(r R ,y R ) and eLij(r L ,y L ), will represent the bulk of 
the initial data. The remaining data will ultimately be represented by fields living on 
S - fields equivalent to </>, r^o, and pa, and x an d to corresponding to a particular 
parametrization of the generators. 

These quantities are not entirely unconstrained. We require that the induced 3- 
metric be consistently defined at So in the sense that the same 2-metric is induced on 
So from the 3-metrics of either branch. Taking into account the transformation from 
Vl to y R coordinates this leads to the conditions 



PR r m = \ det g^\P Lr Lo ( 108 ) 
and 

dyfdyl dy R 
eRij = -^-TTTTldet- — \e Lmn , (109) 
dy R dy R dy L 

and of course is required to be symmetric and unimodular. 

The scheme of two overlapping adapted charts with solution dependent transfor- 
mations between them is rather complicated. Why not base our adapted coordinates 
at S , instead of S R or S^? The coordinates 6 M = (v + , v~, Q 1 , 6 2 ) define such a chart: 
9 1 and 9 2 are coordinates on So which are convected along the generators to all of 
M. v~ and v + are area parameters on the generators of M R and Ml respectively that 
are normalized to 1 on So- v + , 9 1 , and 6 2 are extended to a spacetime region by 
convecting them from Ml along the congruence of null geodesies that are normal to 
the constant v + sections of Ml and transverse to Ml, and v~ is similarly extended by 
convecting it in the same way from Mr. 

The b coordinates will be used to represent data living on So, and to facilitate 
comparison with Sachs' initial data. It is tempting to express all the data in terms of 
these coordinates. In fact the brackets will ultimately be stated for data referred to 
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this chart. However we will not use such data to calculate the brackets. To solve (1531) 
for the auxiliary pre-Poisson bracket we need to be able to identify those variations of 
the metric that vanish (in a fixed chart) in a neighbourhood of dM from their actions 
on the initial data. But even variations that vanish near dM do in general deflect the 
generator of Ma emerging from a given point of So, changing its course all the way 
to Sa, and thus would produce a variation in the b coordinates at Therefore 
data that are local functions of the metric components and their derivatives in the b 
chart generally vary at dM under such variations, complicating their identification. 
It would be simpler if the variations that vanish near dM also left the initial data 
associated with points near dM invariant. 

This is the case for data constructed from the a coordinate components of the 
metric. The variations that leave the metric invariant near dM also leave the a coor- 
dinate components of the metric invariant near dM, because such a variation leaves 
the a coordinates themselves invariant in a neighbourhood of dM: The coordinates da 
are fixed on Sa, so a variation 5 that vanishes in a spacetime neighbourhood W of dM 
clearly leaves the generator corresponding to given y\ invariant in a neighbourhood of 
Sa- Furthermore any generator originating within a sufficiently small neighbourhood 
U of dSo in So remains inside W until it leaves M (See prop. IB. 81) . Such generators 



are entirely undisturbed by 5. It follows that there exists a neighbourhood V of dM in 
M, possibly smaller than W C\M, such that <5a^ = on V hMa- The construction of 
the a a coordinates off Ma then implies that these are in fact invariant in a spacetime 
neighbourhood of cWa H Ma- 

The chief condition defining the bracket on initial data is that the equation (1531) . 



6<p = n[{<p,-} m ,6\, (110) 

holds for all initial data ip and all variations 5 G LP - which is to say all 5 satisfying 
the linearized field equations and vanishing in a spacetime neighbourhood of dM. 
Under such variations 

be-Aij = in a neighbourhood of dM (HI) 

and 

5p A = 0. (112) 



55 The b coordinates can be held fixed near the end Sa of one branch Ma under variations of 
the metric vanishing near the boundary by giving the 9 coordinates fixed values on Sa and then 
convecting them to Sq along the generators, and from there to the other branch. The 9 would then 
be given by a solution independent function of the ua coordinates, or at least by one that is invariant 
under the variations being considered. However, since it is easily shown that the map from the yi, 
to the ur coordinates varies under variations of the solution, even variations that vanish near dM, 
these 9 cannot be invariant at the end of the other branch. Thus there is no choice of 9 which makes 
the b coordinates invariant near all of dM under variations of the metric that vanish there. 
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I shall solve ( 11 101) for the auxiliary bracket imposing only the restrictions (11111) and 
(11121) on the variations 5. This is a priori a weaker restriction than 5 G which 
requires that 5 vanishes in a whole spacetime neighbourhood of dM. But certainly if 
(11101) holds for our larger class of 5s then a fortiori it holds for 5 G L° g . 

On M the a and b coordinates are related in a fairly simple manner. On Mr the 
coordinates 6 M = (v + , t>~, 1 , # 2 ) can be expressed in terms of the a# chart as 

^+ = 1 v- = r R /r Ba {y R ) 9 a = f a R (y R ), (113) 

where /r is a diffeomorphism from H 2 to H 2 . Similarly on Ml the b A are related to 
the al via 

^ + = r L /r L0 (?/L) «" = 1 e a = f L {y L ). (H4) 
Of course Jr = Jl ° <P where <fi is the (solution dependent) diffeomorphism from yR 
to y L coordinates on So- 

One could take to be the identity map (i.e. 9 = y L ) and f R = </>. However it 
will be more convenient not to assume any particular relationship between the 9 and y 
charts. It will be assumed that the 9 coordinates of a point in So depends only on the 
yL and y R coordinates of that point and possibly on the initial data on the generators 
through the point. That way the b coordinates, like the a coordinates, move along 
with the initial data when these are acted on by a homeomorphism in T>^-, and the 
components of the data in the b chart are invariant under such homeomorphisms. 

The parameters v + and v~ on the generators of Ml and Mr respectively define 
tangent vectors to these generators: 

n_ = d v - on Mr (115) 
n + = d v + on Ml (116) 

and define corresponding initial data \ v and uo v which, in a solution matching the 
initial data, are given by 

Xv = n + -n_ (117) 

n+ ■ V a n_ - n_ • V a n+ 

u va = . (118) 

n + ■ n_ 

Another datum, which will be used in place of pL and pr, is the area density on 
So in 9 coordinates: 

Po = r 2 R0 p yR \ det[^)\ =r 2 L0 p yL \det[^-]\. (119) 

The two data pL and pR can be replaced by a single datum because of the constraint 
( 11081) between them, which is expressed by the second equality of (I119p . 

The preceding considerations show that a possible initial data set for general 
relativity consists of 
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• e Lij (r L ,y L ) on Af L and e Rij (r R , y R ) on Af R , 



and on S< 



o 



• r L0 (y L ) and r R0 (y R ) 



• the maps /l : j/l -> and f R :y R ^ 6, 



• Po(0), 



• x„(0), 



These will called the boundary fixed coordinate data (BFC data) since they are 
referred to charts that are fixed in the manifold on portions of the boundary dM of 
the initial data hypersurface under the variations of the solution being considered!^! 
They are subject to the one constraint fllU9p requiring the consistency of the definition 
of the conformal 2-metric at So (and of course e is also required to be symmetric and 
unimodular). 

When using the BFC data we shall drop the subscript v on \ and u. x an d oj 
corresponding to other parametrizations of the generators will have no role to play. It 
will also turn out to be convenient to use instead of fi and f R their inverses, sl and 
s R . Finally, Xi will often be represented by A = — log ] x I ■ T ne sig n °f X is uniform 
on S for the solutions and the M we admit, and does not change under variations, 
so taking the absolute value represents no real loss of information. 



The domains of sensitivity 57 ! of all these data are such that the causality condition 
on the auxiliary pre-Poisson bracket implies that data living on different generators 
have vanishing bracket. These domains are confined either to only one generator or to 
the pair of generators emerging from a single point of Sq\ The domain of sensitivity of 
e at a point p on a generator 7 is the segment of 7 from p to <97V, that of r^oata point 
Po G So is the whole generator of Ma originating at po, that of a (y R ) = f R (y R ) is by 
assumption confined to the pair of generators originating in the point of So specified 
by 9. The domains of sensitivity of all the remaining data at the same values of 9 a 
are also contained in the union of these generators. 

We shall calculate the brackets between the BFC data, but ultimately the brackets 
will be expressed in terms of a completely free data set consisting of 

56 Since the variations of the metric under consideration do not affect generators originating 
sufficiently close to dSo even the b chart is fixed in a sufficiently small neighbourhood of OSq in Sq, 
which is where the data referred to the b chart live. 

57 Recall that we have defined the domain of sensitivity a functional F of the spacetime metric 
to be the support of the functional gradient SF/Sg^ of F. Thus perturbations of the metric in this 
domain, and only in this domain, affect the value of F to first order. 
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• e a b(v,9), the conformal 2-metric on constant v + , v surfaces in b coordinates, 
specified on all of M 

and 

• r L0 , r R0 , s L , s R , p , A, u a specified on S . 

The data r^o and J'a (or sa) define the transformation from the a a chart to the b 
chart on Ma, thus once these data are given specifying e ah is equivalent to specifying 
CAij- Recall that because of the degeneracy of the 3- metric of M the induced metric 
at a point p on a two dimensional spacelike cross section of M depends only on p G A/", 
and not on the tangent plane of the cross section. Thus 

d6 

e a b = d a y A dby 3 A \ det - — \e Aij (120) 
oy A 

in complete analogy with (11091 ) . 

The only constraint on the BFC data, (11091) . simply requires that e a b obtained 
from e^, r L0 and sl agrees with e a b obtained from e R , r R0 and s R at S . It does not 
limit e a b itself. For this reason we shall refer to this new data set as the completely 
free (CF) data. 

The data tao, will turn out to play a somewhat secondary role. Their variations 
can be eliminated from the presymplectic 2-form. This seems natural as they essen- 
tially specify where the generators of the branches of M have been cut off, something 
that we are free to change without changing the solution. 

3.4 Variations of fields in the a and b charts 

The notion of variations of fields, i.e. the derivative of the fields along a one parame- 
ter family of field configurations, is fundamental to the variational principle, and the 
symplectic potential and presymplectic 2-form constructed from it. We have defined 
the variation of a field on a manifold in terms of a chart fixed to the manifold: The 
components of the variation in the fixed chart are the variations of the correspond- 
ing components of the field in this chart. The transformation of the components of 
the variation from one fixed chart to another is straightforward since the map from 
the old to the new chart is field independent. Thus the variations of tensor compo- 
nents transform as do the original tensor components, defining a tensor field, and the 
variation of a connection also forms a tensor field. 

The expression (I39p for the presymplectic 2-form therefore has a straightforward 
interpretation if it is supposed that the field components are referred to a chart fixed 
in the manifold. However we wish to express the presymplectic 2-form in terms of the 
BFC data, which consists of metric and connection components in field dependent 
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coordinate systems. How are variations in such a moving chart related to variations 
in a fixed chart? It turns out that they differ by a Lie derivative. 

Consider first the simplest case - the variation of a scalar field <p. It is clear that 
the variation 6 w tp of the field at a given grid point of a chart w, is related to its 
variation 5 F ip at the corresponding grid point of a coordinate system x M that is fixed 
in the manifold via 

S w ip = 5 F ip + 5 w x t *d ll ip. (121) 

The variations of the x coordinates at aw chart grid point defines a vector v = 5 w x fl d fl 
which may be thought of as the variation of the position of the w grid point in the 
manifold. Thus 

<TV = 5 F cp + d v ip. (122) 

A similar relation holds for the variations of tensor fields, such as the metric g, in the 
w chart: 

5 w g = 5 F g + £ v g. (123) 

Here 5 w g is the tensor field with w components equal to the variation 5 of the w 
components of g at fixed values of the w coordinates. That is 

[S w g] w = 6[g] w (124) 

where [g\ w indicates the w components of g expressed as functions of the w coordi- 
nates. Notice that if w' is a chart which moves together with w, i.e. the mapping 
between the charts is field independent, then 5 W = 5 W on tensors. In other words, 
what is important is not the chart but how it moves when the field is changed. 

These results generalize to any field X that has well defined components in any 
C°° chart. In the possibly field dependent chart w the components of S W X are just 
the variations of the w components [X] w of X. To define the components of S W X in 
another we impose the requirement that, just as for tensors, the variation 5 U X in a co 
moving chart u be equal to 5 W X. Specifically, suppose wo is the w chart corresponding 
to the fiducial field configuration we are varying about, and Uq is another chart. Then 
Uq o Wq 1 is the (field independent) diffeomorphism that maps the w chart to the u 
chart, u = Uq o Wq 1 o w is a chart co-moving with w that coincides with u at the 
fiducial field configuration, and 

[5 W X] U0 = [6 U X} U0 = 5[X] U . (125) 

With this definition 5 W = 5 F + £ v quite generally!^ Suppose 5 is the tangent d/da 
to a family of configurations of the gravitational field parametrized by a e IR, with 

58 The general formalism is presented mostly for the sake of completeness. The reader may find it 
easier to deduce the form of 5 W X for each field X that is needed directly from the definition (|125p . 
The fields of interest, aside from four dimensional tensors, will be Christoffcl symbols and the BFC 
data. 
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a = at the fiducial solution. Let X a be the corresponding family of X fields, w a be 
the corresponding w charts, <j) a the family of diffeomorphisms mapping wq to w a , v 
the vector field generating <p a , and 0° the mapping induced by <p a on the components 
of X, i.e. [X] Wa = ^[[XU]. Then 

[8 W X] W0 = ^-[X a ] Wa \ a=0 (126) 

= ^i[[x a u\ a=0 = rxu + ^ o a {[XoU\a= = r + £ v \xum 

The equation d(j)° a [Xo] m /da\ a = = [£ v X] m is the natural generalization of the defi- 
nition of the Lie derivative for tensors (see [Wald84] Appendix C, but note that <ft* 
there is the inverse of our <p°). It defines £ v to be the generator of the action of the 
family of diffeomorphisms <p a on X. 

The relation between the variations of components in differently moving charts is 
entirely analogous to that between variations in fixed and moving charts. For instance 
the variations of a and b components are related by 

5 b = 5 a + £{, (128) 

where £ is the difference of the v fields associated with the a and b charts. £ can be 
evaluated in the a chart by applying (I128P to the a coordinates themselves: 

5 b a a = + £^a a = £ Q , (129) 

and in the b chart by applying it to the b coordinates: 

To lighten the notation the superscript indicating the chart relative to which a 
variation is taken will often be dropped. When no chart is indicated the variation 
of a field is to be taken with respect to its "natural chart". The natural chart of a 
field component (of a tensor or any other multi-component field) is the chart that the 
component is referred to. For instance 5g^, with the components of the metric 
in a fixed chart, is S F g fllJ . The BFC data, and the variants we have discussed, each 
have a corresponding natural chart: eAij, r^o, and f% all have natural chart a a- Po, 
X, uj a , and s l A have natural chart b. Note that s\ is not an cla vector components, but 
a scalar function of the 9 a . Similarly f% is a scalar function of the ua- The natural 
chart of e a b is of course b. Finally, we adopt the convention that the natural chart of 
a coordinate is the chart it forms part of. Thus 5a a = 5 a a a = and 5b A = 5 b b A = 0. 

Although the b chart is clearly a moving chart, being adapted to the generators of 
A/", one may, when solving (1531) for the auxiliary pre-Poisson bracket, restrict attention 
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to variations of the solution that leave the fixed in a neighbourhood of So, where 
the data referred to the b chart live: (1531) requires that for any initial datum tp on the 
interior of M 

se = ti[{<p,-}„6\ (i3i) 

for all 5 which satisfy the linearized field equations and vanish in a neighbourhood of 
dM. The right side of this equation is unchanged by the addition to either argument 
of Q of a diffeomorphism generator £ T vanishing near dM since such a generator 
is a degeneracy vector of Q. Moreover, if tp is a BFC datum the left side is also 
invariant under the addition of such a term to 5. Now, if tp is a BFC datum living 
at a point p G intM then {tp, •}* acts only on the data on the generator through p, 
thus leaving the b^ invariant in a neighbourhood of dM — Sl — Sr. As explained in 
subsection 13. 3[ 5 also leaves the b^ invariant in a neighbourhood of dM — Sl — Sr. 
The (linearized) diffeomorphisms induced by {tp, •}* and S can thus be cancelled near 
Sq by diffeomorphism generators that vanish in a neighbourhood of dM. 

Note that the pair of charts cir and cannot both be held fixed on So in the 
same manner as the b chart because the transformation from the clr chart to the 
chart depends on the fields. 

3.5 Freeness and completeness of the BFC and CF data 

To show that the BFC data of section [3] is free (apart from the constraint (11091) ) and 
complete it is sufficient to show that any BFC data determines a set of Sachs data 
such that a solution matches the CF data iff it matches the Sachs data. Since we 
assume that any Sachs data determines a unique solution, modulo diffeomorphisms, 
the same would also be true for the BFC data. Furthermore, since the CF data are 
equivalent to the BFC data, this would establish that the CF data are also complete 
and completely free. 

The BFC data, the CF data, and the Sachs data are all explicit coordinate repre- 
sentations of the geometrical initial data defined in subsection 13.21 - the 3-metric on 
M and x an d oj on Sq - but referred to different coordinate systems. The BFC data 
is 

in part referred to the b chart and in part to the a charts, the CF data refers 
to the b chart, and the Sachs data refers to essentially the b chart, but with affine 
parameters of the generators in the solution spacetime replacing the area parameters 
ir^o The key point to be demonstrated is thus that the different data sets contain 

59 Sachs also uses a special chart, adapted to the geometry, on So - essentially the ip chart defined 
in subsection [231 But this specialization of the chart on So plays no role in the results of Sachs and 
Rcndall on the existence and uniqueness of solutions, so we may as well use the 9 chart on So to 
define Sachs' data. 
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sufficient information to define the map to the coordinate system of the other data 
sets. 

The BFC and the CF data each define the transformation (II 131)1 141) from the a 
chart on a branch of M to the b chart on that branch. In particular they define the 
map from the coordinates y l on each branch to the coordinates 9 a . This allows one to 
compute the CF data corresponding to a spacetime metric from the BFC data, and 
vice versa, establishing the equivalence of the these two data sets. 

It remains to show that the BFC data define the map from the b chart to Sachs' 
coordinates, in particular the map from the area parameter r to an affine parameter 
7] (and thus also from v = r/r to rj). In a non degenerate Lorentzian solution 
matching the BFC data the generators of M are geodesies: The degeneracy vectors of 
the degenerate 3-metric defined by the BFC data are orthogonal to all tangents to M 
and null, so M is a null hypersurface and the degeneracy vectors are its normals - from 
which it follows that the integral curves of the degeneracy vectors are null geodesies 
(see prop. IB.3I of appendix [B]) . Now if 77 is an affine parameter of a generator then 
the corresponding tangent vector k = d v is parallel transported along the generator: 
Vkk = 0. In a coordinates this equation is equivalent to 

But on M the field equation G rr = is equivalent to 

r rr = -Ld^dretj. (133) 
o 

(see Appendix 0) Combining these equations one obtains the relation 

which may be integrated using BFC data to obtain 77 as a function of r up to an 
additive constant and a constant factor on each generator in concordance with the 
freedom to change affine parameters by affine transformations. (The "constants" are 
independent of r, but may depend on y.) Of course 77 obtained in this way is not 
actually an affine parameter. Only when the field equation G„ = holds does (11 34ft 
imply f l 1 3 2 f) . making rj is an affine parameter of the spacetime geometry. 

Sachs sets the affine parameters, r] L and r/ R , of the two congruences of generators to 
zero on Sq. To partly fix the normalization of the affine parameters he also demands 
that d VL ■ d VR = — 1 on So- The remaining freedom in the normalization of the 
generators can be fixed in our case by requiring 

dr} R /dv- = 1 (135) 
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on S . Then, by JHZj), 

dt] L /dv + = -Xv (136) 

The BFC data determine the a coordinate components of the 3-metric on Af, and 
also, as we have just seen, the transformation to Sachs' coordinates on M . They 
therefore determine the 3-metric of j\f in Sachs' coordinates. This determines most 
of the Sachs data. 

The Sachs data consists of the conformal 2-metric, e a b(r), 6), as a function of Sachs' 
coordinates on each branch of M and four fields on Sq, referred to the coordinates 9 a : 
the area density p(6) and its derivatives d VL p and d VR p, and ou v , the normalized twist 
uj corresponding to the parametrizations tjl and ri r of the generators]^] All but the 
last of these are subsumed in the 3-metric of A/I 61 ! PI 

The final Sachs datum, u v , is related to, u v , the normalized twist corresponding 
to the parametrizations v + and v~ of the generators via 

to ria (9)=u va (9)-d a \og X v(0). (140) 

(See (EHD, ([136)) and (TT351) .) 

Any solution matching the BFC data necessarily also matches these Sachs data. 
Since the solution is uniquely determined by the Sachs data it is uniquely determined 
by the BFC data@ Thus the BFC data are complete. 



60 Sachs actually takes as his final data a pair of quantities he writes as Ca.i A = 1,2. However 
these are just the 9 coordinate components of the 1-form — Lo n . This can be seen most easily from 
his equation (19), which when a forgotten factor of 1/2 is restored and it is rewritten in our notation 
reads 

\c a>1 =d r , L -V a d VR . (137) 

This equation, the normalization d VL ■ d riR = —1 on So, and the expression (|100p imply the claim. 

61 The explicit relations of these data with the BFC data are not complicated. Since the 3-metric 
on j\f is degenerate the induced 2-metric at a point on a constant -q section of Af is the same as 
the induced 2-metric on the constant r section at the same point. Thus the e a b at a point are 
determined by the e^j at that point and the standard tensor density transformation associated with 
the transformation from y to 6 coordinates. p(9) is both a Sachs and BFC datum, and the derivatives 
of p are fixed by the normalization conditions Q136p and (|135p on the rj parameters: 

d VL p = 2pdv + /d? ]L = -2p/ X v (138) 
d VR p = 2pdv-/dr lR = 2p. (139) 



62 In the coordinate independent description of the data in subsection 13.21 a further datum, x, is 
specified on Sq, but this is identically equal to —1 for the parametrization r/L, r]R of the generators. 

63 It has been proved that the Sachs data determine the spacetime metric uniquely up to diffeo- 
morphisms in a neighbourhood of Sq in 22 |Ren90j , and we are assuming that they in fact determine 
a maximal Cauchy development of j\f. This implies that the components of the metric in Sachs' 



61 



The BFC data are also free, aside from the one constraint (11091) . This follows 
from the fact that, firstly, any BFC data that satisfies (11091) defines Sachs data, and 
secondly, that the solution corresponding to this Sachs data really matches all the 
original BFC data. To verify the second claim consider first the area parameter r', 
defined like r but from the 3-metric induced on M by the solution, r', like r, must 
satisfy (11341) which may be re-expressed as 

5=™<i). d«) 

with k = ^drje^ d v eij a function of r\ that is determined by the Sachs data and is 
thus the same for the original 3-metric as for the 3-metric induced by the solution. 
The solution to this linear second order equation is completely determined by the 
value of d\ogrjdr\ at So and the value of r at the other end of the generator. At Sa 
r' A = ta = 1 by definition, and at So d\ogrA/dr]A = d\ogr' A / ' dr\A = l/2d\og p/dr/A 
since this is part of the Sachs data defined by the original BFC data. Thus r' A = ta, 
and, since r(r]) and the Sachs data determine the 3-metric completely (without use 
of the field equations), the 3-metric induced from the solution is equal to the original 
3-metric. The 3-metric BFC data induced by the solution are thus the same as the 
original ones. Since v' = r'/r' Q = v it is also clear that the remaining BFC data, \ v 
and u> v , induced by the solution are equal to the ones originally specified. 



4 The presymplectic 2-form 
4.1 Preliminaries 

The auxiliary pre-Poisson bracket {•,•}* on functions of the initial data is defined 
by the requirements that it should reproduce the Peierls bracket on observables, be 
antisymmetric, that it be causal in the sense that the bracket vanishes between data 
at points that cannot be connected by any causal curve, and that it respects any 
constraints that the data must satisfy. As shown in section [21 {■, ■}* reproduces the 
Peierls bracket on observables iff for any initial datum 9 on intM 

5 9 = ti[{9,-}*,5 ] (142) 

for all 5q which vanish in a neighbourhood of ON '. 

spacetime coordinates are uniquely determined. Sachs' spacetime coordinates are obtained from the 
b coordinates by replacing the coordinates v^" by tjl and r]n, everywhere given by the same functions 
of v + and v~ respectively as on J\f. Thus the Sachs data define the solution in b coordinates uniquely 
as well. 
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In this and the following section Qf^i,^] will be evaluated explicitly in terms 
of the BFC initial data defined in subsection 13.31 for variations 5i and 62 in a class 
that is wide enough to include all 5o allowed in (11421) and a solution {9, to (11421) . 
Specifically Q[Si, 62] will be calculated for 5i and 82 satisfying the following conditions 
on each branch Ma of the initial data hypersurface: 

a) The variation v\ = —5 F a ( \ of the adapted coordinates a a A corresponding to Ma 
vanishes in a neighbourhood of dM H Ma- 

b) The variation 8pA of the area density pA on the truncating 2-surface Sa (where 
the generators of Ma cross the boundary dM) vanishes. 

These conditions are clearly satisfied by the variations 5q in (1142j) because they 
vanish near dM. On the other hand, the variation {9, ■}* generated by an initial 
datum 9 is not in general expected to vanish on all of dM. However, recall that (I142|) 
determines {8, only up to the addition of an arbitrary diffeomorphism generator 
£ v (a generalized degeneracy vector of Q). It turns out that for any bracket satisfying 
(11421) . and the causality condition, there exist suitable diffeomorphism generators such 
that the modified bracket obtained by adding these to the original bracket satisfies 
a) and b). 

Let us verify this in detail. The requirement that the auxiliary pre-Poisson bracket 
be causal implies that the variation {<fr(p), •}* generated by the value of an initial data 
field 4> at the point p £ M can be non-zero only on the generators passing through 
p. Since all our initial data live on the interior of M this implies that the variation 
generated by any initial datum vanishes in some neighbourhood of dM — Sr — Sl- 
But it would seem that the conditions a) and b) could be violated at some points in 
the interior of Sl and/or Sr. 

In order to work with smooth variations let us smear the initial data fields with 
smooth test functions. Thus, instead of taking 9 to be the value of an initial data 
field at an interior point of M we take 9 to be the integral, weighted by a smooth 
test function with support contained in intM, of such an initial data field. Imposing 
( 11421) for all such test functions implies that it holds distributionally for the initial 
data themselves, which is all we will ultimately require. 

Now suppose that {9, ■}[ is a solution to (1142ft . By hypothesis, the area density 
p transverse to a generator always increases or decreases along the generator - it 
is never stationary - in the spacetime geometry and for the M we are considering. 
Therefore if an infinitesimal variation of the geometry changes the area density at 
a point this change can always be cancelled by a suitable infinitesimal displacement 
along the generator. By adding the diffeomorphism generator £ v to {9, with rj 
taken parallel to the generators and with 
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magnitude such that d v p = — {0, p}* on Sl and Sr one may ensure that the area 
density on these surfaces is invariant under the combined variation. That is, condition 
b) is satisfied. 

What about condition a)? The modified action {9, + £ rj of 9 will in general 
affect the generator passing through a given (interior) point of Sa- Since the generator 
is identified by its point of intersection with Sa this point remains fixed by definition, 
but all other points on the generator will in general be displaced, inducing non- 
zero changes in the adapted coordinates a a arbitrarily close to Sa in violation of 
requirement a). These changes in the a a can clearly be undone in neighbourhoods of 
Sl and Sr by adding a second diffeomorphism generator, £ v i, with 7/ chosen to vanish 
on Sl and Sr so it does not affect p and b) remains satisfied. Furthermore, the original 
variation of the generators must vanish outside the causal domain of influence of the 
support of the smearing function of 9, and thus in a neighbourhood of ON — Sl — Sr, 
which implies that the action of 9, with the two added diffeomorphism generators, 
leaves the adapted coordinates invariant in a neighbourhood of all of dj\f. That is, 
{9, -}^ + £ 1) + £ v i is a solution to (11421) satisfying conditions a) and b) above. We 
may therefore solve ( 1142ft for {9, assuming a) and b) hold (and then recover all 
the remaining solutions by adding diffeomorphism generators to the first solution). 

One consequence of condition b) is worth noting. r 0L , r 0R , and p are a priori 
independent BFC data but, because Pa(u) — Po(fA(y)) r o(y)~ 2 det[<9j/^] is invariant 
under the class of variations being considered, the variations 5vaq m ay always be 
reduced to variations of other BFC data. We shall do this and variations of r$A will 
not appear in the expression for the presymplectic 2-form finally used to solve for 
the auxiliary bracket between the data. This does not mean that the brackets of the 
r Q A are undetermined. Since the variation {9, •}* satisfies b) by the very argument 
we have just given {9,Tao}* can be expressed in terms of the the brackets of 9 with 
other data. It is not difficult to do so explicitly and the reader may do so using the 
brackets between the remaining BFC data given further on. 

In the next subsection the symplectic potential Q[8] will be evaluated in terms of 
the BFC data for 5 satisfying conditions a) and b). Then, in the following subsection, 
the expression for obtained will be used to compute Cl in terms of the same data for 
variations satisfying a) and b). The solution of the resulting explicit form of equation 
( 11421) . which yields the auxiliary pre-Poisson brackets between the initial data, is left 
to section [5j 
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4.2 Calculation of the symplectic potential in terms of the 
almost free BFC initial data 

We shall calculate explicitly only the part of the symplectic potential corresponding 
to the branch Mr, denoted Or. Since the two branches of M enter the problem on 
an exactly equal footing, Ml contributes an entirely analogous expression. 

Or is just the contribution of Mr to the boundary term (I19p in the variation of 
the action integral (on a region of which Mr forms part of the boundary). That is 

e*0 = I r ^ta^V^g (i«) 

o7TLt JNr 

in terms of field components and their variations referred to manifold fixed coordi- 
nates. 

Recall that the BFC data on Mr consist of bulk data, namely the conformal 2- 
metric components in the a coordinates, specified on all of Mr, and surface data 
associated with the 2-surface Sq. The contributions to Qr[5] of the variations of the 
bulk data and of the surface data on Sq can be neatly separated by expressing 
the variation of the metric as the sum of a contribution due to the variation of its 
components with respect to the a coordinates corresponding to to the R branch, and 
another due to the variation of these a coordinates themselves (with respect to fixed 
coordinates). 

That is, we decompose Or according to 

eR[S}=e R [6 a }-QR[£ v }, (144) 

where v is the vector field —6 F a a d a . The first term can be written immediately in 
terms of a coordinate components, and reduces to an integral over Mr of a functional 
of Cij and its variation. The second term turns out to be a surface integral over So 
and depends on the data on So and their variations (as well as on the derivative of 
eij along the generators at Sq). 
In direct analogy with (1791) 

Q R [£v] = -77^ / Wv^G^V (145) 

16nG JdN R 

By condition a) on the variation we are considering v vanishes in a neighbourhood of 
dM R - Sq — dM n Mr, so 

Qr[£v] = 77^77 / WV^tEV (146) 
(The sign reflects the orientation of Sq which is taken to be opposite to that of 8Mr.) 
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As discussed in subsection 13.41 one may restrict attention to variations such that 
the b chart is fixed in the manifold in a neighbourhood of Sq in M. Then v = — £ 
on Sq, where £ is defined by S a = 5 b — That the substitution of v by — £ does 
not affect the symplectic potential may be verified directly using (11461) : The entire 
symplectic potential, obtained by summing the contributions of the two branches of 
M, is 

e[6] = e R [s a «] + e L [6 a >] + / v> L - vtfV^gdz,*,. (147) 

107rG J So 

As in (11461) the orientation of Sq is opposite to that of BMr, and thus equal to that 
of ONl- Now since 5 b = 5 + £ Vb and S aA = 5 + £ VA , £a = —va + ^6- That is, £4 differs 
from — va by a vector field Vb that does not depend on the branch A. It follows that 
v L — vr = — (£x — We will therefore take 

&r[S] = ®R[8 a ] + ®R[£t] = Q R [6 a ] + t^77 / WCV^gdZ^. (148) 

167TG J So 

as the contribution of Mr to Q[5]. We will also use the conceptual crutch of supposing 
the b chart to be fixed near Sq under the variations being considered, so 5 b = 5 on 
fields there. This condition does not contradict conditions a) and b) on the variations, 
since these conditions refer only to what happens near Sl and Sr, whereas the dif- 
feomorphism generator that must be added to the variations to keep the b chart fixed 
near Sq may be chosen to have support in a closed neighbourhood^ of So excluding 
Sl and Sr. 

In terms of the b coordinates the pullback to Sq of s/—g gEat is 

v^rfS AT = ^V^gzArsndb 3 Adb n (149) 
= -2 pox d [A v + d T] v~ M 1 AcW 2 . (150) 

Since v + and v~ are constant on Mr and Ml respectively dv + and dv~ must be 
proportional to the normals n_ and n + respectively] 65 ! In fact, since d n+ v + = d n _v~ = 
1 it follows that n± = xdv T . Hence 

2 

V^gEat = n r A n +T] e (151) 

^ X 

where e = po dO 1 A d8 2 is the area 2-form on Sq. 
Substituting this form into (I146P yields 

0fi[A] = T7T^/ e-{n_-V n+ e-n + -V re _£}. (152) 

167TG J So X 



64 Thc closure of an open neighbourhood of Sq. 

65 Recall that the tangents to the generators n_ = d v - are normal to AT by prop. IB. 21 
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The derivative along n + may be eliminated in favour of a variation of x- 

5 X = 6g + _ = [8 a g} + ^ + £^ (153) 
= d v+ a a d v -a p 8g a p + 2n*+n v _V {j £ v) . (154) 

But dv-a 13 = r 5P, and (by (11031) ) g ar = —8%, which is constant, so [5 a <?]+- = and 

S X = n + ■ V re _£ + n- ■ V re+ £. (155) 

The integrand of (11521) is thus equal to 

e-Px-2n+.V n _£]. (156) 

X 



Under the variations being considered u remains zero on A/"ro £ is thus tangent to 
Mr and can be expanded as 

£ = £j. + fra-, (157) 

where £~ is the v~ component of £ and £j_ is tangent to the v~ = constant cross 
sections of Mr. 

generates the variations of the y coordinates at fixed 9 on the v~ = constant 
surfaces induced by the variation of the metric. That is, £^ ± y l = 5 b y l . Now by 
definition the curves y % = constant and 8 a = constant, always coincide with the 
generators, a fact which the displacement of the a chart with respect to the b chart 
induced by a variation of the metric must respect. That is, the constant y l curves 
must always remain tangent to the vector field d v - = n_. Thus £ n _y l = d v -y l = 
and £ n _5 b y l = 5 b d v -y l = 0. It follows that 

= [£ n -£u - £a £ n-]y i = o. (158) 

Clearly also £[ n _^ ± ]V ± = 0, so 

[n-,£ A ] = £ n _Z ± = 0. (159) 

The second term in (11561) may therefore be expanded as 

2n + -V n _£ = 2n + -V n _[£ L + f~n_] = 2xd n _C + 2£~n+- V n _n_ + 2n + - V a n_ (160) 

Let us consider each term in turn, beginning with the last. From the definition ( 11 18[> 
of UJ 

2fi + .V ix n_ = d £± x + X^ . (161) 



66 In the definition of the space of solutions S and the space of solutions to the linearized field 
equations it is required that M always remains null. 
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The middle term is a little more subtle. Since n_ is tangent to a geodesic and 
n_ = d v - = r d r 

n + ■ V n _n_ = n + ■ r 2 Q V ' 9r d r = n + ■ r T r rr n_ = x^oKr- ( 162 ) 

T r rr can be expressed in terms of the BFC data by means of one of the field equation 
(we are working in the space of solutions) - see (I133p and Appendix [Dl However this 
will not be necessary since the term we are considering, together with another that 
we will find further on, constitute a total variation which may be subtracted off from 
without altering the presymplectic 2-form. 

It remains to express £~ and d n _^~ in terms of the BFC data and their variations. 
£~ turns out to measure the variation of r (y). Define 

5 y = 5 h -£^. (163) 

(Since 5 b y l = £]_ this is the variation associated with the hybrid moving coordinate 
system (v + , v~, y 1 , y 2 )-) On Mr the variation of r at fixed v~ and y % is then 

5 y r = 6 b r - £^r = 5 a r + £^ n _r = £~d v -r = £~r . (164) 

Hence £~ = 5 y r/r = v~5 y logr . On Sq therefore 

r = 4_r = ^logr = log p v , (165) 

or equivalently 

er = ed n _C = -5*6, (166) 

where 5 y p = 0, which follows from condition b) on 5, has been used. 
Substituting all these results into (11521) we find 

Qr[£<-] = ~Ttr~n I < 6VX + + 5V < 1 + r o r -)' ( 167 ) 

167TG J Sq 

where A = — log |%|. 

Now let us turn to the first term in (|144j) . 0#[5 a ]. In working with this quantity it 
is convenient to use the a coordinates, because in these coordinates the a variation, 5 a , 
of a field reduces to simply the variation, at fixed a a , of the a coordinate components 
of the field. For instance 

= 5V%{a) (168) 

Thus we may write 
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In a coordinates the components of the metric on Mr are (by (1103jl ) 



9a/3 



-1 
-1 



a/3 





-1 



1 





h ij 



;i7o) 



with h l i the inverse of h^. Thus yj—g = ydet[hij] = p y and the pull back to Mr of 
\f—g dE/3 is 

\f—g cE/j = p y 5p dy 1 A dy 2 Adr = S^^Adr. (171) 
( A beneath a form indicates that the form is pulled back to M.) 



The integrand in (116911 thus reduces to 



(172) 



arfc^eAdr = ^{-5T a ar + g^ST ray }eAdr 

-{5d r log + 28Y rur - h ij 5T rij }e Adr ( 1 73) 



2 

--{5<9 r logpy + Sd u g rr + -h lJ Sd r h i:j }e A dr 



(174) 



Now recall that p y (y,r) = r 2 p(y), where p is the area density on Sr. Thus for the 
variations we are considering, which leave p(y) invariant, Sp y = 0. As a consequence 



and 



5d r log p y = 



-5[9 r /iy]/i JJ = <5<9 r logp<, - -d r [p y eij}5 — 



(175) 



(176) 



(Here e* J is the inverse of and the fact that e^Se 1 ^ = — <51ogdet[ejj] = 0, has been 
used.) 

The remaining, middle, term in (11741) is proportional to 5T r rrl since on Mr 



da = 2T r 



;i77) 



Substituting (11751) . (11761) and (11771) into (1 174ft yields an expression for 6fi[5 a ] in 
terms of the BFC data and Y r rr , and their variations. Adding to this the expression 
(11671) for R [£^] we obtain 



I 1 -r^e^e*' - 2r 2 5T r rr dr 



16nG { JSr Jro 2 
J So ) 



(178) 
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Here the identity 

[ feAdr= f e f r 2 fdr, (179) 

J Mr JSr Jro 

valid for any function / on Mr, has been used to turn the integral over Mr into an 
integral over each generator followed by an integration of the result over the set of 
generators. 

The dependence on T r rr can be eliminated. Adding the variation of a functional 
of the data (an exact form on T) to the symplectic potential does not affect its curl, 
the presymplectic 2-form. Thus we are free to subtract from Qr[S] the variation 

= -isRj {I,. f C - 2rHr " dr + L s " e{1 + r ° ry } • (180) 

To obtain the last equality note firstly that Stq is 5 y r evaluated at Sq (i.e. v~ = 1), 
and secondly, that es R ro5 v r = \p(y) dy 1 Ady 2 8 v Tq = |5 y es @ We will therefore take 
as the symplectic potential contributed by Mr 

Adding the contribution of Ml we obtain the entire symplectic potential: 

16ttG I Js R Jr 2 
Jro 2 

+ / e (^ A + ^A + ei^ a -ei^a)j. (182) 

-'So J 

All quantities appearing in the integral over Mr are referred to the cir coordinates and, 
similarly, those appearing in the integral over Ml are referred to the coordinates. 
In the integral over S the term L u; a appears with a minus sign because under the 
interchange of L and R, which interchanges n + and n_, the field uo a = [n + ■ V a n_ — 
n_ ■ V a n + ]/n_ ■ n + is replaced by — o; . 

In (11821) there appear, aside from the BFC data and their variations, only the two 
area forms es^ = pAdy\ A dy\ and the two vector fields £±a- The area forms we know 
how to express in terms of BFC data (see the end of subsection 14. ip . and they will 

67 Strictly speaking the three 2-forms in the equation live in different spaces. But they are equal 
when they are all pulled back to the cotangent space of H 2 using the y chart. 
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be left as they are. The vector £±a is the variation of the map sa that gives the i/a 
coordinates of a generator as a function of its 6 coordinates: C±a = ^Va = $ S A or 5 
transforming to 9 components, 

f)f) a 

= W ^ A . (183) 

sa is the inverse of the BFC datum Ja it may be taken as an element of the BFC 
data, in place of /a- (Notice that the 9 a are the natural coordinates of sa-) 

It turns out to be helpful to make explicit the dependence of the S integral in 
(11821) on the variations of the sa'- 

$y\ = 8X- C±d a \ (184) 

so 

gu*x + gUL X + £ lRUJa _ £ LUa = 26X _ ^ ±RU} _ a _ £l LU+a , (185) 
where oj+ = u + dX and uj- = —uj + dXr3 With the definitions 



86 a , 86 a , 

urn = poU- a -— and uj Lm = p uj +a ——, (186) 
dy l R ayf 

the contribution of Mr to the symplectic potential may be written as 



f [podX-CjRids^de 1 Ad6 2 \. (187) 

J So J 



'So 

and the contribution Q' L [8] of Ml is given by a completely analogous expression. 



4.3 Calculation of the presymplectic 2-form in terms of al- 
most free initial data 

The first term in (11871) is a "bulk" term, an integral over Mr, while the second term, 
an integral over So, is a surface term. There are no contributions from the boundary 
of Mr except that from S . 

The contribution of the hypersurface Mr to the presymplectic 2-form is &r[8i, 82] = 
5iO^[5 2 ] — (1 <-> 2). [For convenience the commutator [^i,^] of the vector fields 81 
and 82 on the space of field configurations T has been set to zero. As discussed in 
footnote EU this may always be done without restricting the action of 8\ and 82 on 

68 These quantities were introduced by Epp in |Epp95| . 
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the fields at the fiducial solution, which is a single point in T\ The bulk term in Qr 
is obtained by varying the bulk term in 0'^ with ro held fixed. It is 



it ej 1 \r 2 b^d r b 2 e l0 dr - (1 ^ 2)1 . (U 



16nG {Jsr 

The surface term at So comes both from the surface term in Q' R and the variation of 
ro in the bulk term of O'^. The latter contribution is 



{Is e5 i r 4 r2 ^ e ^'], - (1 <- 2)} 



167rG usjj 

= i^{W*^-< l ~ 2 >}- (189) 

(The subscript r on the bracket indicates that the bracketed expression is to be 
evaluated at r = r . That is, the integral ranges over Sr with the integrand at each 
point of Sr given by an expression evaluated at r = ro on the generator through that 
point.) 

In (I189p the derivatives by r along the generators may be replaced by derivatives 
by v~ = r/r using r <9 r = d v -, and the variation 5 a may be replaced by 5 y : 

5\r 5 a 2 e 13 - (1 <-> 2) (190) 
= Sfro S y 2 e lj - Sfro 5lr Q d r e i] - (1 2) (191) 
= 5\r 5 v 2 e lj -(1^2). (192) 

The first equality is evidently correct if b\e % i is interpreted as the variation, at fixed 
(v + , v~, y 1 , y 2 ), of the a components e l K The result is stronger than that however. 
Because the 3-metric on TV is degenerate along the generators it takes the form 
hijdy 1 <g> dyi (see I103p . Therefore the metric induced on any two dimensional cross 
section of M will have y components hij. The e u are both the a components of the 
inverse conformal 2-metric, e _1 , on a constant (r, u) surface as well as the compo- 
nents of e~ l on a constant (v + ,v~) surface in the hybrid chart (v + , v~ , y 1 , y 2 ). As a 
consequence the 5 y e* J may also be interpreted as the (v + , v ~, y l , y 2 ) components, or 
the a components, of <5 y e -1 . 

The contribution (11891) to tin can thus be written as 

{J So ^\ d v-e^ - (1 <- 2)} . (193) 



16ttG 

The contribution to Qr of the surface term in Q R [S] is 

- yLg{I s 6lpo62X ~ 5 ^R^ S R dQl A d ° 2 -t 1 ^ 2 )} ' ( 194 ) 
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Summing (TTHHjl . (TT931) . and (TTMD one obtains 

- I S l p 5 2 X - Sxoor^s^ dO 1 A dfl 2 - (1 <-> 2)1. (195) 

J So J 

This is the needed expression for the presymplectic 2-form in terms of the BFC data. 

5 The pre-Poisson bracket on initial data 

With the expression ( 1195ft for the presymplectic 2-form in terms of BFC data in hand 
we are ready to solve for the auxiliary pre-Poisson bracket satisfying ( |53i) : 

8 <p = n[{<p,-}. t 5 ] V5 eL° g . (196) 

Actually, in order to preserve the symmetry between the L and R branches of A/", we 
will solve for an equivalent bracket, {•, -} , on a slightly larger "extended" phase space, 
from which the * brackets between the BFC data may be read off. In subsection 15.21 
the resulting bracket is stated explicitly and discussed. The calculation of the bracket 
is given in subsection (15. 3p . 

Neither the * nor the o bracket satisfies the Jacobi relation. (For this reason they 
are called pre-Poisson brackets and not Poisson brackets.) But recall that (1531) and 
the skew symmetry, linearity and Leibniz rule properties of the * bracket suffice to 
ensure that the * brackets of observables constructed from the initial data are equal 
to the Peierls brackets of these and thus satisfy the Jacobi relation. 

5.1 Preliminaries 

In the following the extended phase space and the o bracket will be defined. In 
addition a useful parametrization of the degrees of freedom encoded in the conformal 
2-metric in terms of a single complex valued field, /x, will be presented. 

5.1.1 A helpful extension of the phase space 

We have introduced the b chart in order to make it easier to treat the two branches 
of M democratically - quantities associated with the 2-surface Sq common to the 
two branches are referred to the b chart rather than the a chart adapted to either of 
the branches. However, to be able to solve ( l53il it is necessary to fix some specific 
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relation between the coordinates 9 a and the y l R and the y° L . One could attempt to find 
a suitably symmetrical relation but this does not seem to be simple. The simplest 
choice is to equate the 9 a with either the y R or the y 3 L , which would of course break 
the symmetry between L and R. 

A way to maintain the symmetry is to maintain the 9 coordinates independent of 
the y coordinates and use a somewhat different bracket. Maintaining the 9 coordinates 
independent introduces a new gauge degree of freedom, namely the freedom to change 
the 9 chart by a diffeomorphism. The value of the presymplectic 2-form on a given 
pair of variations 

does not depend on the 9 chart that has been chosen. Indeed reference to the 
9 chart can be eliminated altogether by choosing 9 = y^ (for instance) and thus 
expressing Q[5\, 62} entirely in terms of the variations of the fields with respect to the 
cll and an charts. Thus a variation of the BFC data which corresponds solely to a 
diffeomorphism of the 9 coordinates, leaving the relationship between the and a R 
coordinates and the fields referred to these invariant, is necessarily mapped to zero 
by Cl. For this reason it is impossible to invert Cl to obtain a bracket on a phase space 
of initial data including fields referred to the 9 chart, unless either the 9 a are taken 
to be fixed functions of the y R and y™ or the phase space is extended by including 
momenta conjugate to the 9 a . 

We shall take essentially the second option, although we will not actually specify 
variables conjugate to the 9 a . What is really necessary is that suitable functions on 
the extended phase space generate the 9 diffeomorphisms. The original phase space 
is then recovered as the submanifold on which these generators vanish. 

Such an extension can be realized in a quite natural manner. In the expression 
(11951) for the R contribution to the presymplectic 2-form and the corresponding for- 
mula for the L contribution the field uj appears in the guise of two distinct quantities, 
CjRi and u)Lm, derived from uj. In the required extended phase space these two are in- 
dependent fields on So- In the original phase space they are of course not independent. 
They are related by the equation 

"RiQ^ + VLmdfit = P^-* + = 2p df- X) (197) 

which serves as a constraint that defines the original phase space of initial data as a 
submanifold of the extended phase space. As we shall see shortly this constraint also 
generates diffeomorphisms of the 9 chart. 

The symplectic potential on the extended phase space, denoted 0, is taken to be 
equal to the expression 1 1 187ft for 0' but with uj Ri and LOLj independent fields. The 
presymplectic 2-form Cl is the curl of 0. Explicitly Cl is the sum of a contribution 
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from Af R , Q R = A R + B R + C R with 

A R [S U 8 2 ] = -j^l 8 lPo 8 2 X-8 1 u Ri 8 2 s R d9 1 A d9 2 - (1^2) (198) 
107rG J So 

Br[5 1 ,5 2 ] = -i- / 5\ed v -e l3 5 y 2 e^ - (1 <- 2) (199) 

647TG J5o 

Cr[8iM = I e f 1 r 2 5 1 e^d r 5 2 e i] dr-{\^ 2), (200) 

OZTlU JS R Jr 



and an entirely analogous contribution Cl^ from Ml (obtained by substituting L for 
R, and t> + for v~, in Or). 

The auxiliary bracket on the extended phase space, {•, -} , is defined by the re- 
quirement that 

6<p = n[{tp,-} ,8\ (201) 

for all extended phase space initial data tp on the interior of M and all variations 5 of 
the data compatible with the metric consistency constraint fllQ9p . and the condition 
that Seij = at the two ends S R and Sl of M . These variations 5 will be termed 
"admissible" . 

Note that we have weakened the condition on the variation Se^ of the a chart 
components of e. We are not requiring it to vanish in a neighbourhood of dj\f, only 
on DM itself. This does not lead to any actual change of the implications of f)20ip for 
the specific Cl we are using. Restricting the admissible S to ones for which Se^ vanishes 
in some neighbourhood of dj\f does not enlarge the set of brackets solving (12011) . as 
the reader may easily verify when this condition is solved for {-,-} in subsection 
15.31 But a simpler definition of the admissible 5 seems preferable, and may facilitate 
understanding of the results, especially in connection with the Jacobi relations. In 
any case it is clear that enlarging the class of admissible 5 does not weaken (12011 ). 
and so our key result, that the brackets satisfying this relation reproduce the Peierls 
brackets of observables (when the constraint (11971) is imposed) continues to hold. 

Let us return to the constraint ( 1197jl . Imposing this constraint on the interior of 
So is equivalent to requiring that 

k[v] = -1— I v a [2pd a X — uj Ri d a s R — LOLmdaS^dO 1 A d6 2 (202) 

1D7TG J So 

= -&[£v], (203) 

vanishes for all smooth vector fields v tangent to So and vanishing on dSo- k[v] acts as 
a diffeomorphism generator, it produces a variation S y 6 a = —v a in the 9 coordinates 
at fixed y R and y^, and no change in the relation between the a R and ai charts or 
in the fields as functions of either of those charts. Equivalently it acts as £ v on the 
fields as functions of the 9 a . 
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To see this note that 

Cl[£ v , 5} = £ V Q[5] - 8G[£ V ], (204) 

where 5v has been set to zero, so that [5, £ v ] = £$ v = 0. £ v acts only on the Sq 
surface term in Q[5] since the bulk term, in which no fields referred to the 9 a intervene, 
is invariant under diffeomorphisms of the 9 chart. But since the integrand of the So 
term is a scalar density under diffeomorphisms of the 9 and v vanishes on 8Sq the 
Sq contribution to Q[S] is also invariant under the diffeomorphisms generated by v. 
Thus <£„9[5] = and 

6k[v] = &[£„, 6\. (205) 

Note that 5 need not be an admissible variation. Admissibility is not assumed in (I204jl 
nor in obtaining the explicit expression (11981 - [2001) for Cl from O, via the calculation 
of subsection 14.31 (Indeed, in 14.31 the variations are only assumed to satisfy the 
conditions a) and b) of subsection 14.11 ) We may therefore take 5 to be the action 
{</?, -} of an initial datum ip, in which case 

{k[v],<p} = -Sk[v] = Q[{p, -} , £ v }. (206) 

Since £ v is an admissible variation (in fact, it leaves invariant) equation (12011) 
implies that 

{n[v],ip} = £ v ip. (207) 

Here £ v acts only on the 9 dependence of the fields. That is, it is equivalent to a 
generator of diffeomorphisms of the 9 coordinates. 

Because it preserves the symmetry between the branches of M the bracket {•, -} 
rather than {•,•}* will be given explicitly here. Note however that the bracket {•,•}* 
on the original phase space corresponding to the choice 9 1 = y\, 9 2 = y\ (i.e. sl = id) 
can be recovered easily from {•, -} as a Dirac bracket: 

{-, •}* = {■, -}o - / jtt [{; s l L } {K a , -} - {-, K a } {s l L , -} ] d9 l A d9\ (208) 

where n a = n[d a ] = j^[2pd a \ - u Ri d a s l R - uj Lm d a sf}, and with the fixing of 
9 adopted is a Kronecker delta. Of course the * bracket may also be obtained by 
solving the original equation ( 1531) with this fixing of the 9 chart. 

The explicit form of the o brackets, given in subsection 15.21 shows that the only 
datum that has non-zero o bracket with Sl is ujl (see (124 ip and the subsequent 
discussion). {•, •}* thus differs from {•, -} Q only when one of the arguments is u>£. This 
datum is not independent on the constraint manifold so the * brackets of a complete 
set of phase space coordinates on the gauge fixed constraint manifold, namely the data 
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we have been considering excluding sl and ujl, are identical to the corresponding o 
brackets. The * brackets involving Cjl can be obtained from brackets between these 
other data by substituting the expression for &l in terms of these provided by the 
constraint (I197p and the gauge condition sl = id. The * brackets of of course 
vanish. 



5.1.2 A useful parametrization of the conformal 2-metric degrees of free- 
dom 

We turn now to the definition of a useful parametrization of the conformal 2-metric 
e. The unimodular symmetric matrix field e a t, is free initial data in the sense that it is 
unrestricted by the field equations, but of course the components of e are algebraically 
restricted by the requirements of symmetry and unimodularity. It is thus convenient 
to parametrize the two degrees of freedom of e by an unconstrained complex number. 
We shall do this as follows: The degenerate line element on Af can always be expressed 
in terms the complex coordinate z = 9 1 + i9 2 as 

ds 2 = h ab d9 a d9 b = p(l - \fi\ 2 y l \dz + fidz][dz + pdz]. (209) 

with a complex number valued field of modulus less than 1. \x provides the desired 
parametrization of e. In component form (with respect to the 9 coordinate basis) 



e a b 



2 



|l + /i| 2 2Inu\i 
2Im.fi |1 — jj\ 2 



(210) 



Under a transformation of coordinates 9 — > 9' in which the 9' coordinates, like the 9 
coordinates, are constant on the generators of A/", the components e a b are transformed 
to 

d9 a d9 b , , dff . / * 

^ = W^ det m Kb (211) 

(a fact that was used in the definition of the constraint (I109P ). The corresponding 
transformation of fi is 

, . dz dz , , r dz dz . . 

"^"=fe + "a?l/[a? + "a?l' (212) 

The \x field referred to the ur chart will be denoted by fiR and that referred to 
the dl chart by fj,£,. 



5.2 The brackets 

In this subsection the o brackets will be given between the data u)m, Cjim, Sr, s£, 
A, po, and \x and p, (or equivalently e a b), where the conformal 2-metric data \x and 
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Jl are, like the other data, referred to the b chart. A subset of these form a free and 
complete data set. Moreover, the o brackets among the data in this subset are equal 
to their * brackets. 

As pointed out earlier, the constraint (11091) that the or and components of e 
must satisfy at S is equivalent to the requirement that e ab at So obtained by trans- 
forming eRij from the clr chart to the b chart is the same as obtained by transforming 
ZLmn from the ol chart, i.e. it is equivalent to requiring that e a b can be consistently 
defined. The constraint (11091) thus disappears if the b components, e a b, are used to 
specify the conformal 2- metric throughout M , in place of eRij on Mr and e^ mn on 
Ml- 

The other constraint, (jl97j) . is still present, and so is, necessarily, a gauge fixing 
condition defining the 9 coordinates. However, as shown in subsection 15.1.11 if the 
simple gauge fixing sl = id is adopted then (1197ft can always be solved for Cjl, and 
the remaining data e a b, Cjri, s r , A, and p form a free and complete set. Using of the 
explicit expressions for the o brackets given here it is easy to check that the * 



coincides with the o bracket between these data. To maintain the symmetry between 
L and R in the formalism, and in order to not commit ourselves to any specific gauge 
for the 9 coordinates, the o brackets between all the data, including ljl and sl will 
be given. 

For 1 and 2 points on Sq (or more precisely values of the b coordinates correspond- 
ing to points on S ) 

{A(l),p (2)} o = -8nG5 2 (9(2)-9(l)), (213) 
{^(1),4(2)} = -167rG^ 2 (fl(2)-fl(l)), (214) 
{u R [f],X(l)} = 7rG[(d v -e ab -d v+ e ab )£ f e ab ] 1 (215) 

{w*[/i],£fl[/ 2 ]}o = -2ttG / £ h e ab d v -e ab £ h e-(l++2) (216) 

J So 

{ou R [f],u L [g}} = 1ixG\ [£ f e ab d v+ e ab £ g e- £ g e ab d v -e ab £ f e], (217) 

J So 

where Cj R [f] = J So LJ Ri f l d9 1 A d9 2 with / a test vector field tangent to So, and Cji[g\ 
and g are defined similarly. Note that the jjr components f l of /, and similarly the 
Ul components g m of g, are taken to be independent of the dynamical fields and thus 
o commute with everything. This is important to remember in order to interpret the 
given expressions for the brackets correctly, as can be clearly seen when the smearing 
fields are eliminated. 

Expressions analogous to ( 12141) . ( 12151) . and (12 16f) give the brackets {d>£ m (l), s£(2)} , 
{\(l),u L [g]} , and {dj L [gi},u L [g 2 }}o respectively. 

For 1 and 2 points on the same branch of M 

{//(l),/2(2)} = 2TxG5\9(2)-9(l))sgn(l,2) 
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where for 1 and 2 lying on the same generator (i.e. 6(1) = 6(2)) 

(1 2) — { ^ ^ * S ^ ur ^ ner fr° m &o along the generator than 1 (219) 
S ^ 1 — 1 if 2 is closer to Sq along the generator than 1 

and the integral in the exponential is evaluated along the segment of the generator 
from 1 to 2. Note that this defines the bracket between /i and /i smeared with smooth 
test functions on J\f. 

When 1 and 2 both lie on So {/x(l), /i(2)} = 0. This is not unambiguously implicit 
in f)218p . It is an additional specification of the bracket compatible with the defining 
equation (12011) . 

For 1 and 2 any two points in intM 

= {Ml),M2)}o = {/i(l),/i(2)} o , (220) 

Notice that when 1 and 2 both lie on So (12201) . together with the vanishing of 
{/x(l), /i(2)} in this case, implies that {e cd (l),e a b(2)} = 0. 

These brackets are not uniquely determined by (I20ip . In particular 5 2 (6(2) — 
6(1)) sgn(l, 2) in (J2IED may be replaced by 5 2 (6(2) - 6(1)) sgn(l, 2) + iT(6(l), 6(2)), 
where T is a real distribution on 6, and possibly a function of the data, resembling 
5 2 (6(2) — 6(1)) in that it has support only at 6(2) = 6(1) and that it transforms like 
5 2 (6(2) — 6(1)) under diffeomorphisms of the 6 chart. 

This indeterminacy arises for the following reason: Denote the a chart \i and p, 
fields on the truncating surface Sa C dM of a given branch by fi and /2. The bracket is 
only defined on data on the interior of A/", but we may define {fi(y), -} and {ft(y), -} 
to be the limits of {/i(p), -} and {fi(p), -} , respectively, as p tends to Sa along the 
generator specified by y. For any admissible variation 5 (which we always take to be 
smooth) 

Cl[{fi(y), -} , 5} = lim Cl[{n(p), -} , 5} = lim 5(i(p) = 0, (221) 

P^Sa p^Sa 
and similarly for {fl(y),-} . (12011) thus permits 

the addition of multiples of {fi(y), -} and {fi(y), -} to any solution {(pj, -} . This 
freedom in the brackets is greatly reduced by the requirements of causality and anti- 
symmetry, but not entirely eliminated. 

Another possible indeterminacy is associated with the bracket {e cd (l), e a b(2)} for 
1,2 G So- The vanishing of this bracket is shown to be consistent with (12011) . but 
it is not shown to be necessary. On the other hand, no candidate for an alternative 
solution for {e cd (l),e ab (2)} was found. In other words, this bracket has not been 
demonstrated to be uniquely determined, but it may be so. 
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The freedom in the bracket {/i(l), /i(2)} propagates to some other brackets. 
However, if it is assumed that {/i(l), /i(2)} is given by the expression (12X81) and 
{e cd (l), e ab {2)} = for 1, 2 e So all other brackets are uniquely determined by (12011) . 
Here I will make this choice, as it is the simplest and provides the most natural 
seeming solution to (1201 p . 

For 1 e S and 2 e Mr - S 

{A(l),e afe (2)} = ^[^ + e cd ] 1 {e cd (l),e a6 (2)} - 4vrG[— d v -e ab ] 2 5 2 (9(2) - 9(1)), 
4 Po 

(222) 

{u R [f],e ab (2)} = -\ I e£ f e cd {e cd ,e ab (2)} 

Z J So 

-lQnG[£ f e ab - \d v -e ab ^\ 2 , (223) 
£ Po 



{u L \g],e ab (2)} = -f [e£ g e cd - £ g e d v +e cd }{e cd ,e ab (2)} . 

£ J So 



'So 

(224) 

Of course analogous expressions with the roles of Mr and Ml reversed give the cor- 
responding brackets when 2 G Ml — So- 
When 2 lies on So 

{A(l),e Qfc (2)} = -27iG-[d v -e ab + d v+ e ab ]5 2 (6(2)-6(l)), (225) 

Po 

{u R [f},e ab (2)} = -^G[£ f e ab -\d v -e ab ^] 2 . (226) 

^ Po 



1 £ 

{u L [g],e ab (2)} = -87rG[£ g e ab --d v+ e ab ^} 2 . (227) 

^ Po 



These brackets for 2 e So are the averages of the expressions for the same brackets 
valid in the cases 2 e Mr - S and 2 e M L - S . (Recall that {e cd (l), e a6 (2)} = 
when both 1 and 2 lie on So.) 

All other brackets vanish. In particular 

= {A(1),A(2)} . (228) 

and sr, sl and po commute (that is, they have vanishing o brackets) with everything 
except (2>r, ul, and A respectively. 

The requirement that k generates diffeomorphisms of the 9 chart - specifically 
that <p} = —£v i P for all data (p on intM - provides a non-trivial check of the 

brackets. The reader may easily verify that the expressions for the brackets given 
here satisfy this requirement. 
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Notice that by (1218j) and (12201) the variation {//(l), -} does not preserve the reality 
of the metric. This corresponds to the fact that \x is complex. The action of a real 
functional of fi and p does preserve the reality of the metric, as can be seen easily 
from the relations 



{/j(l),/2(2)}o = 02(l),M2)}o, (229) 
{/i(l),M2)}o = {/2(l),/2(2)} 0) (230) 

which follow from (EIBjl and (12201) . (12291) and 

( 12291) can be interpreted as asserting that the bracket in itself, {•, -} , is real, so 
that complex conjugation acts only on it's arguments. 

Since the coordinate systems yn, y^, and 9 labelling the generators of M are 
essentially arbitrary it is important to verify that the brackets are invariant under 
diffeomorphisms of these coordinates. All the data except \x and p, transform as tensor 
or tensor density fields under such transformations. For all the brackets except those 
involving \i and p it is therefore easy to check that the brackets and the expressions 
proposed for them transform in the same way under the various possible coordinate 
transformations. Moreover, the transformation law ( 12 121) of /i shows immediately 
that {/i(l), /i(2)} = is a coordinate invariant statement. The only non-trivial case 
is therefore the bracket ( 12 18f) . 

To verify the invariance under coordinate transformations of the equation (1 2 1 S |) 
note first that 5 2 (8(2) — 6(1)) sgn(l, 2)-^=-^= is invariant. From (I212j) one shows that 

a change of coordinates transforms the ratio of a variation A/x and 1 — |/i| 2 by a phase 

factor: 

A// a + fif3 A/i 



1 — \fi'\ 2 a + /jl(3 1 — |/i 

dz 7 I-' ~ W- 



where a = -fr and /3 — Thus 



WW (2)} 



[^§]i[^4] 2 [ r ^]i{Ml),M2)}o[-^] 2 . (232) 

To complete the demonstration of the invariance of ( 12 181) we only have to show that 
the phase exp(/ 1 2 1 _^| 2 [fid/i — /idp]) transforms in the same way. 
( 123T1) and (12T2D show that 

fi'd/jf (3 + p,a dfi 



|/i'| 2 a + fi/31 — |/i| 2 
j3dfi \id\x 
a + fi/3 1 — \/i 



(233) 



+ -r^frr, (234) 
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It follows that indeed 




(236) 



The equation (12181) is therefore invariant under transformations of the 9 chart. 

( 12181) has a further invariance. It is invariant under reparametrization of the 
generators. In fact (12181) does not even refer to a particular parametrization of the 
generators. This makes it easy to take the flat spacetime limit of this bracket, and 
it suggests the conjecture that {fi, p,} is given by the expression (12 181) even when 
the area distance is not a good parameter along the generators. For instance, in flat 
spacetime with M swept out by parallel generators (neither diverging nor converging), 
or in a generic spacetime along a generator with neighbours which initially diverge 
and then re converge. Checking this conjecture is left to future work. 

Invariance of the expressions for the brackets under transformations of the 9 coor- 
dinates is necessary because these coordinates are arbitrary, but this is not so for the 
parametrization of the generators, because v = r/r is not arbitrary. It is completely 
determined by the geometry of N '. Indeed the expressions for several of the brackets, 
involving d v -e and d v +e at So, are not invariant under such reparametrizations of 
the generators. Of course one may still re express these brackets in terms of another 
parametrization, A, of the generators by replacing d v e a b at So by 



One sees that the brackets cannot all be straightforwardly extended to hypersurfaces 
M with vanishing expansion of the generators at So- The area distance is not a good 
parameter at So in this case, and when it is replaced by a good parameter, such as 
an affine parameter, several of the brackets become undefined. The comparison of 
the present theory with canonical formulations of general relativity linearized about 
Minkowski space |AG73] in terms of data on null hyperplanes is thus a delicate matter 
which is also left to future investigations. 

5.3 Calculation of the brackets 

How were these brackets obtained from (12011) ? In the following (120 ip will be solved 
step by step, yielding the brackets given above and at the same time demonstrating 
that {■, -} defined by these brackets really satisfies ( 12011) . 




'\e a b = 




d\e ab 



(237) 
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We begin by considering ( 12011) for variations S such that = Se = 5sl = Ssr = 
5po- Then 



n[{<p(l),-} ,6\=A[{<p(l),-} ,S\ = 

2{cp(l), po}o6\ + Ml), s R } Su Ri + Ml), s?} 6u> Lm dO 1 A d^38) 



167rG J So 

Substituting for ip each of the fields coordinatizing the extended phase space one finds 
that the only non-zero bracket involving po is (12131) . 

{X(l),p (2)} o = -8nG5 2 (9(2) - 9(1)), (239) 

the only non-zero bracket with sr is (12141) . 

{(b Ri {l), 4(2) } = -1^G515 2 (6(2) - 6(1)), (240) 
and the only non-zero bracket with sl is 

{iu Lm (l), s n L (2)} = -lQnG5l5 2 (6(2) - 9(1)). (241) 

Moreover, it is easy to see that, given these non-zero brackets and the fact that all 
other brackets involving p , sr, and sl vanish, (12011) is satisfied for <p = p , s R or 
and all S. 

Now let us consider (I20ip with <p = pr evaluated at fixed values of the cir coor- 
dinates corresponding to a point 1 G Mr — So, and admissible variations 5 such that 
= 5sl = 5s r = 5po but 5e not necessarily zero (except of course on dj\f). This will 
give us the o brackets of p with p and pr 9 \ 

To solve this condition we first express fl[{p(l), -} , 5} in terms of p, p, and their 
variations. Since we are only considering variations 5 that satisfy = 5sr = 5po and 
since = {p, s l R } = {p, p } o 

{Ir[{p(1),-} ,5}=Cr[{p(1),-} ,5} = ^ J SR ej\%e^d r 6 2 e lJ dr-(1^2), 

(242) 

with 5\ = {p(l), - } and 62 = 6. To this we must add an entirely analogous expression 
contributed by Ml- 



69 To lighten the notation the subscripts R and L indicating the branch to which variables are 
referred will be dropped when the context makes this information redundant. In particular \x will 
be used to denote variously fj,R, ^ referred to the 8 chart, or several of these indistinctly, as the 
context indicates. For example, in the following paragraphs (until equation (|263p ). concerned with 
the brackets of )J,r(1) for 1 £ Mr — Sq, /x(1) will denote /j,r{1) and plain /1 denotes (xr, hl, or either 
depending on the context. 
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[Note that while 62 is smooth, admissible, and real in the sense that 62/2 is the 
complex conjugate of 621*1, 61 = {/i(l), -} need not have any of these properties (and 
indeed, it turns out to have none of them). In particular we shall see that 6\jx is 
not the complex conjugate of 61 fi, so this variation takes us outside the space of real 
metric fields. All the same there is no ambiguity in the definition of 6iey. Even for 
real variations, such as 6, 6fi and 5ft are linearly independent, that is, no non-zero 
pair of constants a, b exist such that at a given point a6fi + b5ft = for all possible 
real 6. eij thus has uniquely defined partial derivatives by fi and ft, which in turn 
define 6\e^ unambiguously, even when 5\ft 7^ 6\n] 

In terms of fi 

6 ie ij d r 6 2 e i3 - (1 2) 

= _4 7i r~m? {8iftd r 5 2 ft + 6iJid r 6 2 iJi} 

(1 - H 2 ) 2 

_8 7^ i — to \p9 r ft - /i9 r /i]5i/x5 2 /i - (1 2). (243) 

(1 - \fi\ 2 y 

[This can be obtained as follows: By a direct calculation from (12101) one sees that 

Ae tJ Ae ij = -8— j-— -A/x Aft, (244) 
(1 - |/i| 2 ) 2 

for any variation A. (This is valid even when A/2 is not the complex conjugate of 
A/i.) Now, substituting Ai + A2 for A in (12441) and taking the part linear in Ai one 
obtains 

A ieij A 2 e ij = -4- —-{A 1 fxA 2 fx + Ai/2A 2 /x}. (245) 



In particular, with Ai = d r and A2 = 61 



r ey<Jie tJ = -4- — {d r ft5 x ft + <9 r /i<5i/i}. (246) 
(1 - |^| 2 ) 2 



Finally, by acting with 62 on (I246p and anti-symmetrizing with respect to interchange 
of 61 and 62 one obtains (12431) .] 
If we let 

f r 1 

a ( r ^y)= 1 7-r-lfidrfi - (idrfi] dr, (247) 

Jr 1 - |/i| 2 

where the integral is taken at constant y\ that is, along a generator of M , then 
6 1 e lj d r 6 2 eij — (1 2) 

~ -^^i^ 8 0- (1 " 2) ' (248) 
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(Note that a is pure imaginary, so e a is a phase) 
This invites us to define, for any variation A, 



□ = ^^A/i □ = f^Ap. (249) 

□ is a formal complex conjugate of □, obtained from □ via the substitutions p <-> p 
and Ap <-» A/2. When A is a real variation, so that Ap = Ap, □ really is the complex 
conjugate of □, but for variations, such as Si, which are not real this is not the case. 
With these definitions 

C[{p(l),-} ,5] = --^—[ [ dy R Ady 2 R f U Rl d r n R2 + U Rl d T U R2 dr 

87TG JSr Jr R0 

+ / dyiAdylf a L1 d r a L2 + n L1 d r a L2 dr 

JS L Jr L0 

- (1 <-> 2)]. (250) 

(□l and O r signify □ formed according to (I249p from fields on Ml, referred to the 
(ll chart, and fields on Mr, referred to the a R chart, respectively.) 

The contribution of either branch to f)250p may be integrated by parts in two 
different ways: 

f □ 1 (9 r .D 2 + Uj^Da - n 2 d r Ui - U 2 d r U x dr = 2/ Uid r n 2 + U x d r U 2 dr 

Jro Jro 

+[n 1 n 2 + n 1 n 2 ] ro , (251) 



or 



-2 f n 2 d r a 1 + u 2 d r u x dr 
-[□^a + niDa]^. (252) 



In both these expressions the boundary term at r = 1 (i.e. at dM) vanishes because 
Seij and thus n 2 and 2 vanish on dM. 

C[{/i(l), -} ,S] can be put in a form that facilitates the solution of (12011) by in- 
tegrating the contribution of Mr by parts one way, according to (12521) . and the 
contribution of Ml the other way, according to (12511) . For variations 5 satisfying 
= Spo = 5s R = Ssl the surface terms of the two contributions then cancel: By 
(12491) and ([2451) 

□ 1 U 2 +U 1 D 2 = P ? ^J StfiS^ + S^fi} (253) 



I/' 



—pySie^S^K (254) 
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Since both 5\ = {/^(l), -} and 5 2 = 5 leave po, Un and Ul invariant Si R = 5^ = 5i L 
and 5 2 R = ^2 = $2 L ■ The fact that 5\ and 62 both respect the constraint (11091) thus 
implies that 

/ PyR6l R e Rij 6Z*e R ij dy l R Ady 2 R = f PyL ^e Lmn ^e L mn dy\ A dy 2 L . (255) 

J So J So 

That is, the surface terms cancel and 

C[{/2(i),-} ,5] = -^—[ +( dy R Ady 2 R f 1 n R2 d r n R1 + n R2 d r n R1 dr 

47TG JS r Jr 

-f dy\ Ady 2 L [ n L1 d r n L2 + U L1 d r n L2 dr]. (256) 

JS L Jr 

(Further on we shall need to evaluate C[{/i(l), -} , 5} for arbitrary admissible 5. It is 
therefore worth noting that (l23U|) . (gHTJ , (ES2D, and (12511 hold for all admissible 5.) 
We are now ready to solve 

fy(l) = 6[{//(l),-}o,tf] (257) 

for /x(2)} and {/x(l), /x(2)} . Since 1 G jVr — our causality requirement 

implies that both these brackets vanish when 2 £ A/z, — So, because 2 is then outside 
the causal domain of influence, J [I], of 1. (See appendix iBl prop. IB. 711 . Thus in the 
second term of (12561) . the contribution of Ml, n Li and D L1 can have support only on 
S itself. Could they be singular distributions supported there? 5/Il and SJjll (and 
therefore d r O L2 and d r O L2 ) can be varied freely in a neighbourhood of S in A/l@ 
Moreover, both are independent of 5fi(l). It follows that n L1 and D L1 must vanish 
as distributions on Ml — dM, and that the contribution of Ml to ( 1256j) vanishes. 

On Mr R2 and D R2 are freely variable except that they must vanish on dM. 
Thus ( I257P requires that 

d r a R1 = (258) 

and 

-V^la r n i?1 = 4nG5 2 (y - y(l))S(r - r(l))a (259) 
1 — luf 

as distributions on Mr — dM. Here a is +1 if r increases away from So (i.e. if Tq < 1) 
and — 1 if r decreases away from So (ro > 1). This factor must be inserted because 
the integral over an interval [a, b] of a delta distribution with support between a and 
6 is +1 if b > a and — 1 if b < a. 



70 5/il and Sfl^ are not independent since one is the complex conjugate of the other, but they 
are linearly independent, that is no non-zero pair of constants a, b exist such that at a given point 
adjiL + bSpi = for all possible variations. This is sufficient for our purposes. 
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(12581) and (Egg}) are satisfied by 



{//(l),//(2)} = (260) 
Oi(l),0(2)} o = 2vrGc%(2)-y(l))^(r(2)-r(l)) C r 



for 1 G A/# — So and 2 G jVr, where 

S9n( r(2)-r ( i))={+; i|$j<$) (26 2) 

Notice that when 1 and 2 lie on the same generator sgn(r(2) — r(l))a = sgn(l, 2) 
and that fzp2 \fid r fj, — fj,d r p]dr = J 2 x _f , 2 [/trf/t — /id/7] , so (12611) differs from (I218p 
only in that the chart is used to label the generators instead of the # chart. 

These brackets are not the only possible solutions to ( 12571) . ( I258p and ( 12591) allow 
the addition of functions independent of r to O r1 and to O ri . The requirement 
that the brackets be antisymmetric reduces this freedom. In the case of {/i(l), /a(2)} 
it eliminates the ambiguity all together - (12601) is the only possibility. In the case 
of {/i(l), /i(2)} one is still free to replace 5 2 (y(2) — y{l)) sgn{r{2) — r{l))a in 12611 
by 5 2 {y{2) — y[l)) sgn[r{2) — r(l))<r + iT(y(l),y(2)), where T is real. Causality 
requires that T is a singular distribution with support only on y(l) = y(2), and 
covariance under transformations of the y coordinates requires that it transforms 
as 5 2 {y{2) — y{l)) does. Some possibilities are T = constant x 5 2 {y{2) — y{l)) or 
T = \{y)x5 2 {y{2)-y{\)). 

The freedom in the bracket {//(l), //(2)} propagates to other brackets. Here the 
simplest bracket, with T = 0, will be adopted since it seems the simplest and most 
natural seeming choice. 

What about the brackets between \x R on Mr — So and /ix, and ^ on Ax? Causality 
and the constraint (11091) imply that 

{ yUR (l),/i i (2)} = (263) 

The case 2 G So requires explanation. In this case the constraint (I109p implies that 
11l(2) may be expressed as a function of Hr{2), and the coordinate transformation 
from the complex coordinate z R = y R + iy 2 R to Zl = y\ + iy\'- 

i \ \® Zr , dz R dz R dz R 
Afc(Ain) = [ Wl + ^g-]/[^- L + A**^], (265) 
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(see (1212jl ). Similarly lil may be expressed in terms of //r, and the coordinate trans- 
formation, via the formal complex conjugate of (12651) . Since /zr(1) commutes with 
Vr, Si> Sij and p 

{li r {1),li l (li r ;2)} o = (266) 

and 

{li r (1), fi L (ii R - 2)} D = [|^] 2 {^(1), fl R {2))} . (267) 

(Note that we are using the fact that by (12651) ill depends only on ll r and the 
coordinate transformation, and not on p R .) The fact that {/i#(l), l1l{2)} does not 
vanish on all of Ml does not contradict the requirement that n^i vanishes as a 
distribution, for if L;1 as defined by (jZMj) , (I265j) . and (126 If) is integrated against a 
smooth function on Ml the result vanishes. 

Note that although the variation of e produced by {/i#(l), -} is discontinuous at 
S (I264p assures that {/i#(l),-} respects the constraint (11091) . for the constraint is 
a requirement of consistency rather than one of continuity. The constraint requires 
that eRij and eimn define the same e a b on So, and {(i R (l), -} respects this condition. 

Instead of parametrizing e by ji R on M R and by /il on Ml, one may parametrize e 
on all of M by \x referred to the b chart. When solving (12011) we wished to use variables 
with vanishing variation on dM under admissible variations of the solution. However, 
now that we have solved (12011) for the brackets between the /i and Jx fields nothing 
stops us from transforming these brackets into the b chart. Indeed this transformation 
simplifies the brackets somewhat. 

To carry out the transformation first note that since the brackets of ji (and Jx) with 
s l R , s™, and po vanish the transformation from the a to the b coordinates, although 
really field dependent, can be treated as field independent when transforming the 
brackets of fi and /L Second, recall that the only difference between (12611) and ( 12181) 
for points on M R is that in the former the y R chart is used to label the generators 
while in the latter these are labelled by the 9 coordinates. The same is true for the 
two versions, (I260p and (12201) . of the bracket between /i(l) and /i(2) on M R . But we 
have shown that (12201) and (I218j) are invariant under arbitrary diffeomorphisms of the 
9 coordinates, and thus in particular under the transformation from the 9 chart to 
the dr chart and vice, versa. Transforming (12601) and (12611) - and their analogs on 
M L - to the 9 chart thus yields exactly (12201) and (12181) . 

{/i(l),/i(2)} = {/i(l),/2(2)} o = 0, (268) 
{/i(l),/2(2)} = 2nG5 2 (9(2)-9{l))sgn(l,2) 

[llJd!] 1 [illH!] 2 Jl iq^ [ ^-^ ] , (269) 
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for points 1 and 2 belonging to the same branch of M with at most one of these lying 
on So- (12631) and (12641) show that the brackets between these fields at points 1 and 
2 belonging to distinct branches, with neither lying on So, vanishes. If (12681) and 
fl2HD hold then the formulae (12611]) and (l2"Mj) for the brackets of ^ R {l) with fi L (2) 
and /ii(2) for 2 £ Sq follow as automatic consequences. 

( 12681) and ( 1269[) also suggest values for the brackets when 1 and 2 both belong to 
Sq. They suggest that 



for 1,2 £ So. In applying (I268H and (12691) to the case 1,2 £ So we are of course 
extrapolating. These equations have not yet been shown to be consistent with (12011) 
in this case. Moreover, (1269}) . which is unambiguous when /i(2) is smeared with a 
smooth function on A/", is ambiguous when 1 and 2 are restricted to So because the 
discontinuous function sgn(l, 2) has not been defined except when 1 and 2 lie on the 
same generator. If this function is discontinuous on 5*0 (which sgn(r(2) — r(l))a is) 
then the product 5 2 (8(1) — 6(2))sgn(l,2) is ambiguous. We therefore need to check 
the compatibility of the proposed brackets (I270j) with (12011) . 

As a first step let us verify that ^[{^(1), -} ,<5] = 6/j,r(1) for 1 £ 5*0 is satisfied 
for admissible 5 such that Ss^ = 5s™ = 5p = provided all the expressions (12601) . 
(12611) . (12631) . and (12641) for the brackets among the /i and ft fields obtained so far, as 
well as the proposed brackets (I270p . hold. 

To do this we integrate both terms in ( 12501) by parts according to ( 12511) . Then 



(Recall that Di corresponds to 5i = -} and D 2 to <5 2 = S.) 

By (1270p the So integrals vanish. ( 12601) and its analog for Ml imply that 



o = {Mi),M2)}o = Mi),M2)}o = 



o 



(270) 




(271) 







(272) 



and (1261]) implies that0 




O 



(273) 



Recall that □ is need not be the complex conjugate of □. 
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2ttG 



1 — 



S 2 (y R -y R (l)) 



(274) 



IB 



on M R — S , and (I270p implies that it is zero on So itself. Therefore 
- f dy R A dy\ f U R1 d r n R2 + U R1 d r U R2 dr 



1 



i/' 



2 



-2ttG 

2ttG 
2ttG 



''o 



d r a R2 dr] y -_ 



-I 1R 



Similarly, at any point of Ml 



□ 



L 1 



p y e 



{lI R (l),II L }o 



(275) 

(276) 
(277) 

(278) 



J L 



The bracket {/i R (i), //£,(2)} with 1 G So an d 2 G Ml — So may be evaluated using 
the complex conjugate of the analog for Ml of (12641) . yielding 



□ 



L 1 



p y e 



dfiR 
dp L 



{fi L (l),fi L } . 



(279) 



It follows, by the analog of the calculation carried out in (12741 - 12771) . that 

-/ dyiAdyl / 1 n Li a r U i2 + U il 9 r D L2 rfr = 2ttG fe 1( ^ L (l) (280) 

J St. Jrn Ulll. 



dp L 
2ttG 5 aL jj R (l] 



(281) 



5 aL fi R (l) is 5 aL fi evaluated in the a R chart at the values a R (l) of the a R coordinates. 
By (12121) fi transforms autonomously under change of coordinates, i.e. \i in the new 
chart depends only on [i in the old chart and on the transformation between the charts 
- it does not depend on other fields, such as fi, in the old coordinates. S aL //r(1) is 
thus defined by (11251) . This definition implies that 



5 aL fi R = dfi R /dp L 5 aL PL, 



(282) 



where in the partial derivative fi R is /i referred to to the fixed chart a R o which agrees 
with a R at the fiducial solution, and similarly \xl is referred to the fixed chart a^o 
defined by to at the fiducial solution. The partial derivative is of course evaluated 
at fixed qlo (and thus fixed a R o). 
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The sum of all the terms in (1271j) is thus 

U a ^ R {l) + U a ^ R {\). (283) 

Note that this evaluation of C[{/xr(1), -} , 8} is valid for all admissible 8 - the hy- 
potheses that 8s l R = 8s m = 5po = have not been used. At present however we are 
interested in the case that these hypotheses do hold. Then 8 aR = 8 b = 8 aL and the 
sum f)283p of the terms is 8p R (l) as required. That is, with the brackets we have 
adopted 0,[{p, R (l), -} , 8] = 8p R (l) for all 1 G intAf and all admissible 5 satisfying 

H = Ssf = spo = o. 

But of course this condition must hold also for 5 such that 8s R , 5s™, and 5po are 
not all zero. Considering such variations allows us to determine the brackets between 
e and Co Ri , tiim, and A. 

Let us begin with the case 1 G So. Instead of working with fi R it is somewhat 
more convenient to use the conformal 2-metric in 9 coordinates, e a b- Thus we look 
for brackets that verify the relation f2[{e a h(l), -}o,8] = 8e a b(l). 

To evaluate f2[{e a {,(l), -} , 8} we begin by noting that 

B[{e ab (l),-} o ,8] = (284) 
since e a t(l) on So commutes with e on So and with p , s R , and sl- Furthermore 
A[{e ab (l),-} ,8] 

2{e o6 (l), \} o 8p + {e ab (l), u Ri } 8s R + {e o6 (l), (b Lm } 8s m old 1 A old 2 



16nG J So 

(285) 

and we have essentially already evaluated C[{e a b(l), -} , 8}. The only brackets occur- 
ring in C[{e o6 (l), -} , 8] are of the form {e afe (l), e Ak i} with A = L or R. e ab is a 
function of p R , p R , and s R and p . Therefore, since tAki commutes with the latter 
two fields, 

{e ab {l), e Akl } = e Aki}o + complex conjugate (286) 

and thus 

C[{e ab (l), -} , 8} = [- — ]iC[{// R (l), -} , 8} + complex conjugate. (287) 
Now by flglaD 

C[{fi R {l),-} ,8] = ^ [ 8 a «p R + 8 aL p R } 1 . (288) 
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Hence 



But 



1 c) 

C[{e ab (l),-} ,S] = -[^(i a H^)Hi + complex conjugate. (289) 

l OfiR 

= \{8 aR + 5^)e ob {l). (290) 



S a e ab = 5 y e ab -5 y rd r e ab (291) 

= 5 b e ab - £ (± e ab - —[5 b - £ (± ]p d v e ab , (292) 
^Po 



and 



£t x Po = d c \po&], (293) 

with £j_ = 5s l and v = r/r the v coordinate that parametrizes the generators of the 
branch under consideration. Cl[{e ab (l), -} ,S] = 5e ab {l) therefore requires 

n r 1« l d c [p g R } a lS b p 

U — [— T.i-^n&ab + t o v -e ab — o v -e ab \i 

A 4 p 4 p 

+ / {e a 6(l),A} o 5p + {e a b(l),u;iii}o54 del A ^ 

1d7tG Jso 

+ corresponding L branch terms. (294) 



From this one obtains the brackets (I225P and (12261) . and the L analog of the latter: 

{e afc (l),A(2)} = 2nG-[d v -e ab + d v+ e ab ]5 2 (e(2)-e(l)), (295) 

Po 

{e ab (l),Mf]}o = STrG^/e^-i^-e^^]!, (296) 

A Po 

{e ab (l),uj L [g]} = SKG^geab — Tidvteab—^]!. (297) 

* Po 

Now let us consider the case 1 G N R — S again. We require 

£l[{e Rij (l),-} ,S\=5e Rij (l) (298) 



for all admissible variations 5. Reviewing the derivation of (12561) one verifies that for 
general admissible 5 

C[{e Rij (l),-} ,6] = 6e Rij (l) - ^ J e{e Rij (l), e ab } [5 a « - 8^]e ab . (299) 
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1 

16nG J So 



The other terms in Q are 
A[{e Rij (l),-} ,5] 

2{e Rij (l), \} 5po + {e Rij (l),u Ri } 8s R + {e Rij (l) : u Lm } 8s^ dO 1 A d6 2 

(300) 

and 

B[{e Rij (l) } -} ,5] 

= -y^Q J s {ZRH (!)> e° 6 }o[^- eo6 ^ + d v+ e ab 5 yL ]e. (301) 
Thus, by fl^j) . (gSHD requires 

0= / rf^A^ 2 (2{e R y(l) J A} */o 

J So 

+{e jR y(l), w Rfc } 5s^ + {e Rij (l),cj Lm } 5s™ 
+\po{e Rij (l), e ab } [£^ R - £z ±L }e ab ). (302) 



From which it follows that 

I. A(2U„ = 



{e flii (l), A(2)} = -{e^(l), e afe (2)} [^ + e ab ] 2 , (303) 



{e Rl ,(l),^[a R ]} = -- / e{e Rij (l),e ab } £^ R e ab (304) 

£ J So 

{e Rij {l),u L [£ ±L ]} = -f {e Rij (l),e ab } (e£ i±L e ab - £ ixL ed v +e ah ) (305) 

Z J So 

[Using 

£^e ab = 2D {a [e b)c e ± ] - e ab D£ c ± (306) 

(where D is the covariant derivative defined by the metric po e on So) an d t ne f ac t 
that {e R ij(l), e ab } vanishes on dSo for all 1 e intN R because <9So lies outside the 
causal domain of influence J[l] of the point 1, the last two brackets may be shown to 
be equivalent to 

36 
dy R 



{e Rij (l),u> Rk (2)} = [D a [{e Ri3 (l),e ab } o )p e bc ^\ 2 

de c 



{e Rij (l),u} Lm (2)} = [-D a [{e Rij (l),e a } o ]p e h 



1 f)f) c 
^9 c [{ ei?l ,(l),e afc } ^ + e ab ]po^] 2 . ] (308) 
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The brackets of e a b(l) and with A(2), cjr(2), and 0)1,(2) are easily obtained 
from (SOD, fMl . and (EffiHl) . They are given in (122211 . (12231) . and (12241 . But they do 
not seem to be more illuminating than the brackets found above. It is reassuring, given 
the non-obvious expressions for the brackets, that they ensure that emj commutes 
with k[v], the generator of diffeomorphisms of the 9 coordinates, as it should. 

The brackets obtained so far ensure the validity of ( 12011) for tp equal to po(l), 
Sr(1), or Sl(1) at any point 1 G intSo and for tp = e ab iX) (or any other representation 
of the conformal 2-metric) for any 1 G intM. It remains only to impose ( 120 ip for tp 
equal to A, Cjr, and Cjl- 

Let us impose 

fi[{A(l),-} ,5] = 5A(l). (309) 
The brackets that have already been determined imply that 

A[{X(1),-} ,5] = 5X(1) 

+ T7T7^ / 2 {A(1), Mo^Po + {A(l), Lj Ri } 5s R + {X(l),u Lm } 8s™ d6 l A d6 2 

l07TG J So 

(310) 

and 

B[{X(l),-} ,8] 
= -l5 m e ab d v -e ab -]-5^e ab d v+ e ab 

o o 

-JTn I W),e ab }o(d v ~e ab 5^e + d v+ e ab 6^e) (311) 
o47rG Js 

= -\^e ab d v -e ab - l -5y^e ab d v+ e ab 

o o 

+^(dv-e ab + d v+ e ab )(d v -e ab °—^ + d v+ e ab °—^). (312) 
32 po Po 

Evaluating C requires a little more thought. No derivatives of Di appear in the 
expression ( 12711) for C[5\, 62}. The bulk term in this form of C[{A(1), -} , 5} may thus 
be evaluated using the expression ( I303p for {A(l), ejnj(2)} , valid when 2 G Mr — So, 
and the analogous expression corresponding to the L branch - ignoring the jump 
discontinuity in these brackets at So. Thus 

C bulk [{X(l),-} ,5] = ^[<WliCW[{e ab (l),-} ,5] + ^[d v -e%C Lbul k[{eab{l),-}o,S\. 

(313) 

But by (1277P and (j25Tj) 

C Rbulk [{e ab (l),-} ,5] = U aR e ab , (314) 
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C Lbu i k [{e ab (l),-} ,5] = ^5 aL e ab , (315) 

so 

C bulk [{X{l), -} , S] = \d v+ e ab 5 a *e ab {l) + \d v -e ab 5 a ^e ab {l). (316) 

o o 

The Sq surface contribution to C in the form (I27ip is 

C So [5 1 ,S 2 ] = -—— I dy^AdyKn^u^ + u^n^ + dyiAdyKn^n^ + n^n^) 

07rG Jso 

(317) 

= ^ / ^i OR e o6 ^e^ + ^e a6 ^e^] (318) 

ilTXKj J So 

Now if 6i = {A(l), -} then by 



<5?e afe = {A(l), e ab } - r^{A(l), p }o d v e ab (319) 

^Po 

(because {A(l), s % } = 0). Substituting the expressions (12251) and (12131) for these 
brackets one obtains 

5t R e ab = 2nG-[d v -e ab - d v+ e ab }5 2 (6(2) - 6(1)) = -8 a ^e ab . (320) 
Po 

Hence 

C So [{X(l), -} , 5] = hd v -e ab - d v+ e ab ](5 a « - 5 a ^)e ab . (321) 
lb 

Summing all the contributions, taking into account (12921) and various cancellations, 
one finds 

fi[{A(l),.} ,a] = SX(l) + -^f {X(l),X} o 5p d6 l Ad6 2 

07TG JSo 

+1 p{A(l),^[^]}o + T ^{A(l),^[^]} 

+ T6^ e<lb ~ d - +e ^ £ ^ eab - £ ^ab)x (322) 

(since Ss^ = C ±R and Ssf = £f L ). 

The requirement (13091) thus implies that 

{A(1),A(2)} = (323) 
{X(l),u R [f]} = -nG[(d v -e ab -d v+ e ab )£ f e ab ] 1 (324) 
{X(l),uj L [g]} = -nG[(d v+ e ab -d v -e ab )£ g e ab ] 1 (325) 
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It remains only to impose 

Cl[{u R [f},-} ,6] = 6u; R [f], (326) 

and the analogous equation for cDxfg]. Our method will be closely analogous to that 
used to solve the condition ( 13091) corresponding to A(l): First note that by ( 13241) 

A[{u R [f],-} ,5] = 5u R [f] +1 f £ f e ab (d v -e ab -d v+ e ab )5e 



So 

{u R [f],u R [^ R }}o - 7^{£*[/],£L[£iiJX327) 



16ttG l "■ LJ1, ~ L ^ JJ " 16ttG 
To evaluate the B contribution recall that 5 y = 5 b — £^ ± with £j_ = 5s\ Thus 

B R [5 1: 5 2 ] = -— !— / 5 y e ab d v -e ab 6 y e-(1^2) (328) 

047TG JS 



^ / 5 l e ab d v -e ab 5 y 2 e-£^ Rl e ab d v -e ab 5 y e 

iTTLr JSo 



64vrG J So 

-^ed v -e ab ble ab + £^ R1 ed v -e ab 5 y e ab (329) 

and, by (12961) . (I2T41 . and (l292|) . 

5[{^[/],"}o,5] 

= -10Tn S &Rlf},e ab Ud v -e ab 6 y «6 + d v+ e ab 6y^) 

bAirG JSo 

+ - J (£ f ed v -e ab 5 yR e ab -£ f e ab d v -e ab 5 yR e) (330) 
= -I [ [£ } e ab -\d v -e ab ^}{d v -e ab 5 y *6-d v+ e ab 5 y ^) 

o JSo £ Po 

+\ f £ f ed v -e ab 5 a "e ab (331) 

4 JSn 



The contribution from C is obtained as in the evaluation of C[{A(1), -} , 5], by 
expanding C in bulk and surface terms according to (12711) and noting that the bracket 
and the analog of (13051) with L and R interchanged imply 



C bu ik[{u) R [f], -}o, 5] — — — I t£ fe ab C Rbu ik[{e ab , -} , 5] 

Z J So 



I So 

and therefore taking into account ( 13 14ft and ( 13151) 



+- f (e£ f e ab - £ f ed v -e ab )C Lbu i k [{e ab ,-} ,5], (332) 

Z J So 
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WML-Jc^-j / e£ f e ab (S a « - 5 aL )e ab + £ f ed v -e ab 6 aL e ab (333) 

There remains the surface term Cg [{&#[/] -}o, 8] to evaluate. Cs [Si, 62] is defined 
by (I318p . If 5\ = {ujR[f],-}o then, because lor commutes with p and sl, K l = $i 
and thus, by ([296]) 

^e afe = {wfl[/],e a6 }o = -%nG[£ f e ab - -^-^-d v -e ab ]. (334) 

* Po 

5" H , on the other hand, differs from <5j. By ( 125211 . (12T4D . and (125511 

5?*e«& = {^[/J.e^ + ieTrG^/eoft-i^^-eat] (335) 

^ Po 

= -^e a , (336) 

Thus 

C So [{Mfl '}o, 5] = ^ / 5o Kf e ^ - ^/e^-e ob ](^ - <^)e afe . (337) 
Summing this surface term and the bulk term one finds 

C[{u Ri (l),-} ,S] = -\( £ f ed v -e ab (6 a « + 5^)e ab . (338) 

o J So 

Summing A, B, and C one is left, after some cancellations, with 

n[{u R [f],-} ,5] = ^[/] + -i-i^f/],^^^ + -L^t/],^^^ 

107rG 1o7tG 

+1 f £ f e ab {d v -e ab £^ R e-d v+ e ab £ UL e) 

o J S 

-- I £ f ed v -e ab (£ (±R e ab -£ i±L e ab ). (339) 

o JSo 

Thus 

{wfl[/i],Wfl[/ 2 ]}o = -2ttG / i? /l e afe 9 1) -e a6 i? /2 e-(1^ 2) (340) 

{u R [f],u L [g]} = 2nG f [£ f e ab d v+ e ab £ g e- £ g e ab d v -e ab £ f e]. (341) 

Interchanging the roles of Mr and Ml gives us the L branch analog of (I34UI) and the 
same expression (13411) for u>£,[<7]} again. 

With these brackets fi[{<p, ■}<>,$] = Sip for all initial data <p and all admissible S. 
That is {■, -} satisfies (12011) . 
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6 Concluding Remarks 



We have found a fairly complete canonical formalism for the gravitational field on 
the domain of dependence of a single pair of intersecting null hypersurfaces. The 
phase space has been defined and a simple representation of this space in terms of 
metric fields on H 4 has been given; the Poisson bracket has been defined on a class of 
"nice" phase space functions, called "observables" , and most importantly, a bracket, 
compatible with the Poisson bracket on the observables, has been defined on a set of 
free initial data variables. 

Several things remain to be done. The most important are: 

• Resolve the problem of the Jacobi relations for the bracket on the initial data. 
The bracket on initial data proposed here does not satisfy the Jacobi rela- 
tions, which precludes a conventional Dirac quatization of the initial data. The 
requirement that the bracket on initial data reproduce the (Peierls) Poisson 
bracket when applied to the observables leaves some freedom in the choice of 
bracket on the initial data. One could try to exploit this freedom to find a new 
bracket on the initial data which does satisfy the Jacobi relations. A hope is 
that this would also lead to brackets that are simpler than the ones found here. 

I am currently working along these lines, and at the time of writing it looks 
like it will work. If, on the other hand, no such bracket exists, this would also 
be interesting. That the quantum commutators should reproduce the Peierls 
bracket on the observables seems an eminently reasonable requirement. If there 
is no corresponding bracket satisfying the Jacobi relation on the initial data, 
then this would indicate that the initial data are not quantized in the usual 
sense, and some generalization of Dirac quantization is called for. 

• Define a canonical theory for the whole Universe. The current form of the the- 
ory provides a canonical description of the domain of dependence of a single 
hypersurface A/", which never constitutes all of an inextendible solution of Ein- 
stein's field equations. One way to describe all of spacetime is to use an "atlas" 
of interrelated phase spaces corresponding to a covering of the spacetime by 
domains of dependence of the type we have considered. A state of this global 
theory would consist of points in each of a collection of phase spaces, and a spec- 
ification of how the corresponding initial data hypersurfaces, all of the form of 
A/", intersect. Of course there would then be many ways to describe the same 
spacetime corresponding to differently positioned initial data hypersurfaces, and 
one would need to define the action on the phase points corresponding to a dis- 
placement of one of the initial data hypersurfaces with respect to the others. 
This action seems to be generated by a sort of Hamiltonian, which is currently 
being studied. 
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The multi-phase space formalism outlined seems very attractive and naturally 
adapted to the local nature of general relativity. In particular it does not seem 
to impose the global causality restriction that a theory based on global Cauchy 
surfaces does. This restriction does correspond to conventional ideas about 
time, but it is alien to general relativity and perhaps simply realizes a prejudice 
we have acquired in our (nearly) Minkowskian environment. 

Of course there are a number of other, in my opinion lesser, loose ends: 

• The meaning of gauge invariance it the present canonical formalism should be 
completely cleared up. At present it is unclear to me whether diffeomorphisms 
that map the end surfaces Sl, Sr C dM to themselves, and thus permute the 
generators of Ml and Mr, should be treated as gauge or global symmetry trans- 
formations. 

• One expects that the causality condition, which we have imposed as an ansatz 
on the bracket on the initial data, really is a consequence of the form of the 
presymplectic 2-form. This should be demonstrated if true. 

• The bracket on initial data has been found assuming that the corresponding 
solution has no Killing vectors. The theory should be extended to solutions 
with Killing vectors, to the extent that this is possible. This would close a 
gap in the theory and make it possible to study the theory in concrete exact 
solutions. If no such extension is possible this would also be an interesting fact. 

Finally, we have assumed at several points that our initial data on M define 
unique maximal Cauchy developments, and that the metric of this development is 
differentiable in the parameter of any smooth one parameter family of initial data. 
To my knowledge this has not been proved, although an heuristic argument has been 
given by Sachs |Sac62j . and slightly weaker claims have been demonstrated by Rendall 
|Ren90] . and it seems virtually certain to be true. To obtain the proofs an extension 
of Rendall's method seems viable. 

The main motivation for the present work has been the desire to understand the 
holographic principle |Bec73] [tHoo93] [Sus95] . Susskind |Sus95j has given compelling, 
yet inconclusive, arguments in favour of holography based on semiclassical black hole 
theory. He also gave arguments supporting the conjecture that string theory is holo- 
graphic, using a light front formulation of string theory. My aim has been to provide 
a framework in which some of these claims could be stated precisely and proved. 

Given the present formalism (and the yet to be found bracket on initial data 
that satisfies the Jacobi relations) one is in a position to try to demonstrate a holo- 
graphic limit on the dimensionality of the Hilbert space quantizing the phase space 
corresponding to M in either a full or semi-classical quantization of the initial data. 
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To see get an idea of how the theory might be quantized one might try to see how 
(and if) the AdS/CFT duality works in terms of a null canonical formulation of the 
bulk theory. Of course this requires an extension of the formalism, to supergravity and 
string theory. Indeed the present results on gravity on a null hypersurface might give 
insight into string theory in the light front gauge, and thus help to verify Susskind's 
|Sus95] conjecture that string theory is a holographic theory. 
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A Definitions regarding manifolds and causal struc- 
ture. 

Here a number of definitions and a few propositions important to their meaning are 
gathered together for convenience of the reader. The definitions are mostly taken 
from [Wald84] or [HE73] , or equivalent to the ones given there, but the definitions of 
Cauchy surfaces and developments are slightly modified. The definitions given here 
do not form a self-contained logical structure. A complete exposition can be found in 
|Wald84j and in [HE73] . 
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A manifold will be denned as in [Wald84] p. 11. A manifold equipped with a non- 
degenerate metric will be called a metric manifold, and it will be called smooth or C°° 
if both the atlas and the metric are C°° . Unless otherwise specified all manifolds will 
be assumed to be Hausdorff and paracompact. (see |HE73j and references therein). 

A manifold, according to the (standard) definition of [Wald84j . has only interior 
points. However the concept of a manifold can be generalized to that of a manifold 
with boundary. If such a manifold is differentiable one can distinguish in the boundary 
differentiable components as well as edges, corners, and in general lower dimensional 
"strata" . 

Definition A.l A manifold with boundary is defined in the same way as a manifold 
without boundary, except that the charts, instead of being maps to open subsets o/!R n ; 
are maps to open subsets of the half space ±IR n consisting of the halfoflR n on which 
the first component, x 1 , of each n-tuple is > 0. The boundary dY of a manifold with 
boundary Y is composed of the points that are mapped to x l = by the charts. 

Differentiable manifolds generalize to stratified manifolds in which the charts may 
be maps to open subsets of the 1/4 space, 1/8 space, or l/2 m space, consisting of real 
n — tuples of which the first two, three, or m components, as the case may be, are > 0. 
The nm dimensional strata consist of the points mapped to = x 1 = x 2 = ... = x m 
in the l/2 m space by these charts. M is a stratified manifold. Note that according to 
the present definition, the boundary dY of a stratified manifold includes the union of 
all the lower dimensional strata. 

In the present work "manifold", without further qualification, will denote a C°° 
manifold with or without boundary (or lower dimensional strata). Note that a mani- 
fold without boundary can be the interior of a manifold with boundary. The absence 
of a boundary does not imply that a boundary could not be attached to the manifold, 
only that it is not included in it. 

Now we turn to definitions associated with the causal structure defined by a 
smooth Lorentzian metric on a manifold X without boundary. We shall assume that 
X is time orientable with this metric - there is a continuous choice of the future half 
of the light cone throughout X. In fact, unless otherwise specified, we shall always 
assume that Lorentzian metric manifolds representing spacetime are time orientable. 
Furthermore we shall assume that spacetime is 4 dimensional. However most defini- 
tions and results generalize in an obvious way to spacetimes of any dimensionality. 

Definition A. 2 The causal future, J + [S;X], of a set S C X in X is the set of 
points of X that can be reached from some point in the set S via a future directed 
differentiable causal curve in X - a curve with tangent everywhere null or timelike. 
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Definition A. 3 The chronological future, I + [S; X], of a set S C X in X is the set of 
points of X that can be reached from some point in S via a future directed differentiable 
timelike curve in X of non-zero duration - a curve with tangent everywhere non-zero 
and timelike. 

Definition A. 4 The future domain of dependence D + [S;X] of S C X in X is the 

set of points p G X such that all inextendible causal curves in X from p intersect S 
to the past of p. 

(Often, when the ambient spacetime manifold with respect to which these sets are 
defined is clear from the context, it will not be specified explicitly, and we shall write 
I+[S], J+[S], and D+[S}.) 

The definitions of the causal past J~ , chronological past I~ , and past domain of 
dependence D~ are the obvious time reversed analogs of those of J + , I + , and D + . 

Uniting the past and future parts one may define the full domain of dependence 
as D[S] = D + [S] U D~[S], the causal domain of influence as J[S] = J + [S] U «/~[iS], 
and the chronological domain of influence as I[S] = I + [S] U D[S] may be 

defined directly as the set of points p G X such that all inextendible causal curves in 
X through p intersect S (to the past or to the future of p). While J[S] and I[S] are 
the unions of the causal and timelike curves through S respectively. 

Definition A. 5 A subset S of X is achronal iff no timelike curve in X crosses S 
more than once, i.e. at two or more values of the curve's parameter. 

Equivalently, S is achronal iff S, I + [S], and I~[S] are disjoint. 

The achronality of a set requires that its interior, suitably defined, be an embedded 
3-manifold: If S is achronal then S is disjoint from -/ _+ [S'], but of course S C 
Thus S C i + [S). But by prop. 6.3.1. of [HE73] I + [S] is a C^H embedded 3- 
manifold. If the interior of S, intS, is defined to be the topological interior of S 
in I + [5'], then intS is also a C l ~ , embedded 3-dimensional submanifold (without 
boundary). Note that intS is not the interior of S in the ambient spacetime X, 
which is empty. 

This makes possible a more intrinsic definition of intS: 

Definition A. 6 The interior, intS of an achronal subset S C X is the largest subset 
of S that is an embedded 3-submanifold of X without boundary. (Equivalently it is 
the union of all such subsets). 

The equivalence of this definition with the earlier one can be deduced from the 
following topological lemma: 

72 It satisfies a Lipschitz condition, and a fortiori is C°. See |HE73j . 
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Lemma A.l If A is an n-manifold with boundary embedded in a n-manifold without 
boundary B, then A — dA, the interior of A according to its structure as a manifold 
with boundary, equals A — A, the interior of A in the topology of B, and dA C A. 

Proof. A point p G A lies in the domain of a chart ip : V — >§IR n of A and a chart 
: V — > IR n with V and V open neighbourhoods of p in A and B respectively. If U 
is an open neighbourhood of p in A which is contained in both V and V then (poifj~ l 
is a homeomorphism from i/j[U], an open subset of |H n , to the subset <j)[U] of H n . 

Suppose p G dA. ip[U] C|lR n cannot then be open in JR n , because it contains ^(p) 
which lies on the boundary of ±lR n in the topology of IR n . On the other hand, since A 
is an embedded submanifold with boundary the open sets of A are the intersections 
of A with open sets of B. Thus there exists a set U' open in B such that U — U' DA. 
4>[U'] is open in H n . Brouwer's invariance of domain theore m | Brol2| requires that 
the homeomorphic image in IR n of an open set in IR n be openo The open set 4>[U'] 
thus cannot be homeomorphic to ip[U], and therefore not to <j)[U] = <f)[U' R A] either. 
U must contain a point outside A. It follows that p G A, and hence that dA C A. 

Now suppose that, on the contrary, p G A — dA. Then ip(p) does not lie in the 
boundary plane of ±IR™ and there exists a neighbourhood R of ^(p) in ^[Vfl V] C§IR n 
which does not intersect the boundary of |IR™ and which is open in ±IR n and also 
in H n . Define U = il)~ l [R\. U is a neighbourhood of p open in A and a subset of 
V D V . Since R — ip[U] is open in IR", <f>[U] is also. U is thus open in B, implying 
that p lies in the interior A — A of A according to the topology of B. Together with 
the previous result that no point in dA lies in A — A this implies that A — dA = A — AM 

From the lemma it follows that any subset U C S that is a 3-manifold without 
boundary (dU = 0) must coincide with its interior in the topology of I + [S] , and thus 
be open in this topology. This establishes the equivalence of the two definitions of 
intS. 

A useful notion of boundary for S is the following: 

O — 

Definition A. 7 The embedding boundary of an achronal subset S C X is S= S — 
intS, where S is the closure of S in the topology of X. 

o . o 

S is simply the topological boundary of S in Note that S does not in 

general coincide with the causal edge of S defined in |HE73j p. 202. It is easy to show 

that SC. edge[S}. However not all points of the edge need to lie in sE 

73 A homeomorphism maps the open sets of its domain to open sets in its range. This does not 
by itself imply that the range is open in the topology of an ambient space whenever the domain is 
open in an ambient space. 

o , , 

74 S does equal edge[S) if strong causality holds. See |HE73j p. 196-197 on Alcxandrov topology. 
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Lemma I A. 1 1 provides us with a further useful and obvious corollary. 



Proposition A. 2 If is a 3-manifold with boundary achronally embedded in space- 
time then intH = X — <9£, the interior o/E according to its manifold structure, and 

<9£ CE. 

o 

Conversely, if the embedding boundary, S, of an achronal set S C X is embedded 

— — o 

in X then S is a manifold with boundary, with dS =S- [If a manifold is embedded 
according to our definition (that of |HE73j ). which is to say regularly embedded in 
the definition of [CDD82j . then according to |CDD82j p. 242 it is a submanifold of 
the ambient manifold according to the definition of [CDD82j p. 239. The result then 

_ O 

follows immediately] Note that if X is closed prop. IA.2I implies that £= <9E. 
Now we turn to Cauchy surfaces and Cauchy developments. 

Definition A. 8 A Cauchy surface of an open set A C X is an achronal subset C of 
the closure A of A such that any causal curve in A which is inextendible in A crosses 
C or has an endpoint on C in A. 

Definition A. 9 An open subset Ti, C X is globally hyperbolic in X iff Ti contains 
a Cauchy surface ofTL. 

Whether H is globally hyperbolic in X depends to some extent on the properties 
of the ambient spacetime X. Cauchy surfaces must be achronal, and whether a given 
subset is achronal depends on the manifold in which it is embedded. However several 
important consequences of global hyperbolicity hold for a subset Q provided only that 
it is globally hyperbolic regarded as a manifold by itself, i.e. provided that it's Cauchy 
surfaces are achronal in Q. For instance, this weaker property, which might be called 
intrinsic global hyperbolicity, is sufficient to establish the existence of a foliation of Q 
by homeomorphic spacelike hypersurfaces [HE73J prop. 6.6.8. 

Def. IA.9I of global hyperbolicity has been adopted because it is equivalent to the 
definition used in |HE73j . There a subset Ti of X is globally hyperbolic in X iff strong 
causality holds on the subseiQ and J + {p) PI J~(q) is compact and contained in Ti for 
all p,q G H. When H is open this can be shown to be equivalent to the definition 
IA.9I using Geroch's theorem prop. 6.6.8. |HE73j and prop. 6.6.3. of |HE73] . (See 
|Wald84j p. 209.) 

The following result, which extends the relation between Cauchy surfaces and 
global hyperbolicity established by def. IA.9I to Cauchy surfaces outside the globally 
hyperbolic set, justifies our definition IA.8I of Cauchy surfaces. 

75 Strong causality holds at p if every neighbourhood of p contains a neighbourhood which causal 
curves cross at most once. See [HE 73] . 
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Proposition A. 3 An open subset H of X is globally hyperbolic iff intTL = H, and 
7i posses a Cauchy surface. 

Proof. Let C C 7i be a Cauchy surface of 7i. 7i is therefore a subset of intD [S] which, 
by [HE73j Prop 6.6.3, is globally hyperbolic. (S is achronal because S is.) Thus strong 
causality holds at all points of H and J + (p) fl J~(q) is compact Vp, q G H. All that 
remains to show is that J + (p) fl J~(q) C 7i. 

Let p' G J~(p) n and g' G I+(g) n 7i, then J+(p) n J _ (g) C I + (p') n I~{q'). It 
is therefore sufficient to show that I + {p') fl I~{q') C Vp', g' G 7Y. 

Suppose this were not so and J + (p') fl I~(q') contains a point r outside 7i. Let 7 
be an inextendible timelike curve through p', r, and q' . Since q 1 and p' lie in 7i but 
r does not, 7 fl H consists of at least two segments that are inextendible in 7i, one 
before r, the other after. Each segment must cross C or have an endpoint on C in 
H, so 7 crosses C at least once at r or later (according to the parameter of 7), and 
at least once at r or earlier. On the other hand C is achronal, so 7 can in fact cross 
C only once, implying that the crossing occurs at r, and therefore that r G C C 7i. 
We may conclude that I + (p') fl I~(q') C 7Y. But I + (p') H I~(q / ) is open, so it must 
actually be contained in mt7Y = 7Y. 

Conversely, if is globally hyperbolic it contains a Cauchy surface. What must 
be shown is that infH = H, in this case. Suppose p G intTL, then 7i is dense in 
any sufficiently small open neighbourhood U of p. Let p' G I~(p) D £/ H 7i and 
g' G (1 f/ fl H, then p G I + (p') fl I~{q'). The argument used to show that 

I + (p') fl I~(q') C 7i applies, but since the Cauchy surface now is contained in 7i one 
may conclude that I + {p') fl I~(q') C and thus p G 7Y. ■ 

Prop. IA.3l is essentially a generalization of prop. 6.6.3. of [HE73j which states that 
the interior of the domain of dependence of a closed achronal set is globally hyperbolic. 
Here we reproduce this very useful result, and also show that the assumption that 
the achronal set be closed is unnecessary. 

Proposition A. 4 The interior of the domain of dependence of any achronal subset 
of X is globally hyperbolic in X . 

Proof. The claim can be deduced from prop. IA.31 but it is easier to prove it directly. 
Suppose S C X is achronal and let B = intD[S], B is a subset of intD[S] which, by 
|HE73j Prop 6.6.3, is globally hyperbolic. Thus strong causality holds at all points 
of B and J + (p) fl J~(q) is compact Vp, q G B. All that remains to show is that 
J + (p)nJ'(q) C B. 

Let p' G J-(p) H B and q' G I + {q) fl B, then J+(p) n J~(q) C I + {p') fl /"(</)• An Y 
point r G I + {p') fl I~(q') lies between p' and g' on a timelike curve. If this curve is 
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extended until it is timelike and inextendible then it must cross S, and it must do so 
precisely once since S is achronal. Any point between p' and the crossing point, or 
between q' and the crossing point, must lie in D[S] so r G D[S]. Since I + (p') fl /"(</) 
is open it must in fact be contained in the interior, £>, of -DfS 1 ]! 

A peculiar feature of our definitions is that while the interior of the domain of 
dependence of an achronal set S is globally hyperbolic S is not necessarily a Cauchy 
surface of intD[S], because it may not be contained in intD[S}. For instance suppose 
a is a spacelike disk in Minkowski space, and let S = I~[o], S H intD[S] = a. 
Nevertheless, it will be proved further on (prop. IB. lip that M is a Cauchy surface of 
V = intD[M]. 



Definition A. 10 The future Cauchy horizon of a set S is H + [S] = D[S] - I~[D[S\] 
The past Cauchy horizon is defined analogously. 

We may now give the definition of Cauchy developments used in the present work 

Definition A. 11 A smooth Cauchy development of initial data specified on a 3- 
manifold E is the domain of dependence of 'E embedded achronally in a C 00 ]^ time 
orientable, Lorentzian metric J^-manifold without boundary such that the data specified 
on E match those induced on E by the J^-metric. Moreover the 4-metric is required to 
satisfy the field equations in the domain of dependence. 

In general Cauchy developments defined in this way are not manifolds (with or 
without boundary). Nevertheless some work might yield a more intrinsic definition 
which does not refer to the ambient boundaryless manifold in which the domain of 
dependence exists. On the other hand we make use of this ambient manifold on 
occasion, so the definition seems adequate for our purposes. 

Given initial data on a given 3-manifold can have many different Cauchy develop- 
ments. For instance if some of the points in the domain of dependence are eliminated 
from the ambient manifold the new domain of dependence in the manifold thus re- 
duced would be a proper subset of the original one. More generally for any given 
Cauchy development there are many others that are isometric to a proper subset. 
Moreover the Cauchy development might be isometric to a proper subset of a larger 
one. This leads naturally to the concept of a maximal Cauchy development. 

Definition A. 12 A maximal C°° Cauchy development is a Cauchy development 
which is not isometric to a proper subset of any other C°° Cauchy development of 
the same data on the same 3-manifold. 



76 If no hypothesis is made about the degree of differentiability of the 4-metric then it becomes 
difficult to define what is meant by a solution to the field equations. 
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Zorns lemma asserts that maximal Cauchy developments exist, while uniqueness 
theorems assure that they are uniquely determined up to diffeomorphisms by the 
initial data, for the classes of initial data to which these theorems apply (see |Wald84j 
p. 263-264). 

B Some facts about the null initial data hypersur- 
faces used in the present work 

This appendix collects together results about the topology, geometry, and causal 
structure of the initial data hypersurface Af (defined in subsection 12. ip and its Cauchy 
developments. 

We shall begin with the local aspects of the geometry of Af. It is shown that 
the tangents of the generators of Af are normal to Af. As a consequence all two 
dimensional cross sections of Af that are transverse to the generators are spacelike 
and normal to the generators. A converse result is also given that establishes that 
if the normal vectors of a hypersurface are null, making it a null hypersurface, then 
these normals are tangent to the hypersurface and integrate to null geodesies that 
sweep out the hypersurface. 

Next we turn to the issue of caustics and crossings of generators. Typically the 
generators, if extended far enough, will form caustics. But (see fig. [3]) a generator may 
well cross another generator before reaching a caustic and thus enter the chronological 
future of Af. If this occurs Af is not achronal. The question therefore arises whether 
it is really true that there always exists a solution in which Af with given data is 
achronal, or whether the initial data must satisfy some additional conditions in order 
for such a space to exist. Moreover, if there is no such condition, and any data 
corresponds to an achronal embedding, can a solution in which Af is not achronal be 
represented by a Cauchy development of initial data? 

In a Cauchy development Af is achronal by definition, but we have at present 
no proof of the existence such developments. There is only a theorem showing the 
existence of Cauchy developments of a neighbourhood of Sq in Af, and we have been 
assuming that developments of all of Af exist. In the following it will be shown, quite 
independently of any existence theorem for Cauchy developments, that if Af is not 
achronal in a spacetime X then there exists a locally isometric covering manifold V 
of a neighbourhood of Af in X in which the future directed null geodesies normal to 
So sweep out an achronal hypersurface covering Af. This answers both questions of 
the previous paragraph: Achronality of Af requires no additional restrictions on the 
initial data (beyond those that ensure the absence of caustics on Af). Furthermore, 
if we suppose that initial data determines unique maximal Cauchy developments so 
that (portions of) solution spacetimes can be represented by corresponding initial 
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data, then solutions in which M is not achronal are represented by locally isometric 
covering manifolds with geometry corresponding to the initial data on N '. Therefore 
such solutions are represented by initial data, but the identifications that differentiate 
the covering manifold from the original solution spacetime, and make M non-achronal, 
leave no trace in the initial data. 

Adopting the covering manifold as our spacetime we can be certain that M is 
achronal. Using this we prove a number of further results: Points on M are causally 
related iff they lie on the same generator; All Cauchy surfaces of D[M] share the same 
boundary dM; and others. 

Finally two useful results regarding Cauchy developments of M are established: 
Firstly, that any Cauchy development of M may be extended to a boundaryless glob- 
ally hyperbolic manifold containing M, and secondly, that all Cauchy developments 
of M, with all possible data, are diffeomorphic to each other (though of course not 
necessarily isometric) . 

Let us first revisit the definition of M. Suppose Sq is smoothly embedded as a 
spacelike 2-surface in a 4-manifold X with smooth Lorentzian metric and let O be 
the normal bundle of Sq (i.e. the bundle of normal vectors). This bundle contains the 
subbundles Nl and Nr of future and inward directed null normal vectors and future 
and outward directed normal null vectors respectively (inward and outward referring 
to the arbitrarily chosen orientation of So). Ml and Mr are obtained by suitably 
truncating each ray of Nl, and Nr and applying the exponential map. Here a ray is 
the set of positive multiples of a vector, and the exponential map is the map that 
takes a point p E X and a vector v in the tangent space of X at p to the point of affine 
parameter 1 on the geodesic with tangent v originating at p. The exponential map 
is well defined as long as v is not so large that the geodesic leaves X before reaching 
affine parameter 1, and then it is a smooth map because the metric is smooth (see 
|HE73j p. 33). Restricted to the normal bundle O it defines a smooth map / from 
a neighbourhood of Sq in O into X. The image of the neighbourhood is not open in 
X but a slight modification of our setup makes it so. Since So is smoothly embedded 
it may be extended to a (boundary-less) spacelike 2-submanifold S' D So of X. [The 
smoothness of the embedding and the Whitney extension theorem |AR67j imply that 
the embedding map defined on So, diffeomorphic to the unit disk in H 2 , may be 
extended to a boundaryless manifold containing Sq. An open subset S' D So of this 
extension will be spacelike.] / extends to the normal bundle O' of 5" and maps an 
open neighbourhood of So (or indeed of S') in O' to an open neighbourhood of the 
same 2-surface in X. 

f is not necessarily globally invertible, but if it is restricted to an open subset V 
of its domain in O' that excludes all points at which the Jacobian df is degenerate, 
then by the inverse function theorem 
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Proposition B.l / : V — > X is a smooth immersion ofV in X. That is, for each 
point P e V there is a neighbourhood U of P such that f restricted U is a smooth 
diffeomorphism of U to f[U]. 

Note that the images in X of points in O at which df is degenerate are caustics 
of the congruence of geodesies normal to Sq. That is, they are points at which 
neighbouring normal geodesies cross, or nearly cross in the sense that the coordinate 
separation of points of equal affine parameter on a family of nearby normal geodesies 
is zero to first order in the difference between their initial data - in the differences 
between the coordinates of their starting points on S and of the components of their 
initial tangents. 

The truncations, Nl and Nr, of Nl and Nr which, when exponentiated, give rise 
to Ml and Mr are chosen so that there are no caustics on Ml or Mr. V may thus be 
chosen to so that it contains Nl and Nr. V will furthermore be chosen so that any 
ray of O' leaves V at most once. It is then simply connected, prop. IB. II implies that 
V is a covering manifold of f[V] C X. Let us denote by V the manifold V equipped 
with the metric pulled back from X via /. V is then a locally isometric covering 
manifold of f[V]. The exponential map / : V — > V of V into V is simply the identity 
map from V to V and thus evidently one to one. 

In the following we will mostly study the hypersurface M = f[N L U Nr] swept 
out by the future directed null geodesies normal to So in the covering spacetime V. 
Since this M maps to the old M in X when the isometric identifications are made 
that reduce V to f[V] C X, any local property of M in V is shared by TV in X. In the 
corresponding Cauchy problem data is given on N = Nl U Nr and the embedding of 
TV in the Cauchy development is diffeomorphic to M in V. 

Note that in our original definition of M we only required that the generators of 
M not form any caustics on M. Here we require that the congruence of geodesies 
normal to Sq, which includes also timelike and spacelike geodesies, not form caustics 
on M . In prop. IB.4I the two conditions will be shown to be equivalent. 

Now let us establish the key local property of M, namely that it is a null hyper- 
surface and the null normal vectors are the tangents to the generators: 

Proposition B.2 The tangents to the generators of a branch Ma of M are normal 
to M A - 

Proof. Suppose ( is a tangent vector of Ma at a point p. Let A be an affine parametriza- 
tion of the generators of Ma such that the corresponding field of tangent vectors 
I = d/dX of these is smooth (this requires only that d/dX be chosen smooth on So), 
and define the vector field r\ on the generator through p by Lie dragging ( along the 
generator. That is, define rj to be the field of vectors tangent to Ma that satisfies the 
differential equation £\r\ = and the initial condition rj(p) = C- 
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Since t] is tangent to Ma a t 5*0 it can have only components tangent to So and to 
/. This implies that I ■ r\ — there, for the component tangent to So is orthogonal to / 
by the definition of the generators, and the component tangent to I is orthogonal to / 
since I is null. But I ■ r\ is constant along the generator through p: £in = means that 
V;?7 = and, since I is the tangent to an affinely parametrized geodesic, V;/ = 0, 
so 

-jL[l. v ] = l-V l ri = l-V v l= \^l 2 = 0. (342) 

Thus I • j] — at pM 

A converse result is worth noting. 

Proposition B.3 On an arbitrary null hypersurface JC (a hypersurface with null nor- 
mal everywhere ) in a smooth Lorentzian metric manifold Y the normals are tangent 
to the hypersurface and 

integrate to a congruence of null geodesic curves that sweep out the hypersurface. 

Proof. Let w be a function on Y that is constant on JC and has non-zero gradient 
there. Then k = dw, obtained by acting with the inverse metric on dw, is a normal 
vector field on /C, which by assumption is null and thus has null integral curves. The 
torsion freeness of the connection implies that these integral curves are geodesies: 

= 2hPV ]li V v] w = W k k v - VW u k^ = V k k v - -V^ 2 , (343) 

so Vfc/c oc k, since k 2 = on JC. That is, k is parallel transported, modulo modulus, 
along it's integral curves which are therefore geodesies.! 

Note that any 2-surface o C JC that cuts the integral curves of the normals is of 
course orthogonal to these. (This means that in any local Lorentz frame in which 
an element of a is purely spatial the integral curve emerging from this element is 
spatially normal to it.) 

Note also that if the level sets of w are null hypersurfaces in a neighbourhood of 
JC, or indeed if only the gradient of (dw) 2 = k 2 vanishes on JC, then it follows from 
( 13431) that the integral curves of k are affinely parametrized geodesies. 

We are now ready to show that the restrictions we place on the initial data imply 
that df is non-degenerate on N , which is to say, the generators of M and neighbouring 
geodesies normal to Sq form no caustics on M. What are these restrictions? Suppose 7 
is a generator of a branch Ma of M, and thus the image of a ray of N, and suppose P is 
a point on this ray. Then the area parameter at P is defined to be r(P) = Ja(P)/A, 
the root of the cross sectional area A at f(P) of an infinitesimal bundle of generators 
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neighbouring 7, over the cross sectional area A of the same bundle at the truncating 
surface Sa- (See subsection 13. ip . In the initial data sets we admit the area parameter 
along each ray starts at a finite, positive, non-zero value ro at the beginning of the 
ray on So and increases or decreases monotonically to a value of 1 at the truncating 
surface where the ray leaves N. Thus a bundle of generators having non-zero cross 
sectional area on Sq necessarily has non-zero cross sectional area everywhere along 
its trajectory to Sa- 

This restriction on the initial data together with the following proposition implies 
that df is non-degenerate on N. 

Proposition B.4 df is non-degenerate on N provided r/r does not vanish on J\f. 

Proof: If df is degenerate at a point P G O then / maps any basis of the space Tp 
of vectors tangent to O at P to a linearly dependent set of tangent vectors to X at 
p = f(P) G X. To show that df is non-degenerate on a branch Na of N we will choose 
a basis of Tp for all P G Na and show that the 4-volume spanned by the image of 
this basis is non-zero, implying that the four image vectors are linearly independent. 

Coordinates can be defined on Na as follows: Choose coordinates on Sq and 
convect these along the rays to all of Na (i.e. set them to be constant on the rays), 
then choose a linear parameter A on each ray so that A = on Sq and the A = 1 
surface is smooth. The associated coordinate basis consists of L = d/dX, tangent to 
the null rays of Na, and two vectors, E\ and E 2 , corresponding to the coordinates on 
Sq. This can be completed to a basis of tangents to O on Na by including a timelike 
vector T which is tangent to the fibres of O and constant in each fibre. (T is timelike 
according to the metric on the tangent space of X at the base of the fibre on Sq where 
T can be regarded to reside since it is constant on the fibre). At each point P G Na 
f pushes these vectors forward to vectors I, ei, e 2 , and r tangent to X at p = f(P). 
I is tangent to the generator through p and thus null and normal to Na, while e\ and 
€2 are tangent to Ma and thus normal to I. Since = 1-1 = 1- e\ = I ■ e 2 the 4x4 
matrix of inner products of these vector^] falls into a block diagonal form and the 
4-volume spanned by them is easily seen to be — / ■ tA, where A is the area spanned 
by t\ and e 2 . 

Since L and T are non-zero null and timelike vectors respectively the inner product 
/ • r is necessarily non-zero at So. Less obviously, this inner product is constant along 
each generator: The displacement from a point on a ray of Na, generated by L, to the 

77 When I, a, e2, and r are linearly independent, and thus form a basis, the matrix of their inner 
products is the matrix of covariant components of the metric referred to this basis. The matrix 
of inner products is of course defined even when these vectors are not linearly independent. The 
4-volume is given by the square root of minus the determinant of the matrix of inner products in 
cither case. 
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corresponding point on the timelike ray of O generated by the vector L + T is given by 
the vector field V(A) = AT on the null ray. Ar is thus a geodesic deviation vector field, 
or Jacobi field, and it is subject to the geodesic deviation equation: V;V;At m = A-R^. 
Since V;/ = the inner product XI ■ r satisfies 



which shows that XI ■ r is linear in A, and hence that I • r is constant along the 
generator. In conclusion / • r is non-zero everywhere on the generator. 

Degeneracy of df on would therefore require that the area A spanned by e\ 
and 62 vanish somewhere on a generator. But A = Ao(r/ro) 2 , where Aq is the area 
spanned by these vectors at So, which is non-zero by definition. The condition that 
r/Yo not vanish thus guarantees that A does not vanish and df is non-degenerate.l 

Now let us show that M is achronal in V. This of course does depend on the fact 
that M is embedded, and not merely immersed, in V. The square s — v ■ v (according 
to the metric on S') of the vectors orthogonal to S' constitutes a smooth function 
on O', and thus on V. At a point p G V s is minus the proper time squared elapsed 
along the orthogonal geodesic segment 7 P from 5" to p, or, if 7 p is spacelike, it is the 
length squared of j p . 

s is also the value on these geodesic segments of the action functional a[c] = 
Jq 1 v ■ v dt, when t is an affine parameter. This action functional is defined for an 
arbitrary differentiable curve c : [0, 1] — > V, with t the parameter of the curve and 
v = d/dt the corresponding tangent vector. When the curve is the geodesic segment 
7 p normal to 5" parametrized aflinely so that c(l) = p then, firstly, v(0) is the vector 
at c(0) G S' that the exponential map associates to p, and, secondly, v is parallel 
propagated so v ■ v is constant on the geodesic. As a consequence a = s. 

Using this correspondence we may prove the following useful lemma: 

Lemma B.5 The level set of s passing through p is orthogonal to j p . Indeed the 
gradient vector ds of s ( obtained by acting with the inverse metric on ds ) is tangent 
to 7 P and points away from S' . 

Proof. A direct calculation shows that under a variation of c 



{8c is the vector field on the curve generating its displacement). This allows us to 
evaluate the gradient of s. Let us restrict attention to the geodesic segments 7 P 
and vary p. The base point of 7 P on S' varies, but according to (13451) . since v(0) is 



dX 2 



X[l-r]=l- V^V/Ar = XR Tlll = 0, 



(344) 




(345) 
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orthogonal to 5" this variation makes no contribution to the variation of a. Since 7 P 
is a geodesic there is no contribution from the integral term either. The variation of 
the action thus is 8a = 2<5c(l) ■ v(l) 

ds = 2v(l). (346) 
From this the claim follows immediately.! 

Now let us use this lemma. Let V + be the subset of V consisting of future directed 
causal vectors. It follows from lemma [BT5l t hat s is strictly decreasing along any future 
directed timelike curve in V + = /[V + ]. From this we see immediately that 

Proposition B.6 M is achronal in V. 

Proof: If M were not achronal then there would exist a future directed timelike 
curve in V crossing M twice. Moreover, by |Wald84j theorem 8.1.2., this curve may 
be chosen to be piecewise geodesic, and thus piecewise differentiable. [An arbitrary 
timelike curve may be covered with a sequence of convex normal neighbourhoods so 
that successive neighbourhoods overlap. Then a sequence of points may be chosen 
lying on the curve and in the overlaps of successive neighbourhoods and successive 
points in this sequence may be connected by timelike geodesies, which always exist 
by the theorem.] s is strictly decreasing along the curve at any point in V + — Sq. 
The curve starts on M C V + where s = 0, and if it starts on S it immediately enters 
V + — So because it is timelike and future directed. Since the curve can leave V + — So 
only if s returns to 0, its value on the boundary of this set, the curve can in fact not 
leave V + — So and s always continues to decrease strictly after the starting point of 
the curve. An endpoint of the curve on M would therefore have s < 0, but s = on 
all of M. I 

Prop. IB.6I establishes that the achronality of M requires no restrictions on the 
initial data, save that it make M caustic free. The absence of caustics implies that M 
can always be made achronal by going to the covering manifold V of a neighbourhood 
of A". 

V does not necessarily cover all of the domain of dependence of A/", even in ambient 
spacetimes in which M is achronal. A caustic in the timelike geodesies normal to So 
might well form inside -D[A/] without there being any sort of pathology in D\M] 
that one would wish to exclude. Thus V is not in general large enough to be the 
ambient spacetime of the maximal Cauchy developments of data on A/". Nevertheless, 
the domain of dependence of M in V, while not the maximal Cauchy development 
is sufficient for our purpose of defining a pre-Poisson bracket on the initial data. 
One may define the observables on these, non-maximal, Cauchy developments, and 
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the Peierls bracket on these determines an almost unique pre-Poisson bracket on the 
initial data. 

For the sake of completeness, and conceptual simplicity, I define further on a 
boundaryless manifold M which contains both M as an achronal submanifold and 
the interior of the maximal Cauchy development of the data on M. This construction 
is not necessary to justify the main results of the present work. 

We are now in a position to show that 

Proposition B.7 In V the subset of M that is causally connected to a point p G M 
(i.e. lies in J + [p] U J~[p\) consists of the points on the generators through p. 

Proof. Suppose q and r are points of M that are causally connected, with q to the 
future of r: q G J + (r). Let M e be the extension of M swept out by the generators of 
Af continued to the future as null geodesies until they leave V. s = on M e so M e is 
achronal in V by the same argument used to prove the achronality of M (prop. IB. 61) . 
Now suppose q' lies on the generator through q to the future of q, then a causal curve 
7 from r to q followed by the generator segment from q to q' defines a causal curve 
from r to q' which may be deformed to a timelike one unless it is a null geodesic. The 
achronality of Af e thus implies that 7 is a null geodesic that continues the generator 
segment from q to q' to the past, that is, it implies that 7 coincides with the gen- 
erator through q and q' . It follows that r and q necessarily lie on the same generator. I 

A further property of Af that we will use is 

Proposition B.8 If W is a neighbourhood of dAf in Af then there exists a neigh- 
bourhood U of the boundary in So such that all generators of Af incident on U lie 
entirely in W . 

Proof. Consider one branch Ma of M. The complement C of W in Ma is a closed 
subset of Ma, and Ma is compact. Thus C is compact. The map n : Ma —> So which 
takes the points on a generator to the origin of that generator on So is continuous 
- it is just the composition of the diffeomorphism f~ l (restricted to Ma) with the 
projection to the base manifold in O. Thus n[C] is compact and therefore closed. It 
follows that the complement Ua of tt[C] is open. Since C does not intersect dM and 
7r _1 [<9So] C dM, 7r[C] does not intersect dSo. Thus Ua is an open neighbourhood of 
dS . The same is of course true of U — UrDUlM 

Now we turn to properties of Cauchy developments of M. We begin with a pair 
of lemmas about achronal sets. 

o 

Lemma B.9 If S is an achronal set then its embedding boundary S lies outside the 
domain of influence I[D[S]] of the domain of dependence of S . 
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Proof. Since I + [S] is achronal (by [HE73j prop. 6.3.1.) a timelike curve can cross 
I + [S] at most once. On the other hand all inextendible timelike curves through D[S] 
must cross S C They can thus not cross I + [S] outside S as well. Indeed no 

timelike curve through D[S] can cross I + [S] — S since any such curve could be ex- 
tended until inextendible. /[-DfS 1 ]] n/ + [S'] is thus a subset of S, and since it is open in 

o 

the topology of I + [S] it must be contained in the interior of S and be disjoint from S 1 .! 

Lemma B.10 Suppose S is an achronal set and p £ intS , then I + (p) H D[S] = 
iff there exists a null geodesic in S passing through p without past endpoint in the 
interior of S . 

Proof. If there exists a null geodesic segment 7 in S passing through p without past 
endpoint in intS then either this segment is past inextendible in spacetime or it has 

— o 

a past endpoint s in the closure S of S, and thus in S — intS =S- In the first case 
there exists (by [Wald84] lemma 8.1.4.) through any q £ P^{p) a past inextendible 
timelike curve contained in I + [7] C I + [S] . Since S is achronal this curve never crosses 
S, implying that q £ D[S\. In the second case one may remove the end point s from 
the spacetime making 7 past inextendible. In the spacetime without s there again 
exists a past inextendible timelike curve through q contained in J + [iS]. This curve is 
either inextendible in the full spacetime, indicating that q ^ D[S], or s G S is its past 

O 

endpoint. But by lemma IB~9l no timelike curve through D[S] passes through S, so 
also in this case q D[S}. I + (p) is thus disjoint from D[S}. 

Now to the converse: if I + {p) is disjoint from D[S] then there exists a sequence 
of points qi in I + (p) approaching p and causal, past inextendible curves 7, through 
the qi which do not intersect S. The portions, 7-, of these curves that lie in I + [S] are 
either past inextendible or have past endpoints on I + [S] — S. Removing the closed 
set I + [S] — intS from spacetime leaves a manifold in which the 7- are all past inex- 
tendible. By lemma 6.2.1. of [HE73j there must exist a causal, past inextendible limit 
curve of the 7- through p. Since the 7^ all lie in the closed set I + [S] the limit curve 
does so also. But, because S is achronal the limit curve cannot enter It is thus 

confined to I + [S] — I + [S] = I + [S], which is just intS in this manifold. Moreover the 
limit curve, which is causal, must be a null geodesic, for if it were not then two points 
on the curve could be connected by a timelike curve, violating the achronality of S. I 

Note that when I + (jp) n D [S] is non empty for some p £ S then in fact all of I + (p) 
that lies within a neighbourhood U of p is contained in D[S}: If r £ I + (p)r\D[S] then 
we may take U = I~(r), because p £ I~(r) and past directed timelike curves from r 
that have not yet reached S must be contained in D[S], so I~{r) n I + {p) C D[S]. 
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Lemma [B. 101 allows one to demonstrate some intuitively plausible basic facts about 
the domain of dependence of M '. 

Proposition B.ll For all p G intM ', the part of I + (p) lying in a sufficiently small 
neighbourhood of p is contained in D[M] while I~(p) does not intersect D[M]. Thus 
D[Af\ = D + [M], intD[M\ is non empty, M C intD[M], and M is a Cauchy surface 
ofintD[M]. 

Proof. By prop. IB. 71 the only null geodesies in M passing through a point p G intM 
are the generators through p. These have past endpoints on intSo C intM, and future 

O 

endpoints on Sl or Sr, both subsets of dM cM- Thus by lemma IE. 101 I + (p) must 
intersect D[M], while I~{p) must not intersect D[M]. That U D I + (p) C D[M] for 
some neighbourhood U of p follows from the remarks after the proof of that lemma. 
D[M]= D+[M], intD[M] ^ 0, and M C intD[M] follow immediately. M C intD[M] 
implies that M is a Cauchy surface of intD[M] according to def. IA.81 I 

Since the interior of the domain of dependence of the initial data hypersurface M 
we are using does not cover the whole spacetime, and in fact does not even include M 
itself, it has been necessary to adopt a wider notion of Cauchy surfaces than is usual 
(see def. IA.8I) . In the following we shall see that despite this more ample definition, 
the Cauchy surfaces intD[M] are quite simple: They are all embedded manifolds 
homeomorphic to M with boundary coinciding with dM. 

Proposition B.12 If C and C are both Cauchy surfaces of an open set H C X, 
then C=C 

O — o 

Proof. C is the boundary of C in / + [C] and also of its complement in / + [C]: C= 
[i + [C)-C}°. But 

I+[C] -C = i + [H\ -H, (347) 

the left side being a set depending only on 7i and not on the choice of Cauchy surface. 
To justifyJSlTD note that I + [C]c I + [H] = I + [H], and conversely I+[H] - H C 
J + [C] — 7i, for if p G — "H then there exists a timelike curve from a point 

q G Hr\I~(p) to p which, when extended to an inextendible timelike curve necessarily 
crosses C, and does so to the past of p since any point on the curve between q and 
the crossing point lies in Ti. It follows that 

I + [C] -H = I + [H] -H. (348) 

Now 

i + [H] -h = [i + [H] -n\ -n = [i + [c] -n\ -n = i + [c] - n. (349) 
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In the rightmost expression Ti may be replaced by C: Ti G D[C], hence, by |HE73j 
prop. 6.5.1, for any point r ETi there is a timelike curve passing through both r and 
a point s £ C. Since I + [C] is achronal the curve can intersect I + [C] D C only at the 
one point s, so if r £ -^ + [C] as well as W then r = s £ C*. This establishes (I347p . and 
also the claim of the proposition. I 

O 

Since Af is compact, and thus closed in spacetime, Af= dAf, so prop. IB. 121 implies 

O — 

that for any Cauchy surface C of intD[Af] C— dAf. Indeed, we shall see that C is a 

— o 

manifold with boundary that is homeomorphic to Af ' . It follows that dC =C= dAf. 
To demonstrate homeomorphism we first establish a lemma: 
Let us define the minimal component, C°, of a Cauchy surface C of an open 

hyperbolic set Ti, to be C n I[Ti)& 

Proposition B.13 IfC is a Cauchy surface of an open setTi then C—C° = Tl—I[H]. 

Proof. It is evident that C — C° C Ti — I[H]. What must be shown is thus that 
all points of Ti — I[H\ lie in C — C°. Since C° C I[H] it is sufficient to show that 
7i — I[H] C C. The achronality of C implies that C is also achronal. Therefore, since 
C is a Cauchy surface for H, C is as well. As a consequence 7i C D[C] implying 
that Ti. C -D[C]. By |HE73] prop. 6.5.1, this implies that all inextendible timelike 
curves through H must pass through C. But if p £ H — I[H] then I{p) is disjoint 
from 7i, and indeed from Ti.DC, since I(p) is open. The only point on a timelike 
curve through p not in I(p) is p itself, so this is the only point at which the curve can 
intersect C: p £ C. I 



Proposition B.14 If C and C are both Cauchy surfaces of the open set Hd, 
then C is homeomorphic to C' . 

Proof. Let v be a smooth vector field on X. (By lemma 8.1.1 of [Wald84j such a field 
always exists on a paracompact, time orientable, smooth, Lorentzian metric manifold. 
See also [HE73j pp. 38-40.) Each integral curve of v is inextendible in X, so if it 
passes through Ti it must cross both C and C exactly once. The integral curves thus 

78 C° is the only part of C which causal curves through H actually cross. (Because Ti is open 
J\H] = I[TC\). Thus C° is precisely that part of C that is necessary for it to function as a Cauchy 
surface for Ti.. A Cauchy surface of Ti. contained in Ti is equal to its minimal component. On the 
other hand dAf is not part of the minimal component of Af. Indeed lemmas IB. 91 and IB. 101 (and 
the time reverse of lemma IB.10|) indicate that for any Cauchy surface C° consists of C minus its 
boundary and any null geodesies without future or past endpoint in intC. In particular, by prop. 
ETTJand lemma Af° = intAf. 
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establish a 1-1 correspondence between C° = C nl[H] and C'° = C C\ I[H]. By prop. 
IB. 131 C — C° = H — I[H] = C' — C'° so this 1-1 correspondence actually extends to 
a 1-1 correspondence between C and C'. Moreover, this correspondence is continu- 
ous and has continuous inverse: The smoothness of v implies that the points on the 
integral curves are continuous functions of the curve parameter and of the starting 
point of the curve. Furthermore the achronality of C and C imposes a Lipschitz 
condition on the curve parameter value of the point p' G C that corresponds to a 
point p G C, implying that this parameter is a continuous function on C. p' is thus 
a continuous function of p. Reversing the roles of C and C one sees that the inverse 
is also continuous I 

The preceding results together establish that 

Proposition B.15 The closure of any Cauchy surface ofV> = intD[Af] is an embed- 
ded compact 3-manifold with boundary dAf that is homeomorphic to Af, and thus to 
the closed unit cube in 1R 3 . 

In an open globally hyperbolic region the wave equation, and similar equations, 
such as the linearized Einstein field equation in transverse gauge ([TBI , have unique 
advanced and retarded Green's function (see |CDD82j p. 523). It is therefore relevant 
when using Green's functions in conjunction with initial data on Af, as we do, that 

Proposition B.16 Af is contained in an open subset V of V which, when regarded 
as a spacetime in its own right, is globally hyperbolic. 

The proof consists in extending Af to a larger achronal hypersurface A and showing 
that Af C intD[A] = V. In fact it will not be shown that A is achronal in V, but 
rather in an open subset V, which will serve as the spacetime in which D[A] is defined 

Proof part I. Definition of V and A: V is an open subset of V containing Af, 
that is, like V, an exponential of an open subset (V) of O' from which any ray of 
O' exits at most once. V thus has all the properties that have been demanded of V, 
and could be taken in place of V in all preceding theorems. In addition the closure 
of V in V will be required to be compact. [This requirement is easily fulfilled. For 
instance, define a Riemannian metric on O', and let V' be a subset of the cylinder 
formed by vectors in O' that are shorter than twice the maximal length of the vectors 
in the compact subset N.} As a consequence any geodesic normal to S' crosses the 
boundary of V before leaving V. 

The construction of A proceeds in two steps. The first step is to extend Af to a 
hypersurface Af just like A" but slightly larger. Like A", Af is swept out by future 
directed null geodesic generators normal to a compact spacelike 2-disk that are trun- 
cated before reaching a caustic. The 2-disk Sq on which the generators of Af originate 
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contains So in its interior, and the surfaces Sl and Sr on which the generators are 
truncated are chosen to lie beyond Sl and Sr, so that M C intN '. Recall that O' 
is the normal bundle of a smoothly embedded, boundaryless, spacelike 2-manifold 
5" D S (which is defined before prop. IB. 1 [) . S will be taken to be a compact disk 



with smooth boundary in 5" that contains So in its interior. Furthermore Sl and Sr 
will be chosen close enough to M so that they (and thus all of Af) lie within V'. Thus 
Af = f[N] where N D N is a subset of V, and prop. IB. 61 implies that Af is achronal. 

Having defined Af we now define A to be I~ [Af] . This is an achronal embedded 
3-manifold, by |HE73j prop. 6.3.1. 

To show that Af C intD [A] we use the following lemma: 

Lemma B.17 Suppose X' is the maximal extension (as a null geodesic) of a generator 
of Af in the spacetime V , then X' R I~[Af] C Af in this spacetime. 



Proof. Suppose p G A' PI I~ [Af] — Af. Then there exists a sequence of points {p n } 
in I~ [Af] in V' approaching p, and for each p n a future directed timelike curve 7 n 
in V' from p n to Af. Now consider the problem in the context of the larger ambient 
spacetime V D V . The 7„ are timelike curves in V that are inextendible in V — Af 
so, by [HE73j lemma 6.2.1, they have a causal limit curve 7 through p which is also 
inextendible in V — A/", and thus either inextendible in V or has an endpoint on A/". 

Now note that the proof of prop. IB.6I applies equally well to Af in V and shows 
that it is achronal in V. This implies that 7 coincides with A, the maximal extension 
as a null geodesic of A' in V, to the future of p: 7 lies in I~ [Af] and thus, since is 
achronal, outside I + [A/] . But following A from a point on A/" to p and then continuing 
along 7 one obtains a causal curve that may be deformed to a timelike one, which 
would imply that 7 enters J + [A/], unless 7 coincides with A after p. On the other 
hand the 7„ are all contained in V, and hence 7 is contained in the compact set V. 
7 must thus have an endpoint at An V which does not lie on A/", contained in the 
open set V'. A contradiction has been reached, implying that A' H I~[Af] — M is in 
fact empty in V. I 

proof of prop. \B.16\ part II: Suppose p G N C intN . By prop. IB .71 any null 
geodesic through p lying in A is (an extension of) a generator of M through p. Such 
a generator has a past endpoint on intSo C intA and leaves M at the boundary of 
M. Now recall that A = I~[Af] C I~[Af]. Lemma [B . 1 71 thus implies that the exten- 
sion A' of a generator as a null geodesic cannot remain in A outside A/". A' PI M thus 
has a future endpoint on the boundary of A/", which is contained in intA. since M is 
compact. By prop. IB. 101 and its time reverse, I + (p) and I~(p) intersect -D[A], so 
p G I + (q) fl I~(r) with q, r G D[A]. But it is easily shown that the whole open set 
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I + (q) fl / (r) must be contained in -D[A], implying that p G intD[A]. I 

V = int.D[A] is thus a globally hyperbolic boundaryless manifold which is a lo- 
cally isometric covering of the original embedding of M in the arbitrary, smooth, 
Lorentzian, spacetime X. V is sufficient spacetime to develop our whole theory in 
terms of the domain of dependence of M. In particular we may define the domain 
of dependence of M in V, observables on this domain of dependence, and a Peierls 
bracket between these, which in turn induces our auxiliary pre-Poisson bracket on 
the initial data. Nevertheless V does not in general contain the maximal Cauchy 
development of data given on N '. Even in a perfectly regular solution to Einstein's 
field equations a caustic may develop in the timelike geodesies normal to S within 
the interior of the domain of dependence of M . It is therefore interesting, though by 
no means necessary for the definition of the phase space or the Poisson structure on 
it, to show that a boundaryless manifold M extending V can always be constructed 
that is globally hyperbolic and contains M and the interior of its maximal Cauchy 
development. 

M is constructed by gluing two manifolds: Choose a Cauchy surface Eo of V lying 
to the future of The chronological past / [E ] of £ in V, denoted A, is one of 
the components of M, which of course contains M . The other component is obtained 
as follows. Choose a spacelike Cauchy surface £ of the interior of the domain of 
dependence of M in A, V = inW[N] A}. (Note that V = intD[Af; V] n J~[E ; V].) 
If the metric on T>q satisfies the field equations then it induces admissible, spacelike, 
initial data on E, which in turn define a maximal Cauchy development T> m . T> m is 
the other component of M. 

T> must be isometric to a subset of T> m . Indeed 

Lemma B.18 There exists an isometry uj : V U cr — > T> m where o"o = E D 
intD\N;v\. Moreover, a' Q = u[o~o] is the border u>[T> ] in T> m . 

Proof. The definition of the maximal Cauchy development implies that there exists 
an isometry uj : T>q — >• T> m . This isometry extends to boundary points of V = u[Dq] 
in D m , for all such boundary points are naturally images of points on the boundary 
of V in V: Suppose p' G T>' , and 7' is a (future or past directed) timelike curve in V 
which is inextendible in V Q but has an endpoint in the direction of increasing curve 
parameter at p'. A timelike curve through T>q that is inextendible in V must pass 
through the Cauchy surface Af of V Q to the past and the Cauchy surface E , of V, to 

79 By |HE73] prop. 6.6.8. there exists a continuous "time" function t on V which increases strictly 
to the future along any causal curve, and ranges from —00 to 00. Since Af is compact t has a finite 
maximum on N '. Any level set of t with t greater than this maximum will serve as the desired 
Cauchy surface. 
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the future, it cannot remain in T>q. Thus 7 = cj -1 ^'] must either have an endpoint 
in the sense of increasing parameter value inside T>o (but it cannot because 7' does 
not end inside T>' ), or at a point p on the boundary of T>q in V. 

Now choose a convex normal neighbourhood U p of p in V and a convex normal 
neighbourhood U p > of p' in T> m . 7 is contained U p for sufficiently large curve parameter 
and 7' = a; [7] is contained in for sufficiently large parameter values. Thus Y' = 
U p i PI fl X>o] contains 7' beyond a minimum parameter value. Let r' G 7'n 
K' and construct Riemann normal coordinates on f/ p / based at r'. Then define 
y to be the chart on U p obtained by first pulling back y', restricted to a convex 
normal neighbourhood of r' contained in Y', to U p via uj, and then extending the 
resulting Riemann normal coordinates about r = u;~ 1 (r / ) to all of U p . This defines a 
diffeomorphism <fi = y'^ 1 o y from part of U p to part of U p >. 

The boundary of V in V consists of the achronal hypersurface V — I~T>o to the 
future and the achronal hypersurface V — I + V = j\f to the past. Define 1° to be 
I + (r; U p ) if p lies on the future boundary, and I~(r; U p ) if p lies on the past boundary. 
Since the geodesies from r to a point in J°nDo are entirely contained in T>q <p coincides 
with uj on 1° H TV j/07 and y' o 7' thus define the same curve (for sufficiently large 
parameter value) in H 4 , and, since the limit points p and p' lie in the domains of y 
and y' respectively, y(p) = y'{p')- p therefore lies in the domain of <j) and 4>(p) = p'. 
Indeed a whole neighbourhood W p of p lies in the domain of <p, which is open. Since 
p G 1° this neighbourhood may be taken to lie inside 1°. On W p 4> is then a smooth 
extension of uj. Moreover, since the metric is smooth in both of the charts y and y' 
it is an isometry on the closure of Dq in W p . 

Covering V' with sets of the form W p ' = uj[W p ] a unique extension of uj^ 1 to V is 
obtained by patching together the maps = y~ l oy'. The extension is well defined 
and unique because any two diffeomorphisms extending uj~ x must agree on t>' in the 
overlap of their domains, since they agree on T>' Q in the overlap of their domains. The 
extension is an isometry, and its inverse, which extends u to part of T>q, will be called 
u as well. 

What part of V lies in uj^ 1 [D / q ]7 A boundary point p' of V in V m that lies to the 
past of £' = cj[E] must be the image of a point p G jV. (By prop. IB. Ill ) Through 
p' there must therefore pass a null geodesic which is the isometric image of a portion 
of a generator through p. This null geodesic would be future inextendible but would 
not cross since the generator crosses E only at its future endpoint on dj\f, which 
is not contained in E C T> , by prop. IB. 91 However E' is a Cauchy surface of V m . 
Thus there can be no boundary point to the past of E'. 

A boundary point p' to the future of E' must be the image of a point p lying either 
on the future Cauchy horizon i7 + [E;.4] of E in A, or on the future boundary, So, of 
A itself. If p lies on ff + [E; V] D if + [E;^4] there would be a past inextendible null 
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geodesic from p' which does not cross E', for through p there must exist, by theorem 
8.3.5. of |Wald84j . a null geodesic contained in if + [E;V] C J + [E;V] which is past 
inextendible or has past endpoint on the causal edge of £@ But the causal edge of 
E cannot be included in E since £', being a Cauchy surface contained in J) m , has 
empty causal edge (by |Wald84j . corollary to prop. 8.3.6.). p therefore does not lie 
mH+[E;V}. 

A future directed timelike curve from T>q cannot cross Eo without crossing H + [T>; V] 
- and thus leaving intD[Sf; V] - first, except by passing through <Jo = Y>QC\intD[Sf] V]. 
cr = u)\<Jq\ thus contains the boundary of T>' in T> m . Moreover, all points of o~q are in 
fact mapped into V m : Consider using a different Cauchy surface, E , of V, lying to the 
future of E , in place of E in our construction. The set V = intD[Sf; V] n/~[£ ; V], 
which contains T>$ and in particular <7o, would then be mapped by an isometry u into 
T> m . The uj may be chosen so that it agrees with u on the overlap of their domains, 
including ovl^l uj thus maps all of a into V m . I 

M is the union of A and D m with points in P an d T> m C\I~ [E ] identified according 
to the isometry uj. The open sets of M are those sets which have open intersections 
with A and T> m in the topologies of those manifolds. (This is the standard "gluing 
topology"). 

Proposition B.19 M is Hausdorff 

Proof. Suppose M is not Hausdorff. Then there must exist points p and q in M 
such that all neighbourhoods of q intersect all neighbourhoods of p. If both p and 
q lie in T> m , or both lie in A this cannot occur because T> m and A are Hausdorff. 
Suppose therefore that q G V m and p G A. V m is thus a neighbourhood of q and A 
a neighbourhood of p, so in order that the pair violate the Hausdorff condition it is 
necessary that q G A = T>' = o~' Q (by lemma IB. 181) in T> m and p G T> m = T>q in A. To 



show that M is actually Hausdorff it is thus sufficient to find disjoint neighbourhoods 
of q and p in M for such p and q. 

Let go = ^ _1 (<?)- By lemma IB. 181 qn G cr C intD[J\f;V] which is necessarily 

distinct from p lying in the disjoint set T>q C D[M\V}. Since V is Hausdorff there 
exist disjoint open subsets of V U qo 3 qo and U p 3 p. Indeed we can, and will, require 
that U p C A and U qo C intD[hf] V] fl J~[E , V], where, as in the proof of lemma IE. 181 
E is a Cauchy surface of V lying to the future of E . 



80 The causal edge is defined in |HE73j p. 202. In V, which is globally hyperbolic, and thus strongly 
causal, ec?ge[E] — dAf. Sec IIIT.il p. 196-197 on Alcxandrov topology. 
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uj o lo^ 1 defines an isometry of T>' to another open subset of T> m . Since T>' contains £' the 
unicity of the maximal Cauchy development implies that this isometry extends to an isometry of 
T> m to T> m . If this isometry is non-trivial one may simply compose its inverse with uj to obtain a 
new a;, satisfying all requirements, which agrees with uj on T>q. 
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The desired neighbourhoods of p and q in M are constructed from U qo and U p . An 
open set in M is a pair of open sets, one in A, one in T> m , such that identified pairs 
of points are either both in the sets or both outside. Let u be the isometry of lemma 
IB. 181 corresponding to the case in which £ is replaced by X , and recall from the 
proof of that lemma that uj may be chosen so that it agrees with u on the intersection 
of their domains. With this choice of Co &[U qo ] is an open neighbourhood of q in D m , 
and (U qo C\Vo, u[U go \) and (U p ,u[U p r\Vo\) are disjoint open neighbourhoods of p and 
q are respectively. I 

The paracompactness and time orientability of M follows from that of A and T> m 
and the fact that V m fl A has only one connected component. (See |HE73] p. 14.) 

Proposition B.20 A (and thus M) is achronal in M , and M is globally hyperbolic. 

Proof. Suppose p £ A and 7 is a future directed timelike curve from p. If 7 remains 
in A C V it cannot cross A again because A is achronal in V. If 7 leaves A without 
leaving M it must cross the boundary, £ , of A fl T> m in T> m . If 7 is then to return 
to A it must again cross this same boundary, this time entering A. But 7, which is 
future directed, cannot cross the boundary in both directions since So is a level set 
of a time function which increases along any future directed causal curve. Thus 7 
cannot cross A C A a second time. This establishes the achronality of A in M. 

To prove global hyperbolicity it is thus sufficient to show that all of M lies in 
the domain of dependence of A. Suppose 7 is an inextendible causal curve which 
passes through A. 7 may be extended to an inextendible causal curve, 7', in V, and 
7' must pass through A and through So- Since So is achronal 7' cannot reenter A 
once it leaves, so 7' fl A = 7, proving that 7 crosses A. If on the other hand 7 is an 
inextendible causal curve passing through T> m then it must cross the Cauchy surface 
X of D m , which lies in A, implying that 7 also passes through A and thus through 
A. All inextendible causal curves in M cross A. I 

The demonstration that only those points in M that lie on the same generator are 
causally related (prop. IB. 71) goes through in M precisely as in V. 

Since M (and V) are globally hyperbolic the transverse gauge linearized Einstein 
field equations (JT3l has unique advanced and retarded Green's functions G ± (p,q) 
supported (as distributions) at points q £ J \v)- ( see |CDD82j p. 523) 

Using the global hyperbolicity of the ambient spacetime M or V it may be shown 
that sources with compact support in the interior of -D[Af] do not perturb (to linear 
order) the metric in a spacetime neighbourhood of dM: 

Proposition B.21 If s is a compact subset of the interior of the domain of depen- 
dence of M in M, then J [s] is closed and disjoint from an open neighbourhood of 
dM. 
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Proof. The fact that J ± [s] is closed follows from [HE73] prop. 6.6.1. To prove that 
J^fs] is disjoint from an open neighbourhood of dM it is therefore sufficient to show 
that it is disjoint from dM itself. Now by prop. IB. 151 prop. IA.21 and lemma IB. 91 
I[dM] is disjoint from D[J\f], and thus /[OA/] D J[dM] is disjoint from intD[M\. It 
follows at once that J[s] D / ± [s] is disjoint from dM. I 

All Cauchy developments of M seem to be diffeomorphic to each other as space- 
times (though of course not necessarily isometric). In the following a slightly weaker 
claim, that they are diffeomorphic everywhere except possibly on their future Cauchy 
horizons (i.e. their future boundaries, see |HE73j ). will be proved. This result makes 
possible the concrete representation we have used of the set of Cauchy developments 
as a space of solutions S to the field equations on a common, fixed manifold P (The 
union of the interior of the hypersurface Mq — {t — \x\, \x\ < l,y 2 + z 2 < 1} and 
the interior of its domain of dependence in 1R 4 equipped with the Minkowski metric 
ds 2 = -dt 2 + dx 2 + dy 2 + dz 2 - see Fig. [2]). 

This representation makes the set of Cauchy developments easier to think about, 
but it is actually not necessary for our main results, which concern the Peierls bracket, 
the symplectic 2-form, and the auxiliary pre-Poisson bracket on the initial data. These 
objects are defined entirely in terms of first order perturbation theory about a given 
solution, and thus in the context of a single, spacetime manifold, fixed once this 
solution is chosen. The use of a fixed spacetime manifold in developing the main 
results of the present work is thus justified quite independently of the diffeomorphism 
between Cauchy developments. 

Suppose M (defined as on p. [5]) is smooth and achronal in a smooth, Lorentzian 
metric 4-manifold. Let us define V = intD [A/] U intM ' . Then 

Proposition B.22 V is the domain of dependence of M minus its future Cauchy 
horizon 

Proof. D[M] = D + [Af] follows from prop. IB. Ill and prop. IB. Ill and prop. IB. 71 to- 
gether imply that intN = r[D[M]]nM. Thus the identities inW[M\ = I~[D[N]] n 
I + [D[N]] and D + [J\f] = D + [M] n / + [A/] U M imply V = D[M]n I~[D[J\f]}, or equiv- 
alently V = D[Af] - H+[M] where H + [N] = D[N] - I~[D[M]] is the future Cauchy 
horizon. I 

Proposition B.23 For any two Cauchy developments, 1 and 2, of N ', X>i and T>_2 

are diffeomorphic as manifolds (though not necessarily isometric). 

Proof. To prove diffeomorphism a smooth chart $ : V — > H 4 is constructed for 
each Cauchy development, such that the range, 3>[P], of the coordinates on V is the 
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same for all Cauchy developments. Then, if $1 and $2 are the charts for Cauchy 
developments 1 and 2 respectively, ^J" 1 o $ 2 provides the desired diffeomorphism. 

The construction of the chart <3> proceeds in two steps. First the chart $ is defined 
on a neighbourhood of So and on M itself. Then this chart is extended to all of V. The 
exponential map / from the normal bundle O' over S' to to the ambient spacetime 
X is a smooth immersion except at degenerate points of df. 

Since df is non-degenerate on S' D So, the generalized inverse function theorem 
( |AR67j . p. 19) implies that / maps an open neighbourhood U of 5" in O' diffeo- 
morphically to a neighbourhood U = f[U] of S' in X. Moreover, M = f[N] where 
N C O' is free of degenerate points of df. Since M is achronal / must be one to one 
on N and thus a diffeomorphism on this set as well. 

Coordinates (x°, x 1 , x 2 , x 3 ) are chosen on the normal bundle O' as follows: (x 2 , x 3 ) 
coordinatize the base manifold S', and (x°, x 1 ) the fibres, consisting of the nor- 
mal vectors at each point of S'. (x°,x l ) are linear, that is, they are the compo- 
nents of the normal vectors with respect to some basis defined at each point of S', 
and ( X ■ X - X - X 1 chosen so that intN corresponds to the coordinate domain 

{(x ,^ 1 ,^ 2 ,^ 3 )!!^ 1 ! = x°, (x 1 ! < 1, (x 2 ) 2 + (x 3 ) 2 < 1}. 

$ maps each point p G U U M to the coordinates ( ) of its preimage 

P = f\u\p) inUUN. 

Note that the definition of the coordinates implies that the vector field 8q on So 
is timelike and future directed. 

We turn now to the extension of Since the continuous vector field do is timelike 
and future directed on So it is also so in an open neighbourhood Bo of So in IA. Let 
C be a compact subset of Bo containing So in its interior, then $ will be set equal to 
$' on a compact subset C of Bo containing So in its interior, and on M itself . To 
extend the chart beyond this set do will first be extended to a smooth timelike vector 
field on all of V by splicing it smoothly with a smooth, non-zero, future directed, 
future complete, timelike vector field on V. Here "future complete" means that the 
supremum of the curve parameter in V is infinite on every integral curve of the vector 
field. 

A field u satisfying these requirements can be found using a suitable Riemannian 
metric on V: Choose a partition of unity {f n } of V subordinate to a locally finite 
atlas {<f) n } such that the domain of each chart <p n has compact closure, and define the 
Riemannian metric r = J2n fn<j>*n£, the sum of the pullbacks via the various charts 
of the standard euclidean metric e on 1R 4 weighted by the partition of unity. Let p = 
J2 n n fn- It is easy to show that the Riemannian metric r' = r(| | Vp| | 2 + 1) is complete 
in the sense that any curve without endpoint in V has infinite length, whereas a curve 
from a starting point to an end point, including starting and endpoints on A/", has 
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finite length@ u is taken at each point to be the future directed timelike vector of 
maximal proper time among the unit vectors of r'. The smoothness of r and the 
Lorentzian spacetime metric guarantee that this u is smooth. 

The splicing is defined using a partition of unity (fo, cpi subordinate to the cover 
Bq, B\ = U — C. That is, <po is a smooth function taking values between and 1 such 
that tpo = 1 on C and outside Bq, and ipi = 1 — ipo. w = (fodo + fiu is the spliced 
field, defined on all of V (for this it is sufficient that do is defined on the support (fo) . 
Since both vector fields being added are timelike and future directed, so is their sum 
w. 

Because w is timelike and non-zero on V the integral curves of w are timelike and 
can have no endpoint in infD = V — M . Every integral curve must therefore cross 
the Cauchy surface Af, precisely once. Indeed lemma [BT91 prop. IB.lll and lemma IB~9l 
show that the set of points of Af at which these integral curves cross is precisely intAf '. 

At a point p G V the coordinate x° of the extended chart $ is the parameter of 
the integral curve of w through p, with zero set so that it agrees with the x° of the 
chart $ on Af. The coordinates x l ,x 2 ,x 3 assigned to p by $ are just the a; 1 ,^ 2 ,^ 3 
assigned by to the point at which the integral curve through p crosses Af. (That 
is, x 1 ,^ 2 ,^ 3 are held constant along the integral curves.) 

The smoothness of w, of the submanifold M — S$ and of the chart on U and on 
M — So, i mply , via standard theorems on the solutions of ODEs, that $ is a smooth 
chart of X>@ By prop. IB. Ill D[N] = D + [Af] so the value of x° on Af ', equal to 
\x l \ is the lower bound of the range of x° on a given integral curve, while the future 
completeness of u implies that the upper bound is oo. 

The range of $ on V is therefore 

= {(a; ,^ 1 ,^ 2 ,^ 3 )!^ 1 ! <x° < oo, (x 1 ! < 1, (x 2 ) 2 + (x 3 ) 2 < 1}, (350) 

which is independent of the particular Cauchy development in question. I 

82 To show this note that the length of a curve exceeds the absolute value of the integral of Vp 
along the curve. The claim then follows from the fact that p diverges along any curve without 
endpoint. 

83 The definition of w, and thus $ extends to an open neighbourhood of V_ in the ambient 
spacetime. The local existence and uniqueness theorem 21.1 of |AR67j implies that $ is one to one. 
Theorem 21.2 (flow box theorem) of [AR67] shows that $ is smooth and d$ is invcrtible. Thus, by 
|CDD82j p. 93 theorem 1, $ is a diffeomorphism of an open set containing P into K 4 . 
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C Solutions to linearized gravity in terms of initial 
data and the relation between the variations of 
observables and the presymplectic 2-form 

The key result of subsection 12.2.11 equation ( 1331) . can be proved in a less abstract 
manner, as a corollary to a Kirchhoff type formula for the solutions of the linearized 
Einstein equations in terms of initial data. This formula, demonstrated in the sequel, 
provides us with some analytic information about the dependence of the metric, and 
the observables, on initial data. 

We begin by examining (133|) . As in subsection 12.2. II we chose a spacetime region 
Q bounded by Af and a Cauchy surface S + of the future domain of dependence of 
A/", so that the domain of sensitivity of A lies in the interior of Q. Then 

5A= f oTSg^yf^gd^x. (351) 

The perturbation A^g of the metric induced by A is determined only modulo a 
diffeomorphism generator, but if we require it to obey the transverse gauge condition 
then, by (jSJ) and (jHJ), it may be written as 

&A9*Ay) = -16ttG / a^(x)A lluaT (x,y)^d^x. (352) 

J Q 

(1331 thus holds for all observables of T> — intV if, for all solutions to the linearized 
vacuum Einstein equations Sg, 

5g^(x) = -16ttG1V[A^(x), 5} + 2V ( ^) (353) 

at all points x of £>, for some vector field £. Here A Atl/ (rr) is the tensor field (depending 
on y) obtained from the two point tensor A^ UIJT (x, y) when x, fi, and v are held fixed. 
The £ term reflects the freedom to change 5g by a linearized diffeomorphism. When 
both sides of (1353 j) are integrated against a to recover (1331 it disappears because 
a is divergenceless and vanishes at the boundary of integration. (Note that 



implies that (l3"3l) holds even when A^ is not in transverse gauge, because, as (I8T1) will 
show, when 5 is a solution to the linearized field equations that vanishes near dM 
5\ = for any diffeomorphism generator 
Equation ( 13531) is in fact valid: Suppose 71 = 5\g is a solution to the linearized 
field equations that satisfies the transverse gauge condition, Xu = ^^l^u = with 
7;uv = lia> ~ ^9nu9 ap la P , and let 72 = S 2 g be the (unique) advanced solution to the 
linear, diagonal, second order hyperbolic system 

V a V a %„ - 2R%/% p = 0^, (354) 
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with 9 a smooth, symmetric tensor field of compact support contained in T>. (This 
equation is just (TT3]) with the source 32nGa replaced by 9.) 72 may be expressed in 
terms of the advanced Green's function as 

= -\ f Q e^(x)G- vaT (x } y)V^d A x. (355) 

Here Q is defined as earlier, but with future boundary chosen so that Q contains the 
support of 9. Note that G~ is regular in the sense that when it is integrated against 
a smooth source 9 the result is a smooth solution |CDD82j . 

Now, for any variation 5 the variation of the Einstein tensor is 

SG^ = V [fM Xu) - l/2g^\/ aX a ~ 1/2(V CT V CT V " 2R%» P % P ), (356) 

so 72 would be a solution to the linearized Einstein equations, with source —^9, if 
X2u were zero. 

However, (13541) and the Bianchi identity imply that 

V%X2„ = V%,, (357) 

implying that \ 2 does not vanish unless 9 is divergenceless - which is not assumed 
here. (13571) also has unique, regular retarded and advanced Green's functions H ± . 
Hence 

X2M = f V^ v (x)H; T (x,y)V=9d 4 x (358) 

JQ 

= - / 9^(x)V tl H- T (x,y)V^gd i x (359) 
Now consider the identity used in (|37|) to prove the hypersurface independence of 



Q: 



= S 1 5 2 I-S 2 5 1 I- pi,«f 3 ]J (360) 

= 5i#2]-MN-#*i^2]] (361) 

+T7T7^ / [hg^iG 1 * ~ 5 ig ^6 2 G^}^g~d 4 x (362) 

-iv^a] - / ^[v^ - ^v ffX 5 - V]v^ 4 z. (363) 



In (I363p the facts that 8iG^ u = 0, that 5 2 g vanishes in a neighbourhood of the future 
boundary of Q, and that the orientation of jV is opposite that of dQ have been taken 
into account. It is easy to obtain (1331) directly from this last expression, by setting 
9 = 32nGa, but it is illuminating to press on to obtain an expression for S\g in terms 
of data on J\f . 
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The last integral of ( 13631) can be written as 

/ 7i^VW=^- J / e^5g, uV ^d 4 x. (364) 
JQ 2 Jq 

Since 71 satisfies the transverse gauge condition, integrating the first term by parts 
reduces it to a boundary integral, in fact an integral over M only: 

- / lT r X2rV ZI g^ = I 9^(x)V^^d 4 x, (365) 
JN JQ 

where 

U%)= [ ^ T H- T (x,y)V^gd^ a . (366) 

JAf 

Now 9 is an arbitrary smooth symmetric tensor field with support contained in T> so 
( 13631) . (13641) . and (I366j) imply that at any point p of V 

5 l9liU (p) = levrGfVlG^A] + 2V ( ^), (367) 

where G~ y {p) is the tensor field obtained by fixing the first argument point in the 
advanced Green's function to be p (and the pair of indices associated with the first 
point to be fi and v). (I353j) follows immediately if one notes that the retarded Green's 
function G + (p, q) vanishes for q G J\f since Af lies outside J + (p)- 

Solutions to the linearized field equations that are not in transverse gauge only 
differ from transverse gauge solutions by a linearized diffeomorphism. Such solutions 
can thus also be expressed in terms of initial data via ( 136 7h or ( 13531) . by using £ other 
than that defined by ( 13661) . 



D G rr — 

In this appendix it is shown that on Mr the Einstein field equation G rr = 0, in the 
adapted coordinates a a = (u, r, y 1 , y 2 ), is equivalent to 

T r rr = - r -d r e rj d r e l >. (368) 

o 

We proceed by direct calculation of G rr . The form fl 1 3 [) that the metric takes in 
the adapted coordinates implies that in these coordinates the Christoffel symbols, 

r ^ 7 = ^9 at {^9^ + d p g n - d e g^}, (369) 
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(373) 
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(375) 

Using these facts, and the adapted coordinate form (11Q3[) of the metric, G„ can be 
simplified as follows: 

G rr = R rr — d a Y rr d r Y ar + r^ r ,r Q( g r^r^ (376) 

— ft r u — f) y u — f) Y i -4- r r Y a -T" r" — r r r r e\77\ 

u rr Ur ur u r>- ; r T i r r L ar L ur ur L rr L rr V 01 1 I 

-Y) r Yl (378) 

2 

\&u9 urr 9 9 U Y rrr [&r9 \^rur 9 9 r Y rU r ~\~ (379) 

+Y r rr Y u ur + h r rr - Y u UT Y u ur - ^ - ^e lk e' l d r e, k d r ea (380) 

= -9 ur [du9rr]T r rr + g ur [d r g ru }Y" ur + Y r rr Y u ur - Y u ur Y u ur (381) 

+-r; r + \d r e i3 d r e ij (382) 
r 4 

= *Y r rr + ~d reij d r e ij , (383) 
which implies the result. 
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